Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 114,510 A→T 5.3% intergenic (+23/+12) guaC → / ← hofC GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913114,5100AT 5.3% 331.1 / 6.4 114intergenic (+23/+12)guaC/hofCGMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base A (62/46);  new base T (4/2);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 118 total aligned reads displayed.

GATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAT  >  NC_000913/114362‑114659
                                                                                                                                                    |                                                                                                                                                     
        gggCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAg                                                                                                                                              <  1:867049/150‑1 (MQ=255)
         ggcggcCTTCCTTCAGCTTGTACATACGTTTGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGc                                                                                                                                             <  2:53753/150‑1 (MQ=255)
          gcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCAGCg                                                                                                                                            <  1:808544/150‑1 (MQ=255)
          gcggcCTGCGTTCAGCTTGTACAGACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCg                                                                                                                                            <  1:1327461/150‑1 (MQ=255)
            ggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGCCCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGAAACACGCCAGCGtt                                                                                                                                          >  1:516898/1‑150 (MQ=255)
            ggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCAGCGtt                                                                                                                                          <  2:592459/150‑1 (MQ=255)
            ggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGtt                                                                                                                                          <  2:936823/150‑1 (MQ=255)
             gcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGttt                                                                                                                                         >  2:270981/1‑149 (MQ=255)
             gcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTa                                                                                                                                         <  1:581148/150‑1 (MQ=255)
               cTGCGGTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                                                                                                                                       <  1:1148513/150‑1 (MQ=255)
                tGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcc                                                                                                                                      <  1:1064481/150‑1 (MQ=255)
                tGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcc                                                                                                                                      <  2:1655977/150‑1 (MQ=255)
                  cGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGAGATatccca                                                                                                                                    >  1:1438330/1‑150 (MQ=255)
                   gTTCAGCTTGTACATACGTTGGGGCTTCCCGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACACCCTGAAATCTCCCAACGCTGGCGTGGAGCAACACCCCAGCGTAaacccaa                                                                                                                                   >  1:361348/1‑149 (MQ=255)
                   gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccat                                                                                                                                   >  2:657763/1‑150 (MQ=255)
                      cAGCTTGTACATACGTTGGGTCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatccc                                                                                                                                <  2:173918/150‑1 (MQ=255)
                            gTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACCCGCCCGCGTTATCCCATCCCACTcct                                                                                                                          >  1:1147372/1‑148 (MQ=255)
                             tACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc                                                                                                                         >  1:965802/1‑150 (MQ=255)
                             aacacacCTTGGTGCTTCACGCCGGAAAGAGCTGACCAAGCGCACCCCGTTTCTACGTGTGCAGGAACAAGAACACCGCATCATCAACAACCTGTAATCTCCCCCCGCTCGCGTGGAGCAACACGCCAGCGTAATCCCATCCCACTcatc                                                                                                                         <  2:337308/149‑1 (MQ=255)
                               cATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgc                                                                                                                       >  1:696493/1‑150 (MQ=255)
                                aTACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgca                                                                                                                      <  1:1613355/150‑1 (MQ=255)
                                 tACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcat                                                                                                                     >  1:1642066/1‑150 (MQ=255)
                                       ggggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTa                                                                                                               <  1:1740275/150‑1 (MQ=255)
                                           cTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATg                                                                                                           >  2:1622479/1‑150 (MQ=255)
                                           cTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATg                                                                                                           <  2:709343/150‑1 (MQ=255)
                                            ttCACGCCTGAAAGATCTGACCAAGCGCACCACGTTTATTCGTCTGCAGGAACAACAAAACCGCAGCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATgg                                                                                                          <  1:98589/150‑1 (MQ=255)
                                            ttCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATgg                                                                                                          <  1:725161/150‑1 (MQ=255)
                                                cGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaaaa                                                                                                      >  1:813/1‑150 (MQ=255)
                                                cGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaaaa                                                                                                      <  1:1397270/150‑1 (MQ=255)
                                                 gCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAAt                                                                                                     >  1:1499985/1‑150 (MQ=255)
                                                 gCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAGGCAAAAt                                                                                                     >  1:1410476/1‑150 (MQ=255)
                                                   cTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACACCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTg                                                                                                   >  1:1262204/1‑150 (MQ=255)
                                                   cTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTg                                                                                                   <  2:1445121/150‑1 (MQ=255)
                                                     gAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCCTCGCATCGCCTAAATGGAAAATTGGc                                                                                                 >  1:794138/1‑150 (MQ=255)
                                                     gAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATAGACTAAATGGAAAATTGGc                                                                                                 >  2:727564/1‑150 (MQ=255)
                                                      aaaGAGCTGACCAAGCGCACCACGTTTATTCGTGTTCAGGAACAATAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCa                                                                                                <  2:357851/150‑1 (MQ=255)
                                                         gagCTGACCAAGCGCACCCCGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTGATCCCATCCCCCTCATCGCATCGCCTAAAGGGAAAATTGGCAGAt                                                                                             >  1:1123806/1‑150 (MQ=255)
                                                           gCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACCCGCCAGCGTTATCCCCTCCCACTCATCGCATCGCCTTAATGGAAAATTGGCAGATAc                                                                                           >  2:483228/1‑150 (MQ=255)
                                                           gCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCAAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATAc                                                                                           >  1:959293/1‑150 (MQ=255)
                                                             tGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACAt                                                                                         >  2:1392623/1‑150 (MQ=255)
                                                               aCCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTg                                                                                       >  2:384812/1‑150 (MQ=255)
                                                                ccAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGc                                                                                      <  2:1671635/150‑1 (MQ=255)
                                                                  aaGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGcca                                                                                    >  2:1435396/1‑150 (MQ=255)
                                                                      gcACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccaccac                                                                                >  2:1826728/1‑150 (MQ=255)
                                                                         ccacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAg                                                                             >  2:582618/1‑150 (MQ=255)
                                                                          cacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGACTGGAGCAACACGCCAGCGTTAACCCATCCAACTAATAGCATCGCGTAAATGGAAAATTGGCAGATACATTGCCAACACCAGc                                                                            >  1:1350694/1‑150 (MQ=255)
                                                                             gTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTa                                                                         <  2:1637041/150‑1 (MQ=255)
                                                                               ttATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTAcc                                                                       <  2:545126/150‑1 (MQ=255)
                                                                                 aTTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaa                                                                     <  1:1229476/150‑1 (MQ=255)
                                                                                     gtgtGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGAGTACCaataat                                                                 >  1:191943/1‑150 (MQ=255)
                                                                                        tgCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTATATGGAAAATAGGCAGATACATTGCCACCACCAGCGTACCAATAATtcc                                                              >  1:971392/1‑150 (MQ=255)
                                                                                           aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCCTCCCCCTCATCGCATCGCCTAAATGGAAAATTGGCAGCTACATTGCCACCCCCAGCGTCCCAATAATtcctcc                                                           >  1:834534/1‑150 (MQ=255)
                                                                                            ggAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTccc                                                          >  1:748692/1‑150 (MQ=255)
                                                                                              aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGt                                                        >  2:1568950/1‑150 (MQ=255)
                                                                                               aCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGtt                                                       <  1:1481691/150‑1 (MQ=255)
                                                                                                cAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCCCCACCAGCGTACCAATAATTCCTCCCGttt                                                      >  1:1359759/1‑149 (MQ=255)
                                                                                                cAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTa                                                      >  2:1823358/1‑150 (MQ=255)
                                                                                                   gAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGa                                                   <  1:1564667/150‑1 (MQ=255)
                                                                                                     aaaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACGCATCGCATCTCCTAAATGGAAAATTGGCAGATGCATTGCCCCCACCAGCGTACCAATAATTCCTCCAGTTATGAta                                                 >  2:1229332/1‑149 (MQ=255)
                                                                                                       cccGCAGCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTCTGATCAg                                               <  2:1123111/149‑1 (MQ=255)
                                                                                                       aCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAg                                               >  1:1284603/1‑150 (MQ=255)
                                                                                                       aCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGGGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAg                                               <  1:469142/150‑1 (MQ=255)
                                                                                                        ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCCCGCATCGCATCGCCTAAATGGAAAATTGGCAGATACCTTGCCACCACCCGCGTCCCAATAATTCCTCCCGGTATGATCAGc                                              >  1:888978/1‑150 (MQ=255)
                                                                                                        ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCCACCGCGTACCAATAATTCCTCCCGTTATGATCAGc                                              >  2:1357544/1‑150 (MQ=255)
                                                                                                              cTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGcc                                        >  1:16130/1‑150 (MQ=255)
                                                                                                              cTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCAATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCCTACCAATAATACCTCCCGTTATGATCAGCAACGcc                                        >  2:1733471/1‑150 (MQ=255)
                                                                                                               ttCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCg                                       >  2:207417/1‑150 (MQ=255)
                                                                                                                tCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCgg                                      >  2:1304388/1‑150 (MQ=255)
                                                                                                                tCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAAAAATTCCTCCCGTTATGATCAGCAACGCCgg                                      >  1:242322/1‑150 (MQ=255)
                                                                                                                 caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                     <  1:657763/150‑1 (MQ=255)
                                                                                                                  aacaacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCTCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGtt                                    <  1:270981/150‑1 (MQ=255)
                                                                                                                    caacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCCCTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTcc                                  >  2:1733057/1‑150 (MQ=255)
                                                                                                                      acCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAg                                <  1:490829/150‑1 (MQ=255)
                                                                                                                      acCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAg                                <  1:1613436/150‑1 (MQ=255)
                                                                                                                        cTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCATCGCCGGGTCCAGTa                              >  1:200823/1‑150 (MQ=255)
                                                                                                                          gTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATAATCAGCAACGCCGGTTCCAGTAAg                            <  2:227087/150‑1 (MQ=255)
                                                                                                                           tAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAgg                           <  1:509132/150‑1 (MQ=255)
                                                                                                                                 cccAACGCTGGCGTGGAGCAACACGCCATCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCTGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGa                     <  2:814775/150‑1 (MQ=255)
                                                                                                                                 cccAACGCTGGCGTGGAGCAACACGCCAGCGTAATCCCATCCCACTCATCGCATAGCCTAAATGGAAAATTGGCAGATAAATTCCCAACACCAGCGTACCAATAATTACTCCCGTAAAGATCAGCAACGACGTTTCAAGTAAGGCTGCGa                     >  2:1785787/1‑150 (MQ=255)
                                                                                                                                 cccAACGCTGGAGTGGAGCAACACGCCAGCGTTATCCCATCCCACTGATCGCATCGCCTATATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGa                     <  2:361348/150‑1 (MQ=255)
                                                                                                                                  ccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAg                    >  2:541149/1‑150 (MQ=255)
                                                                                                                                  ccAACGCTGGCGTGGAGCAACAAGCCAGCGTTATCCCATCCCACTCATCGCATCTCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTAGCAATAATTCCTCCCGTTTTGATCAGCAACGCCGGTTCCCGTAAGGCTGCGAg                    <  1:70750/150‑1 (MQ=255)
                                                                                                                                   cAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCGCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAgg                   >  1:318088/1‑150 (MQ=255)
                                                                                                                                   cAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTTAGGCTGCGAgg                   >  2:1737908/1‑150 (MQ=255)
                                                                                                                                   cAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAgg                   >  1:1574744/1‑150 (MQ=255)
                                                                                                                                    aaCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGt                  <  2:1505107/150‑1 (MQ=255)
                                                                                                                                     aCGCTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCCTCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTACCAGTAAGGCTGCGAGGtt                 >  1:331996/1‑150 (MQ=255)
                                                                                                                                       gCTGTCGTGGAGCAACACGACAGCGTTATCCCATCCCACTCATCGCATCGACTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAt               <  2:516898/150‑1 (MQ=255)
                                                                                                                                       gCTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAt               >  1:425632/1‑150 (MQ=255)
                                                                                                                                       gCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGATAATTGGCAGATACATTGCCACCCCCAGCGTACCAATACTTCCTCCCGTTATGCTCAGCAACGCCGGTTAAAGTAAGGCTGAGAGGTTAt               >  2:1478191/1‑150 (MQ=255)
                                                                                                                                        cTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATc              >  1:1249711/1‑150 (MQ=255)
                                                                                                                                          ggCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCg            >  2:1821538/1‑150 (MQ=255)
                                                                                                                                              tGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAg        >  2:1484147/1‑150 (MQ=255)
                                                                                                                                              tGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAg        >  1:756648/1‑150 (MQ=255)
                                                                                                                                               ggAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAgc       >  1:1789242/1‑150 (MQ=255)
                                                                                                                                                gAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAgcg      >  2:10587/1‑150 (MQ=255)
                                                                                                                                                  gCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCATCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTAATCCGCCCCCGcc    >  1:1157903/1‑150 (MQ=255)
                                                                                                                                                    aaCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAACTGGAAAATTGGCAGCTACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAt  >  2:1181896/1‑150 (MQ=255)
                                                                                                                                                    aaCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAt  >  1:1827094/1‑150 (MQ=255)
                                                                                                                                                    aaCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCATCCCGTTATGATCAGCAACGCCGCTTCCAGTAAGGCAGCGAGGTTATCCGCCAGCGCCAt  >  2:1631551/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAT  >  NC_000913/114362‑114659

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