Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,1910CA28.7% 873.1 / 131.9 470intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (96/232);  new base A (110/22);  total (215/255)
Fisher's exact test for biased strand distribution p-value = 6.42e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 475 total aligned reads displayed.

CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552041‑1552340
                                                                                                                                                      |                                                                                                                                                      
aatatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                        <  2:600498/150‑1 (MQ=14)
 atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                       <  6:99457/151‑1 (MQ=31)
   atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaaa                                                                                                                                                     <  6:506886/151‑1 (MQ=33)
   atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                                                                                                      <  3:578456/150‑1 (MQ=32)
   atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                                                                                                      >  5:424301/1‑150 (MQ=32)
    ttGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                    >  4:366847/1‑151 (MQ=33)
      gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                                                                                                                                   <  4:236055/150‑1 (MQ=33)
        aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAACTGgcg                                                                                                                                                >  7:1120713/1‑151 (MQ=18)
          aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGGTTGTCATCAGTTTCAAATGgcgcg                                                                                                                                              >  4:79528/1‑150 (MQ=18)
          aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAAAAGCGAAAAAGACAAAAAATCTTTTGTCATAAGTTTAAAAAGGCGCt                                                                                                                                              >  2:245916/1‑151 (MQ=11)
            aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTGATCAGTTTCAAATGGCGCTGt                                                                                                                                           >  4:127717/1‑151 (MQ=25)
             aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                                                                                                           >  5:1050038/1‑150 (MQ=35)
               gTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTaa                                                                                                                                         >  6:1186823/1‑150 (MQ=21)
                    cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGCg                                                                                                                                   >  5:1539227/1‑151 (MQ=37)
                     tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGGATGATGCAAGGTTGCGTTTGTCATCAGTTTATAATGGCGCGGTATAAGGCGt                                                                                                                                  >  4:1176699/1‑151 (MQ=11)
                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                               >  4:1511460/1‑151 (MQ=255)
                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                               <  2:491323/151‑1 (MQ=255)
                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                               <  5:129536/151‑1 (MQ=255)
                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                               <  3:1225416/151‑1 (MQ=255)
                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                               <  6:365626/151‑1 (MQ=255)
                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCTATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                                                                                                            <  2:142788/151‑1 (MQ=255)
                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTc                                                                                                                           >  6:185895/1‑150 (MQ=255)
                              gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                         <  4:65075/151‑1 (MQ=255)
                               ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTATAAGGCGTCATTTTCat                                                                                                                         >  5:550924/1‑150 (MQ=255)
                                  gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                     <  3:641636/151‑1 (MQ=255)
                                  gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                     <  3:473268/151‑1 (MQ=255)
                                  gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGACAAAGTAAAAGCCTACATTTTCaaaat                                                                                                                      >  8:1637382/1‑146 (MQ=11)
                                    tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGAGGCATTTTGATATTGt                                                                                                                    >  5:946757/1‑150 (MQ=255)
                                     gcggcgTTAACGCGTTATCCGGCATGAAAACCCTTCAAATCCATTAGATTGCAGTGAACGTGTAGGCCTGATACGCGTAGCGCTTCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                   <  4:1040032/150‑1 (MQ=17)
                                         cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGGGGATCAGGCGATGTGGGGTTAGGTCTAATTTTATTAGGGCGCTGGATATGGCGTCTTTTTCTTTTTGTAGACa                                                                                                               >  4:339113/1‑150 (MQ=2)
                                         cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCATCAGTTTAAAAAGAACCTGTAAAAGTCGTCATTTTCAAAATGTaacca                                                                                                               >  6:1455760/1‑146 (MQ=14)
                                          gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGCGGCATTTTCATATTGTAGACAAc                                                                                                             >  3:829302/1‑151 (MQ=255)
                                          gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                                                                                                             >  3:323576/1‑151 (MQ=255)
                                           tGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACg                                                                                                            >  7:1062938/1‑151 (MQ=255)
                                            gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAATGGGTCATTTTAATATTGTAGACAACGt                                                                                                           >  2:1667815/1‑151 (MQ=14)
                                             aaCGCCTTATCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGt                                                                                                           <  2:208485/150‑1 (MQ=255)
                                              aCGCCTTATCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                                                                                                         <  1:270195/151‑1 (MQ=255)
                                              aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAAGAGACATTTTAATATTGTATAAAACga                                                                                                           >  8:1026828/1‑148 (MQ=11)
                                                gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGCGTCATTTTATTATTGTATACAACGTAGGc                                                                                                       >  5:785109/1‑151 (MQ=255)
                                                 gcTTATCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTATCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCt                                                                                                      <  5:1652214/150‑1 (MQ=255)
                                                           cATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGAGTCATTTTCATATTGTATACAACGTCGGCTTTGTTCATGc                                                                                            >  6:1277620/1‑151 (MQ=255)
                                                            attaaaaTCATTCAAATTCTATATATTGAATTGAAAGTGTAGTAATTATAATCGTAGCTCATCAGTCAATGTTGCTTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                           <  2:1177062/148‑1 (MQ=18)
                                                            aTGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGAAAAAGTAGGCTTTGTTCATGcc                                                                                           >  2:357837/1‑151 (MQ=18)
                                                              gAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCg                                                                                          <  3:1419356/150‑1 (MQ=255)
                                                                aaaCCCTTCTAATCCAATAGATTGCAGAGACCGTGTAGGCATGATAAGCGTAGCGCAGCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                       <  8:847059/151‑1 (MQ=37)
                                                                  aCCCTTCAAATCCAATATATTTCATTGAAGGTGTAGGCCTGATAAGCGTTGGTGTTCTGTCAATGTTGCGTTTTTGTTCTGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                     <  7:157737/151‑1 (MQ=18)
                                                                  aCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGAGTCATTTTCATATTGAAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                     >  1:1056757/1‑151 (MQ=255)
                                                                       tCAAATCCAATAGATTGCAGTGAACGTGGAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                                  <  2:715514/149‑1 (MQ=255)
                                                                       tCAAATCCAATAGATTGAAGTTAACGTTAAGGCCTGATAAGCGTAGAGCATTAGTCAATCTTGCGTTTGTCATCAGTTTCATATTGCGATGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                                  <  8:725506/149‑1 (MQ=25)
                                                                        cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATCGAGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGTCTTTGTTCATTCCGGATGAGGCGTg                                                                               >  2:909104/1‑151 (MQ=255)
                                                                         aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGTCATTTTCTTATTGTCTCAAACGTAGGCTTTGTACATGCACGATgcgtcttg                                                                               >  5:673260/1‑145 (MQ=17)
                                                                         aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGGCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                              <  7:501494/151‑1 (MQ=255)
                                                                          aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                              >  8:1207890/1‑150 (MQ=255)
                                                                          aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGGGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTAATGCCGGATGCGGCGTGa                                                                              >  5:1380987/1‑150 (MQ=255)
                                                                             ccAATAGATTGCAGTGAACGTTTAGTCATGATAACCGTATCTAATCAGTCAATGTTGATTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                          <  8:62704/151‑1 (MQ=18)
                                                                              cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCCTAAGTTTCAAATTCAGCTGTAAAATGAGAATTTTAGATATTGTATACAAGGTACTCTTTGTTAAATCATCATGCGGCGTGAACGcg                                                                         >  4:952568/1‑150 (MQ=9)
                                                                               aaTAGATTGCACTGAACGTCTAGGCCTGATAAGCCTAGCGCATCAGGCAATGTAGCGTTTCTCATCACTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                        <  5:358157/151‑1 (MQ=21)
                                                                                aTAGATTGCCGTGCACGTGTAGGCCTGATAAGGGTAGGGCATGATGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                       <  1:673390/151‑1 (MQ=21)
                                                                                 tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGATCATAAGAGAATCTTGGGTTTGTAATCAGTTTCAAATAGAGATGAAAAAGGGGGGATGTTGATATTGTAGAAAAAGTAGGCTTTGTTCATGCCGGCTGCGGCGTGAACGCCttg                                                                      >  8:1455200/1‑150 (MQ=16)
                                                                                    aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCGTCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGAATACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAACGCCGTagac                                                                   >  6:819699/1‑148 (MQ=21)
                                                                                    aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAAGAGNCATTTTCATATTGTAGACAAAGTAGGCTTTGTTCATTCACGATCAGGCTTGAACGCATTATcc                                                                   >  1:853266/1‑151 (MQ=17)
                                                                                      tGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                  <  3:913789/150‑1 (MQ=255)
                                                                                       gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                <  8:1122081/151‑1 (MQ=255)
                                                                                          gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                             >  4:477894/1‑151 (MQ=255)
                                                                                              aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCGTTATCCGGCATGaaa                                                          >  7:324303/1‑150 (MQ=255)
                                                                                                gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                       >  1:115796/1‑151 (MQ=255)
                                                                                                 tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATCAAAAccc                                                      >  8:1357194/1‑151 (MQ=255)
                                                                                                  gtAGGCCTGATAAGCGTAGCGCATCAGGCTATGTTTCGTTTGTCATCAGTTTCAAATGTCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                       <  8:180932/149‑1 (MQ=255)
                                                                                                   tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCtt                                                    >  3:240716/1‑151 (MQ=255)
                                                                                                    aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                    <  8:517072/150‑1 (MQ=255)
                                                                                                     ggCCTGATACGCGTAGCGCATCAGGCAATGTTGCGTGTGTCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                  <  5:738237/151‑1 (MQ=255)
                                                                                                       ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGGGTCATTTTGATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGATCGCCTGATCCGGCATCACAACCCTTGaa                                                 >  3:1582071/1‑150 (MQ=21)
                                                                                                        cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAt                                               >  3:403300/1‑151 (MQ=255)
                                                                                                         tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATc                                              >  7:868680/1‑151 (MQ=255)
                                                                                                         tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATCCCGGATGCGGCGTGAACGCCGTATCCGGCATGAAAACCCTTCaact                                               >  7:1449492/1‑148 (MQ=255)
                                                                                                          gATAAGCGTAGCGCATCAGGCAATGTTGCGTCTCTCATCAGATTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                             <  7:676602/151‑1 (MQ=255)
                                                                                                             aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                          <  4:1562758/151‑1 (MQ=255)
                                                                                                              atcGTAGCGCATCAGGCTATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                         <  1:1273387/149‑1 (MQ=255)
                                                                                                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                         <  8:402986/151‑1 (MQ=255)
                                                                                                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGAGGCGTGAACGCATTATCCGGCATGAAAACACTTCAAATCCAAt                                          >  1:501336/1‑150 (MQ=18)
                                                                                                                 gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                          >  4:1443114/1‑147 (MQ=255)
                                                                                                                  tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAtt                                     >  7:1332093/1‑151 (MQ=255)
                                                                                                                   aGCGCATCATGTAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCTTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                    <  7:1179416/151‑1 (MQ=255)
                                                                                                                   aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGAttt                                    >  6:1167440/1‑150 (MQ=255)
                                                                                                                        aTCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                               >  2:1282656/1‑151 (MQ=255)
                                                                                                                          gagtcAATGTTTCGTTTGTCATAAGTTTCAAATGGCGCTGTATAAGGCGTCATTTTCATATTGTAGTCTACGTAGGCTTTGTTCATGCCGGATGGGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGGAGTGat                             <  2:1515592/147‑2 (MQ=14)
                                                                                                                          cAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGaa                             <  6:64832/151‑1 (MQ=255)
                                                                                                                          cAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTAATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGa                              >  6:1594405/1‑150 (MQ=255)
                                                                                                                          cAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAAGAAAACCCTTCAAATCCAATAGATTGCAGTGa                              >  7:607690/1‑150 (MQ=255)
                                                                                                                              cAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGACTTATACGGCATGAAAACCCTTCAAATACAATAGATTGAAGTGAACgtg                         >  1:1639522/1‑151 (MQ=17)
                                                                                                                                  gtttcttttGGCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                      <  3:1039196/144‑1 (MQ=255)
                                                                                                                                  gTTGCGTTTCTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                     <  3:1368409/151‑1 (MQ=255)
                                                                                                                                      gctttGTCATCAGTTTCAAATGGCGCTGTAAAATTCGTTATTTTCATATTGTAGTCTATGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                 <  1:370163/149‑1 (MQ=14)
                                                                                                                                      cGTTTGTTATTAGTTTCAAATTGGCTGTTAATACTCTTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                 <  2:1058216/151‑1 (MQ=11)
                                                                                                                                         ttGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt               >  1:1036972/1‑150 (MQ=255)
                                                                                                                                             cataaGGTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           <  5:234638/150‑1 (MQ=9)
                                                                                                                                              atcaGGGTCAACTGGCGCTGTAACAGGCGGCATTTTCATATTGTATACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAACGCCTTATCGGTCATTAAAAGCCTGCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt         <  5:1454399/151‑1 (MQ=11)
                                                                                                                                               tcaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          <  7:370440/149‑1 (MQ=21)
                                                                                                                                               taaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt         <  6:1029686/148‑1 (MQ=17)
                                                                                                                                                    ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc    >  8:342276/1‑150 (MQ=33)
                                                                                                                                                      |                                                                                                                                                      
CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552041‑1552340

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: