Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,033 T→G 5.4% intergenic (+23/‑57) yraH → / → yraI putative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0330TG 5.4% 1145.8 / 43.8 356intergenic (+23/‑57)yraH/yraIputative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (162/172);  new base G (6/14);  total (170/186)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 361 total aligned reads displayed.

CCCGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACG  >  NC_000913/3287884‑3288182
                                                                                                                                                     |                                                                                                                                                     
 ccGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTg                                                                                                                                                      >  2:4799017/1‑150 (MQ=255)
       gAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGAt                                                                                                                                                >  2:3852850/1‑150 (MQ=255)
       gAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGAt                                                                                                                                                <  1:3797292/150‑1 (MQ=255)
        aaTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCTCCCCAGGGGAAGCAGACGCTATCGTTTTCGTAACGCTGGATTACCGTTTATACGTTAGGGCGTTATCTGACCTGTCAGATa                                                                                                                                               <  1:6171905/150‑1 (MQ=255)
        aaTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATa                                                                                                                                               <  1:4489440/150‑1 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAc                                                                                                                                             >  2:5091295/1‑150 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAc                                                                                                                                             <  2:1467425/150‑1 (MQ=255)
             aaTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGCAATCCCCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGcc                                                                                                                                          <  2:3121989/150‑1 (MQ=255)
                 ttATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCTGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTTGGGCGTTATCTGACCTGTCAGATAACGCCCttt                                                                                                                                      <  1:6730713/149‑1 (MQ=255)
                       tCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTTGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttc                                                                                                                                <  2:1180449/150‑1 (MQ=255)
                       tCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttc                                                                                                                                >  2:1268183/1‑150 (MQ=255)
                       tCAATGTTGACCTCAACGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTatccttc                                                                                                                                >  2:2894083/1‑150 (MQ=255)
                          aTGTTCACCAAATCGATAACAACGTAAAAACCGTATTTCCAACACAAATCACCCCAGGGTTAGCAGACGATATCGTTAACGTAACGCTGTATTACCTTTAATACGTTAGGGCTTTTTCGGACCTGTAATATAACGCCCTTTTCCTTCctc                                                                                                                             <  1:4606961/150‑1 (MQ=255)
                            gTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCtt                                                                                                                           >  1:861674/1‑150 (MQ=255)
                            gTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCtt                                                                                                                           >  2:5863891/1‑150 (MQ=255)
                             ttGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCttt                                                                                                                          >  1:6585982/1‑150 (MQ=255)
                             ttGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCttt                                                                                                                          >  1:4710582/1‑150 (MQ=255)
                                aCCTAATCGCTTACTACTTATCAACCGTTTTTACATCACAAATCACCCCAGGGTAAGCAGACGCTATCGTTAACGTACCGCTGTATTACCGTTAATACTTTTGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtctc                                                                                                                       <  1:1708655/150‑1 (MQ=255)
                                 ccTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCg                                                                                                                      <  2:1633938/150‑1 (MQ=255)
                                   tCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                                                                    >  1:4094086/1‑150 (MQ=255)
                                      tCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                                                                                 >  1:4857714/1‑150 (MQ=255)
                                         cTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCa                                                                                                              >  2:4492463/1‑150 (MQ=255)
                                          ttACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAg                                                                                                             <  2:2472985/150‑1 (MQ=255)
                                          ttACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAg                                                                                                             >  1:1308694/1‑150 (MQ=255)
                                           tactacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGTCCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAgg                                                                                                            >  1:1751350/1‑150 (MQ=255)
                                             ctacGTATCAACCTTAGTTCCATCAGAAATCACCCCCGGGGAAGCAGACGCTATCGTTAACGTACCGCTGGATTACCTTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                          <  1:1834031/150‑1 (MQ=255)
                                             ctacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                              >  1:3190490/1‑114 (MQ=255)
                                             ctacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                              <  2:3190490/114‑1 (MQ=255)
                                                 gTATCAACCTTAGTTCCATCATAAATCACCCCTGGGTAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACTTTTGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaa                                                                                                      <  1:2872870/150‑1 (MQ=255)
                                                  tATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaaa                                                                                                     >  2:6485663/1‑150 (MQ=255)
                                                     cAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATa                                                                                                  <  1:387086/150‑1 (MQ=255)
                                                      aaCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATaa                                                                                                 <  2:6575615/150‑1 (MQ=255)
                                                      aaCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTAGAAAAATaa                                                                                                 <  2:1787961/150‑1 (MQ=255)
                                                          gTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTAt                                                                                             <  1:4189038/150‑1 (MQ=255)
                                                          gTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTAt                                                                                             <  2:234275/150‑1 (MQ=255)
                                                                ccATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGATATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGGGGTATCAGGTTGAAAAATAACTATGTCaaa                                                                                       >  2:6185667/1‑150 (MQ=255)
                                                                 cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCaaaa                                                                                      >  1:4369588/1‑150 (MQ=255)
                                                                  aTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGGTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATATCTATGTCAAAAc                                                                                     <  1:3515711/150‑1 (MQ=255)
                                                                  aTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAAc                                                                                     >  1:6602569/1‑150 (MQ=255)
                                                                  aTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGACCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACAATGTCaaaaa                                                                                     >  2:6255221/1‑149 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACg                                                                                    >  2:4663914/1‑150 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACg                                                                                    >  2:423479/1‑150 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACg                                                                                    >  2:6414914/1‑150 (MQ=255)
                                                                           cACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcg                                                                            >  1:6492572/1‑150 (MQ=255)
                                                                            aCCCCTGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGAGTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgc                                                                           <  1:1860362/150‑1 (MQ=255)
                                                                              cccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCGTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGACCATTCGCgg                                                                         >  1:3972894/1‑150 (MQ=255)
                                                                               ccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                                                                        >  2:4424733/1‑150 (MQ=255)
                                                                               ccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                                                                        >  2:4917006/1‑150 (MQ=255)
                                                                                cAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTg                                                                       >  2:5392923/1‑150 (MQ=255)
                                                                                        gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAgg                                                                                                            >  1:536634/1‑105 (MQ=255)
                                                                                         cAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTGAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATAAGGTTGAAAAATAACTATGTCAAAAAGAACATTAGCGGTGATATTAAcc                                                              >  1:1940460/1‑150 (MQ=255)
                                                                                          aGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCt                                                             >  2:2278067/1‑150 (MQ=255)
                                                                                           gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGTCCTGTCAGATAACGCCCTTTTCCTTCGTCTTTCTGGTTGTTTGTGGTTGTAAAATAACTATGTCAAAACGAACATTGGCGGTGATAttttggtt                                                            >  1:6423127/1‑144 (MQ=255)
                                                                                           gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCtg                                                            >  2:5544093/1‑149 (MQ=255)
                                                                                              gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt                                                         >  2:6718701/1‑150 (MQ=255)
                                                                                              gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATaa                                                                                                                                              <  2:1932574/65‑1 (MQ=255)
                                                                                               cTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGCTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg                                                        <  2:5517057/150‑1 (MQ=255)
                                                                                                    gTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAg                                                   >  2:2811255/1‑150 (MQ=255)
                                                                                                    gTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGGTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAg                                                   >  1:5327807/1‑150 (MQ=255)
                                                                                                    gTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTACTTCCTCTTTCTCGGTGTATCAGGTTGAAAAATAACTATGTCAAAAAGAACATTAGCGGTGATATTAAACTTGTTGTGTAg                                                   >  1:1602766/1‑150 (MQ=255)
                                                                                                        acgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTc                                               <  1:5093133/150‑1 (MQ=255)
                                                                                                         cgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCATTTCCTTCCTCTTTCTCGTTGTATAAGGTTGTAAAATAAATATGTCAAAACGAACATTCGCGGTGATATTAACCGTTTTTTGTAGCttat                                              >  1:580511/1‑148 (MQ=255)
                                                                                                           taacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGt                                            <  1:2817333/150‑1 (MQ=255)
                                                                                                           taacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGt                                            <  1:2333989/150‑1 (MQ=255)
                                                                                                             acgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAt                                          <  1:1387498/150‑1 (MQ=255)
                                                                                                              cgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAtt                                         <  2:3438648/150‑1 (MQ=255)
                                                                                                               gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                                        <  2:3064529/150‑1 (MQ=255)
                                                                                                               gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCGTTCTCGTTGTATCAGGTTGTAAAATAACTATGTCAAAACGAACAATCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                                        >  1:990647/1‑150 (MQ=255)
                                                                                                                   gATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCa                                    <  2:6347955/150‑1 (MQ=255)
                                                                                                                         cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTAt                                                                                             <  1:6027318/87‑1 (MQ=255)
                                                                                                                         cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTAt                                                                                             >  2:6027318/1‑87 (MQ=255)
                                                                                                                         cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCt                              >  2:1150821/1‑150 (MQ=255)
                                                                                                                          gTTAATACGTTAGGGCGTTATCTGAGCTGTCAGATAACGGCCTTTTCCTTCCTCTTTGTCGTTGTATCAGGTTGAAATTTAAGTATGTCAAAACGTTGATTCGCGGTGTTATTAAGGTTGTTGTGTAGGTTGTGTATTGGGGAGGCGGtt                             >  2:3558947/1‑150 (MQ=255)
                                                                                                                           ttAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTg                            <  1:1150821/150‑1 (MQ=255)
                                                                                                                             aaTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCATTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACATTGTTGTGTAGCTTCTGTATTGGCCCGGCGCTTGCa                          >  2:1538081/1‑150 (MQ=255)
                                                                                                                              aTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCag                         <  1:5746931/150‑1 (MQ=255)
                                                                                                                               tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagg                        >  2:5226538/1‑150 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                       <  1:2051851/150‑1 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                       <  2:4120063/150‑1 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                       <  1:6251615/150‑1 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                       <  2:4710582/150‑1 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTACTTCCTCTTTCTCGTTGTATCAGGTTTAAAAATAACTATGTCAAAACGAACATTCGCGTTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                       <  1:5056474/150‑1 (MQ=255)
                                                                                                                                 cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCaggag                      <  2:5299521/150‑1 (MQ=255)
                                                                                                                                 cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCaggag                      <  1:5575806/150‑1 (MQ=255)
                                                                                                                                      gggCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATc                 <  2:4968821/150‑1 (MQ=255)
                                                                                                                                       ggCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCg                >  2:2195257/1‑150 (MQ=255)
                                                                                                                                        gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGACTCGt               >  2:5764549/1‑150 (MQ=255)
                                                                                                                                        gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGt               >  2:475671/1‑150 (MQ=255)
                                                                                                                                        gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGt               >  2:3342628/1‑150 (MQ=255)
                                                                                                                                         cGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtt              >  2:3934050/1‑150 (MQ=255)
                                                                                                                                         cGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtt              >  1:1024690/1‑150 (MQ=255)
                                                                                                                                          gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt             >  1:3344839/1‑150 (MQ=255)
                                                                                                                                          gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt             >  2:3669097/1‑150 (MQ=255)
                                                                                                                                          gTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt             >  1:3108109/1‑150 (MQ=255)
                                                                                                                                            tATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTa           <  1:480060/150‑1 (MQ=255)
                                                                                                                                              tCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACa         >  1:2195445/1‑150 (MQ=255)
                                                                                                                                              tCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACa         <  2:261664/150‑1 (MQ=255)
                                                                                                                                                tGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGc       >  1:2603100/1‑150 (MQ=255)
                                                                                                                                                  aCCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGa     >  2:2717022/1‑150 (MQ=255)
                                                                                                                                                     tGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAATAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAAcg  >  2:6195554/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACG  >  NC_000913/3287884‑3288182

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