Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 1,423,597 T→A pseudogene (1627/2646 nt) yfcV → fimbrial assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,423,5970TA97.7% 1458.3 / 9.6 393pseudogene (1627/2646 nt)yfcVfimbrial assembly protein
Reads supporting (aligned to +/- strand):  ref base T (6/3);  new base A (188/196);  total (194/199)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 412 total aligned reads displayed.

GATGCCAGTGACAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACTACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCGACAACAGTACCG  >  DH1/1423450‑1423744
                                                                                                                                                   |                                                                                                                                                   
   gCCAGTGACAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAacaacaaca                                                                                                                                               <  5:831222/151‑1 (MQ=255)
     cAGTGACAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAAcatc                                                                                                                                             >  1:1016140/1‑151 (MQ=255)
       gTGACAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAAcatcat                                                                                                                                           <  1:476238/151‑1 (MQ=255)
         gACAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATg                                                                                                                                          <  7:752244/150‑1 (MQ=255)
          aCAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGct                                                                                                                                        >  5:711121/1‑151 (MQ=255)
            aGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTAAACAGCCACTTATAAACAGAAATTCCGCGAAGCGTGACTTTCGGTTTATCTCAACTATACCAGCGATACGTACTCGGATCGTGAGGAGCAGACAAACAACAACATCATCctct                                                                                                                                      >  3:374963/1‑151 (MQ=255)
                aaaTGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCa                                                                                                                                  <  1:6752/151‑1 (MQ=255)
                  aTGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACt                                                                                                                                >  5:676946/1‑151 (MQ=255)
                   tGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACt                                                                                                                                >  1:861579/1‑150 (MQ=255)
                    ggTCGGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTAt                                                                                                                              <  4:1003634/151‑1 (MQ=255)
                       ccGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTAttt                                                                                                                            <  1:296533/150‑1 (MQ=255)
                       ccGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTc                                                                                                                           <  8:289644/151‑1 (MQ=255)
                       ccGTACGGGCAACGACAAAGAGATGTACACCGCCACTTAAAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTc                                                                                                                           <  6:1325734/151‑1 (MQ=255)
                        cGTACGGTCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCa                                                                                                                          <  8:1766048/151‑1 (MQ=255)
                           aCGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAat                                                                                                                        <  1:229553/150‑1 (MQ=255)
                            cGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAatat                                                                                                                      <  6:1148079/151‑1 (MQ=255)
                            cGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAatat                                                                                                                      <  7:1473008/151‑1 (MQ=255)
                             gggCAACGACAAAGAGATGTACACCGCCACTTATTACCAGAAGTTCCGCGATGCGGGTGTTTGGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAatat                                                                                                                      <  4:299714/150‑1 (MQ=255)
                              ggCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATgg                                                                                                                    <  5:1078312/151‑1 (MQ=255)
                              ggCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATgg                                                                                                                    >  4:1339367/1‑151 (MQ=255)
                              ggCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATgg                                                                                                                    <  6:271616/151‑1 (MQ=255)
                                 aaCGACAAAGAGATGCACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGgca                                                                                                                 <  4:1668578/151‑1 (MQ=255)
                                   cGACAAAGAGATGTACACCGCTACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGgcagc                                                                                                               >  1:743133/1‑151 (MQ=255)
                                    gACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGgcagca                                                                                                              <  3:4461/151‑1 (MQ=255)
                                       aaaGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTc                                                                                                           <  8:843614/151‑1 (MQ=255)
                                         agagTTGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGt                                                                                                         >  5:1667434/1‑151 (MQ=255)
                                         agagATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGt                                                                                                         >  4:841826/1‑151 (MQ=255)
                                           ataTGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTaa                                                                                                       <  4:1549630/149‑1 (MQ=255)
                                           agaTGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGTTTTATCTCAACTATACCTGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTaa                                                                                                       <  6:592373/151‑1 (MQ=255)
                                             aTGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAat                                                                                                      >  5:560570/1‑150 (MQ=255)
                                               gTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAatat                                                                                                    <  7:842009/150‑1 (MQ=255)
                                               gTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAatat                                                                                                    <  3:1600288/150‑1 (MQ=255)
                                               gTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATg                                                                                                   >  7:1138965/1‑151 (MQ=255)
                                                 acacCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGa                                                                                                  <  4:1130686/150‑1 (MQ=255)
                                                 acacCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAg                                                                                                 <  1:1740024/151‑1 (MQ=255)
                                                 acacCGACACTTATAACTAGATCTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGACATACCTTCTGGGATCGTGAGGATCAGACAAACAACAACATCATGCTGTCCCACTATTTCAATATGGGCAGCATTCGTAATATGa                                                                                                  <  6:1618265/150‑1 (MQ=255)
                                                  cacCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGc                                                                                                <  2:938115/151‑1 (MQ=255)
                                                     cGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGtt                                                                                             >  2:1198879/1‑151 (MQ=255)
                                                         aCTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTccc                                                                                         <  2:1817185/151‑1 (MQ=255)
                                                         aCTTATAACCAGAACTTCCGCCATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTccc                                                                                         >  1:518917/1‑151 (MQ=255)
                                                          cTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCt                                                                                        >  4:403931/1‑151 (MQ=255)
                                                          attATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCt                                                                                        <  5:473537/150‑1 (MQ=255)
                                                            tataACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCt                                                                                        <  3:765914/149‑1 (MQ=255)
                                                            tataACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCt                                                                                        <  7:1728516/149‑1 (MQ=255)
                                                            tataACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTg                                                                                       <  7:463830/150‑1 (MQ=255)
                                                               aaCCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGTGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACt                                                                                    <  4:460204/150‑1 (MQ=255)
                                                                  cAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCt                                                                                >  3:870739/1‑151 (MQ=255)
                                                                    gAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTAc                                                                              >  7:971139/1‑151 (MQ=255)
                                                                       cTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGc                                                                           >  3:1352848/1‑151 (MQ=255)
                                                                       cTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGc                                                                           >  2:1551010/1‑151 (MQ=255)
                                                                           cgcgATGCGGGTGTTTCGGTTTATCTCAACTATACCCGACATACATACTGCGATCGTGAGGAGCAGACAAACAACAACAACATGCTCTCCCACTATTTCAATATTGACAGCATTAGTAATATGAGCGTTTCACTTACTTGCTACCGCTACg                                                                       >  2:1709229/1‑151 (MQ=255)
                                                                                  cGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATg                                                                 >  8:745987/1‑150 (MQ=255)
                                                                                   gggTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGAc                                                               >  2:1365236/1‑151 (MQ=255)
                                                                                    ggTGTTTCGGTTTATATCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGCGCGTTTCCCTGAGTGGCTACCGCTACGAGTATGAc                                                               >  6:443628/1‑150 (MQ=255)
                                                                                      tgtTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAAc                                                            >  6:1489938/1‑151 (MQ=255)
                                                                                       gtTTCGTTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAAc                                                            <  8:1228356/150‑1 (MQ=255)
                                                                                       gtTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAAc                                                            <  7:643407/150‑1 (MQ=255)
                                                                                         ttCGGTTTATCTCAACTATACCCGGCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCgg                                                         >  2:69809/1‑151 (MQ=255)
                                                                                           cGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATAATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGc                                                       >  2:1094994/1‑151 (MQ=255)
                                                                                            ggTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACACCCGGGCg                                                      >  7:1388194/1‑151 (MQ=255)
                                                                                            ggTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCg                                                      <  7:1621219/151‑1 (MQ=255)
                                                                                             gTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCgg                                                     >  3:1052371/1‑151 (MQ=255)
                                                                                              tttATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCgg                                                     >  5:1416125/1‑150 (MQ=255)
                                                                                              tttATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCgg                                                     >  4:1418279/1‑150 (MQ=255)
                                                                                               ttATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGa                                                    <  6:1061994/150‑1 (MQ=255)
                                                                                                 aTCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATa                                                  <  3:517829/150‑1 (MQ=255)
                                                                                                     cAACTATACCCGCCATACGTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAggg                                              <  4:909729/150‑1 (MQ=255)
                                                                                                     cAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAgg                                               <  8:574941/149‑1 (MQ=255)
                                                                                                     cAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAg                                                <  3:1802859/148‑1 (MQ=255)
                                                                                                       aCTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCAt                                           >  3:587645/1‑151 (MQ=255)
                                                                                                        cTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATg                                          >  1:1802778/1‑151 (MQ=255)
                                                                                                        cTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATg                                          >  3:849234/1‑151 (MQ=255)
                                                                                                        atataCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATg                                          <  2:464554/150‑1 (MQ=255)
                                                                                                          ataCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGt                                         >  7:596028/1‑150 (MQ=255)
                                                                                                          ataCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGt                                         >  1:1118352/1‑150 (MQ=255)
                                                                                                          ataCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGt                                         >  6:1287864/1‑150 (MQ=255)
                                                                                                          ataCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTa                                        <  1:960299/151‑1 (MQ=255)
                                                                                                          ataCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTa                                        <  2:862803/151‑1 (MQ=255)
                                                                                                            aCCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGt                                         >  4:1734166/1‑148 (MQ=255)
                                                                                                                gCCATACCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTc                                  <  5:393693/151‑1 (MQ=255)
                                                                                                                     aCCTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCa                             <  3:1220973/151‑1 (MQ=255)
                                                                                                                      ccTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTc                              >  8:1562593/1‑149 (MQ=255)
                                                                                                                      ccTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAg                            <  7:1219483/151‑1 (MQ=255)
                                                                                                                       cTACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGt                           >  1:1279384/1‑151 (MQ=255)
                                                                                                                        tACTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGt                           >  2:1756191/1‑150 (MQ=255)
                                                                                                                         aCTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGt                           <  1:374327/149‑1 (MQ=255)
                                                                                                                         aCTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTAt                         <  3:939965/151‑1 (MQ=255)
                                                                                                                         aCTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTAt                         <  8:1418345/151‑1 (MQ=255)
                                                                                                                         aCTGGGATCGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTAt                         <  5:1511952/151‑1 (MQ=255)
                                                                                                                                 cGTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTggg                 <  7:34916/151‑1 (MQ=255)
                                                                                                                                  gTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTgggg                >  7:1028734/1‑151 (MQ=255)
                                                                                                                                  gTGAGGAGCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGAATGCCGTgggg                >  6:909283/1‑151 (MQ=255)
                                                                                                                                    gaggagCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCg              <  3:1692253/151‑1 (MQ=255)
                                                                                                                                     atgagCAGACAAACAACAACATCAAGCTCTCCCACTATTTCAATATGCGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCGCAGTATGCCGTGGGGCGa             <  4:744522/149‑1 (MQ=255)
                                                                                                                                     aggagCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCg              >  2:1318775/1‑150 (MQ=255)
                                                                                                                                       gagCAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATAAGAGCCTATCCATGACTCGCTACCGCTACTAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGAATGACGTGGGTCGac            >  1:1472798/1‑150 (MQ=255)
                                                                                                                                          cAGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCGacaac         >  5:1830247/1‑150 (MQ=255)
                                                                                                                                           aGACAAACAACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCGacaac         >  1:1341959/1‑149 (MQ=255)
                                                                                                                                             acaaacaaCAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGTCGTGGGGCGACAACAGTa     <  7:245153/151‑1 (MQ=255)
                                                                                                                                                aacaacaacATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCGACAACAGTACCg  >  3:668385/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGCCAGTGACAGCGAAATGGTCCGTACGGGCAACGACAAAGAGATGTACACCGCCACTTATAACCAGAACTTCCGCGATGCGGGTGTTTCGGTTTATCTCAACTATACCCGCCATACCTACTGGGATCGTGAGGAGCAGACAAACTACAACATCATGCTCTCCCACTATTTCAATATGGGCAGCATTCGTAATATGAGCGTTTCCCTGACTGGCTACCGCTACGAGTATGACAACCGGGCGGATAAGGGCATGTAGATTTCGCTCAGTATGCCGTGGGGCGACAACAGTACCG  >  DH1/1423450‑1423744

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: