Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,415,3950GA41.7% 99.7 / 45.7 85R664R (CGC→CGT) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base G (14/35);  new base A (35/0);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 2.48e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTG  >  NC_000913/1415246‑1415545
                                                                                                                                                      |                                                                                                                                                       
gCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGc                                                                                                                                                        <  2:251790/151‑1 (MQ=255)
     aaCGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGtt                                                                                                                                                   <  1:56477/150‑1 (MQ=255)
     aaCGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTa                                                                                                                                                  <  1:384827/151‑1 (MQ=255)
     aaCGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTa                                                                                                                                                  <  1:1109789/151‑1 (MQ=255)
          tttGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTCAGGTGCTACCATAAAAC‑GGGTACt                                                                                                                                                >  2:1030800/1‑148 (MQ=255)
           ttGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTATGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGaa                                                                                                                                             <  1:231055/150‑1 (MQ=255)
            tGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCATGTGCTACGATAAAAC‑GGTTTATGaa                                                                                                                                             >  2:719593/1‑149 (MQ=255)
             gATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAAtt                                                                                                                                           <  2:84278/150‑1 (MQ=255)
                  aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCCCGTGCTTCGGTAAACC‑GGGTATTGAAttctt                                                                                                                                        >  2:747312/1‑148 (MQ=255)
                  aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCg                                                                                                                                     <  1:777112/151‑1 (MQ=255)
                  aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCg                                                                                                                                     <  1:596647/151‑1 (MQ=255)
                      gTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑GGGTACTGAATTCTTTCGGtt                                                                                                                                  >  2:270016/1‑150 (MQ=255)
                        tCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCa                                                                                                                                <  2:403277/150‑1 (MQ=255)
                        tCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGATAAATTCCTAGGGTAGCCTAAAAC‑AGTTACAGAATTCTTTCggct                                                                                                                                  >  1:730868/1‑146 (MQ=255)
                          aaTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGt                                                                                                                             <  2:850178/151‑1 (MQ=255)
                          aaTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑GGTTACTGAATTCTTCCGGtt                                                                                                                                  >  1:801130/1‑146 (MQ=255)
                           aTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCGGGCGTTTGTACGGC‑GGTTAAATTCAGGTTGTCCGTTAAAGA‑GGGTACTGAATTTTTCCGGTTCaaa                                                                                                                              >  1:1306916/1‑148 (MQ=255)
                            ttttCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGTGCTACGATAAAAC‑GGGTACTGAATTCTTCCGGTTCAAGTa                                                                                                                            >  1:1233559/1‑150 (MQ=255)
                              ttCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTCAAGGGGTACGATAAAAC‑GGTTACTGAATTTTTTCGGTTCAAGTAccc                                                                                                                         >  1:1169355/1‑151 (MQ=255)
                                 cGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGc                                                                                                                      <  1:447607/151‑1 (MQ=255)
                                 cGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTGCCGGTTCAAGTACCCgg                                                                                                                       <  2:1167990/150‑1 (MQ=255)
                                   gCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑GGGTACTGAATTCTTTCGGGTCAAGTACCCgg                                                                                                                       >  1:1087382/1‑148 (MQ=255)
                                    ccTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGGGCTTCGGAAAAAC‑GGCTACTGAATTCTTCCCGCTCAAGGAACCCGCAg                                                                                                                    >  1:691560/1‑150 (MQ=255)
                                     cttcttCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAAATCAGGTGCTACGATAAAAC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCCGTg                                                                                                                  >  1:683094/1‑151 (MQ=255)
                                       tcttcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTgg                                                                                                                 <  1:808060/150‑1 (MQ=255)
                                        cttcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGTGCTTCGTTAAAAC‑GGTTTCTGAATTTTTTCGGTTCAAGTACCCGGccgggg                                                                                                                 >  2:358183/1‑144 (MQ=255)
                                          tcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGTGCTTCGATAAAGG‑GGTTACTTAATTTTTCCGGTTTAAGTACCCCGccggggga                                                                                                               >  2:1086441/1‑142 (MQ=255)
                                          tcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGGGCCACGATAAAGC‑GGTTACTGTATTTTTCCGGGTCAAATAcaccgccgcgggaa                                                                                                              >  2:230156/1‑137 (MQ=255)
                                             gCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGGGGTACGATAAAAC‑GGGTACTGAATTTTTCCGGGTCAAGTACCCCGCAGTGGgaa                                                                                                              >  1:545110/1‑147 (MQ=255)
                                                  gATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTcc                                                                                                      <  2:1013423/150‑1 (MQ=255)
                                                     aaCCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAAATCAGGTGCTACGATAAAAA‑GGGTTCTGTATTTTTCCGGTTCCAGTACCCGGCAGAGGAAAGCAGTTTCtta                                                                                                   >  1:764523/1‑148 (MQ=255)
                                                      aCCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTTAGGTGCTACGATAAAAC‑GGGTACTGAATTCTTCCCGGTCAAGGACCCGGCAGTGGGAAGCAGGTCCTaaa                                                                                                  >  1:684727/1‑150 (MQ=255)
                                                         gTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAAATTAGGGGCTACCAAAAAAC‑GGTTACCGAATTCTTTCGGTTCAAATACCCCGCAGTGGAAAGGAGGTCCTAAATCg                                                                                               >  1:1070526/1‑150 (MQ=255)
                                                          tttttCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑GGGTACTGAATTCTTCCGGTTCAAGTACCCCGCAGCGGAAAGCAGTTTCTaaa                                                                                                  >  1:924964/1‑146 (MQ=255)
                                                               ccTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGttt                                                                                        <  2:1086197/151‑1 (MQ=255)
                                                                  gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGttttt                                                                                      <  1:1115891/150‑1 (MQ=255)
                                                                  gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGtttt                                                                                       <  1:861071/149‑1 (MQ=255)
                                                                  gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGtttt                                                                                       <  1:1119839/149‑1 (MQ=255)
                                                                   tgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTTACGTGCTAAGATTAAGG‑GGGTTCTGAATTTTTTCGGTTCAAGAAGCCGGCCGGGGGAAACAGTTCCTAAATCGCGGGttttt                                                                                      >  1:1334179/1‑149 (MQ=255)
                                                                       tGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGGGCTACGATAAAGC‑GGGTAATGAATTCTTCCGGTTCCAGTAACCGGCAGTGGAAAGCAGGTCCTAAATCGAGCGGTTTTGTCttt                                                                                >  1:662688/1‑151 (MQ=255)
                                                                       tGCGCATTCCATCAGAAACGCTGTCACTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAGTTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGACCGTTTTTTTCtt                                                                                 <  2:433478/150‑1 (MQ=255)
                                                                         cgcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGTGCTACGATAAAAA‑GGGTACTGAATTCTTTCGGGTCAAGTAACCGGCAGAGGAAAGCAGGTTCTAAATAGAGCGCTTTTTTCttttt                                                                              >  2:1205873/1‑149 (MQ=255)
                                                                         cgcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGGTAAATTGAAGGGCCGCGATAAAAA‑GGCTCTTGAACTCTTACGGTTCAAGGAACCGGCCGGGGAAAGGAATTTCTAAATCGGGCGGTTTTGTCTTtg                                                                               >  1:310934/1‑150 (MQ=255)
                                                                          gcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTCAGGTGCTTCGATAAAACGGGTTATTAAATTTTTCGGGTT‑AAGTAACCGGGCGTGGAAAACAATTCTCAAAATGTGCGTTTTTGTCtt                                                                                 >  1:1015017/1‑147 (MQ=255)
                                                                          gcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTTAGGGCTACCATAAAAC‑GGGTACTGAATTCTTCCGGTTCCAATAACCCGCAGTGGGAAGCAGTTTCTAAATAGAGCGTtttttttttt                                                                                >  2:968151/1‑148 (MQ=255)
                                                                            aTTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTCAGGGGGTAAGATAAAAC‑GGTTACTGAATTTTTTCGGGTCAAGTAACCCGCAGTGGAAAGCCGTTCCCAAATTGAGCGGTTTTTTCTTTttgg                                                                            >  1:260524/1‑150 (MQ=255)
                                                                            aTTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGGGCTTCCATAAAAC‑CGGTACTTAATTCTTTCCGGTCAAGTAACCCGCAGTGGGAAAGAGTTTCTAAATCGAGCGTTTTTGTCTTTttgg                                                                            >  2:792716/1‑150 (MQ=255)
                                                                                    aGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGGTAAAATTCA‑GTGCCAAGGTAAAAG‑GGGTTCTGAATTTTTCCGGTTTAAGGAACCGGGCGAGGAAAGCAGATTCTAAATCGAGCGTTTTTGTCTTTGTGGTGTTCAc                                                                     >  2:1262737/1‑149 (MQ=255)
                                                                                    aGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCTGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgg                                                                   <  1:115778/151‑1 (MQ=255)
                                                                                     gAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑AGGTAATGAATTTTTTCGGGTCAAGGACCCCGCAGTGGAAAACAGTTCCTTAATTGAGCGTTTTTTTCTTTGTGGTGTTCAAggg                                                                  >  1:1241356/1‑151 (MQ=255)
                                                                                       aaCGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCACGCGCTACCATAAAAC‑GGGTACTGCATTCTTCCGGTTCAAGTACCCGGCAGGGGGAAGCAGTTCCTCAATTCAGCGTTTTTTTCTTTGTGGTGTCCACggggg                                                                >  2:125486/1‑151 (MQ=255)
                                                                                         cGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTCAGGTGCTTCGATAAAAG‑GGTTACTGAATTCTTCCGGTTTTAGAACCCCGCCGGGGAAAGCAGGTCCTAAAAAGCGCGCTTTTTTTTTTGTGGTGTTCCCcgggg                                                                >  1:348139/1‑149 (MQ=255)
                                                                                          gCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTTAATTCAGGTGCTACGATAAAAC‑GGTTACTGAAATCTTTCGGGTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTTTGGTGTCCCCGGGGGCAt                                                             >  1:1150389/1‑151 (MQ=255)
                                                                                          gCTCTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCTGGTGCTTCGATAAAAC‑GGGTTCTTAATTCTTCCGGGTTAAGTACCCGGCAGTGGAAAGCAGTTTCTTAATCGAGCGTTTTTTTCTTTGTTGTGTCCACGGGGGCAt                                                             >  2:170140/1‑151 (MQ=255)
                                                                                           cTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGTGTCTCGATAAAAC‑GGGTACTGAAATTTTTCGGGTCTAGTACCCCGCAGTGGGAAGGAGTTTCTAAAATGAGCGTTTTTTTCTTTTTGGGGTCCACGGGGGCAtt                                                            >  2:628093/1‑151 (MQ=255)
                                                                                                 cttcttcttTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCAttttt                                                         <  2:1018201/148‑1 (MQ=255)
                                                                                                  ttcttcttTTCCGGCGTTTGTACGGC‑GGTTTACTTCAAGAGCTACGATTAAGC‑GGTTTCTGCCATTTTCCGGCTTTATTATCCGGGAGCGGGAAACAACTCCTCAATAGAGGGGTTTTGTTTTTGTGGGGTCCACAGGGGGGTTTTTActc                                                     >  1:1228905/1‑149 (MQ=255)
                                                                                                  ttcttcttTTCCGGCGTTTGTACGGC‑GGTTAAATTTAAGTGCTACGATAAAAC‑GGTTAATGAAATCTTCCGGGTCAAGGACCCGGCAGTGGAAAGCAGTTCCTAAATCGCGCGGTTTTGTGTTTGTGGGGTCCACGGGGGGATTTTTTCGc                                                     >  1:1083730/1‑151 (MQ=255)
                                                                                                    cttcttTTCCGGCGTTTGTACGGC‑GGTTAAATTTAGGGGCTACGATAAAAC‑GGGTTCTGAATTTTTCCGGGTCAAGTACCCGGCAGTGGGAAACAGTTCCTAAATCGAGCGGTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGcc                                                    >  1:932951/1‑150 (MQ=255)
                                                                                                       cttTTCCGGCGTTTGTACGGC‑GGTTAAATTCCGGTGCTTCGATAAAAG‑GGTTAATGAATTCTTCCGGTTTAAGGAACCGGCAGTGGGAAGCCGTTCCTAAATTGTGCGTTTTTGTCTTTTTGGTGTTCACGGGGGGATTTTTACGCcaca                                                 >  1:943143/1‑150 (MQ=255)
                                                                                                       cttTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGGGCTACGATAAAAC‑GGGTACTGAATTCTTCCGGGTCAAGGACCCCGCAGTGGGAAGCAGTTCCTAAATCGAGCGGTTTTTTCTTTGTGGGGTCCACGGGGGGATTTTTACGCcac                                                  >  1:499234/1‑149 (MQ=255)
                                                                                                         tttCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCcaca                                                 <  1:39859/148‑1 (MQ=255)
                                                                                                         tttCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACaaa                                               <  1:993530/150‑1 (MQ=255)
                                                                                                         tttCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑GGGTTCTTAATTCTTTCGGGTCAAGGACCCGGCAGTGGAAAGAAGTTTCCAAATCGAGCGGTTTTTTCTTTTTGGGGTCCACGGGGGCCTTTTTTCGCCCCaa                                                >  2:209603/1‑149 (MQ=255)
                                                                                                          ttCCGGCGTTTGTACGGC‑GGTTAAATTCATGGGTTACTGTAAAGC‑GGGTATTTAAATCTTCCGGGTCAAGGACCCGGCAGTGGAAAGCAGTTTCTTAAAACAGCGGTTTTGTTTTTTTGGGGTTTCCGGGGGGATTTTTACGCCACaaaa                                              >  2:1032813/1‑149 (MQ=255)
                                                                                                          ttCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTTTCTTTGTGGTGTCCACGGGGGCCTTTTTACGCCACaaa                                               >  1:944286/1‑149 (MQ=255)
                                                                                                                 gTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAAC‑GGGTTCTGAATTCTTTCGGGTTAAATACCCGGGAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTTTCTTTTTGGTGTCCACGGGGGCATTTTTACGCCCCAAATATAGt                                         >  2:25712/1‑148 (MQ=255)
                                                                                                                       aCGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTa                                <  2:1210480/151‑1 (MQ=255)
                                                                                                                           c‑ggTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATc                              >  1:1281330/1‑149 (MQ=255)
                                                                                                                             ggTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGc                           >  1:414220/1‑151 (MQ=255)
                                                                                                                              gTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCa                          <  1:494006/151‑1 (MQ=255)
                                                                                                                                tAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAAt                        >  2:1271414/1‑151 (MQ=255)
                                                                                                                                 aaaTTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGTGCATTTTTACGCCACAAATATAGTGCCGGAGTATCATCAAt                        <  2:772048/150‑1 (MQ=255)
                                                                                                                                   aTTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTc                     <  2:1180407/151‑1 (MQ=255)
                                                                                                                                      cAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                  <  1:1086500/151‑1 (MQ=255)
                                                                                                                                      cAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                  <  1:317348/151‑1 (MQ=255)
                                                                                                                                      cAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                  <  1:50149/151‑1 (MQ=255)
                                                                                                                                         gttcTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGATTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACACCACAAATATAGAGCCGGAGTATCAGAATTGTCATCGa                >  2:1232438/4‑150 (MQ=255)
                                                                                                                                            cTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTg            >  1:603554/1‑151 (MQ=255)
                                                                                                                                              aCGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTg            <  1:813303/149‑1 (MQ=255)
                                                                                                                                                 aTAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACt       >  2:602928/1‑151 (MQ=255)
                                                                                                                                                   aaaGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACtt      <  1:924024/150‑1 (MQ=255)
                                                                                                                                                   aaaGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACtt      <  1:1030854/150‑1 (MQ=255)
                                                                                                                                                   aaaGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTa     <  1:1200459/151‑1 (MQ=255)
                                                                                                                                                   aaaGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTa     <  1:353090/151‑1 (MQ=255)
                                                                                                                                                      gC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTg  >  1:412260/1‑151 (MQ=255)
                                                                                                                                                      gC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTg  >  2:882882/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGC‑GGTTAAATTCAGGTGCTACGATAAAGC‑GGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTG  >  NC_000913/1415246‑1415545

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