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→CGTrecERac 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

GCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTG  >  NC_000913/1415246‑1415545
                                                                                                                                                      |                                                                                                                                                       
gCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGc                                                                                                                                                        <  2:251790/151‑1 (MQ=255)
     aaCGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGtt                                                                                                                                                   <  1:56477/150‑1 (MQ=255)
     aaCGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTa                                                                                                                                                  <  1:384827/151‑1 (MQ=255)
     aaCGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTa                                                                                                                                                  <  1:1109789/151‑1 (MQ=255)
          tttGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTACCATAAAACGGGTACt                                                                                                                                                >  2:1030800/1‑148 (MQ=255)
           ttGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTATGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGaa                                                                                                                                             <  1:231055/150‑1 (MQ=255)
            tGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCATGTGCTACGATAAAACGGTTTATGaa                                                                                                                                             >  2:719593/1‑149 (MQ=255)
             gATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAAtt                                                                                                                                           <  2:84278/150‑1 (MQ=255)
                  aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCCCGTGCTTCGGTAAACCGGGTATTGAAttctt                                                                                                                                        >  2:747312/1‑148 (MQ=255)
                  aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCg                                                                                                                                     <  1:777112/151‑1 (MQ=255)
                  aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCg                                                                                                                                     <  1:596647/151‑1 (MQ=255)
                      gTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGGTACTGAATTCTTTCGGtt                                                                                                                                  >  2:270016/1‑150 (MQ=255)
                        tCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCa                                                                                                                                <  2:403277/150‑1 (MQ=255)
                        tCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGATAAATTCCTAGGGTAGCCTAAAACAGTTACAGAATTCTTTCggct                                                                                                                                  >  1:730868/1‑146 (MQ=255)
                          aaTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGt                                                                                                                             <  2:850178/151‑1 (MQ=255)
                          aaTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGTTACTGAATTCTTCCGGtt                                                                                                                                  >  1:801130/1‑146 (MQ=255)
                           aTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCGGGCGTTTGTACGGCGGTTAAATTCAGGTTGTCCGTTAAAGAGGGTACTGAATTTTTCCGGTTCaaa                                                                                                                              >  1:1306916/1‑148 (MQ=255)
                            ttttCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGCTACGATAAAACGGGTACTGAATTCTTCCGGTTCAAGTa                                                                                                                            >  1:1233559/1‑150 (MQ=255)
                              ttCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAAGGGGTACGATAAAACGGTTACTGAATTTTTTCGGTTCAAGTAccc                                                                                                                         >  1:1169355/1‑151 (MQ=255)
                                 cGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGc                                                                                                                      <  1:447607/151‑1 (MQ=255)
                                 cGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTGCCGGTTCAAGTACCCgg                                                                                                                       <  2:1167990/150‑1 (MQ=255)
                                   gCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGGTACTGAATTCTTTCGGGTCAAGTACCCgg                                                                                                                       >  1:1087382/1‑148 (MQ=255)
                                    ccTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGGGCTTCGGAAAAACGGCTACTGAATTCTTCCCGCTCAAGGAACCCGCAg                                                                                                                    >  1:691560/1‑150 (MQ=255)
                                     cttcttCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAAATCAGGTGCTACGATAAAACGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCCGTg                                                                                                                  >  1:683094/1‑151 (MQ=255)
                                       tcttcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTgg                                                                                                                 <  1:808060/150‑1 (MQ=255)
                                        cttcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGCTTCGTTAAAACGGTTTCTGAATTTTTTCGGTTCAAGTACCCGGccgggg                                                                                                                 >  2:358183/1‑144 (MQ=255)
                                          tcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGCTTCGATAAAGGGGTTACTTAATTTTTCCGGTTTAAGTACCCCGccggggga                                                                                                               >  2:1086441/1‑142 (MQ=255)
                                          tcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGGGCCACGATAAAGCGGTTACTGTATTTTTCCGGGTCAAATAcaccgccgcgggaa                                                                                                              >  2:230156/1‑137 (MQ=255)
                                             gCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGGGGTACGATAAAACGGGTACTGAATTTTTCCGGGTCAAGTACCCCGCAGTGGgaa                                                                                                              >  1:545110/1‑147 (MQ=255)
                                                  gATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTcc                                                                                                      <  2:1013423/150‑1 (MQ=255)
                                                     aaCCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAAATCAGGTGCTACGATAAAAAGGGTTCTGTATTTTTCCGGTTCCAGTACCCGGCAGAGGAAAGCAGTTTCtta                                                                                                   >  1:764523/1‑148 (MQ=255)
                                                      aCCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTTAGGTGCTACGATAAAACGGGTACTGAATTCTTCCCGGTCAAGGACCCGGCAGTGGGAAGCAGGTCCTaaa                                                                                                  >  1:684727/1‑150 (MQ=255)
                                                         gTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAAATTAGGGGCTACCAAAAAACGGTTACCGAATTCTTTCGGTTCAAATACCCCGCAGTGGAAAGGAGGTCCTAAATCg                                                                                               >  1:1070526/1‑150 (MQ=255)
                                                          tttttCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGGTACTGAATTCTTCCGGTTCAAGTACCCCGCAGCGGAAAGCAGTTTCTaaa                                                                                                  >  1:924964/1‑146 (MQ=255)
                                                               ccTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGttt                                                                                        <  2:1086197/151‑1 (MQ=255)
                                                                  gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGttttt                                                                                      <  1:1115891/150‑1 (MQ=255)
                                                                  gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGtttt                                                                                       <  1:861071/149‑1 (MQ=255)
                                                                  gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGtttt                                                                                       <  1:1119839/149‑1 (MQ=255)
                                                                   tgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTTACGTGCTAAGATTAAGGGGGTTCTGAATTTTTTCGGTTCAAGAAGCCGGCCGGGGGAAACAGTTCCTAAATCGCGGGttttt                                                                                      >  1:1334179/1‑149 (MQ=255)
                                                                       tGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGGGCTACGATAAAGCGGGTAATGAATTCTTCCGGTTCCAGTAACCGGCAGTGGAAAGCAGGTCCTAAATCGAGCGGTTTTGTCttt                                                                                >  1:662688/1‑151 (MQ=255)
                                                                       tGCGCATTCCATCAGAAACGCTGTCACTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAGTTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGACCGTTTTTTTCtt                                                                                 <  2:433478/150‑1 (MQ=255)
                                                                         cgcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGCTACGATAAAAAGGGTACTGAATTCTTTCGGGTCAAGTAACCGGCAGAGGAAAGCAGGTTCTAAATAGAGCGCTTTTTTCttttt                                                                              >  2:1205873/1‑149 (MQ=255)
                                                                         cgcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGGTAAATTGAAGGGCCGCGATAAAAAGGCTCTTGAACTCTTACGGTTCAAGGAACCGGCCGGGGAAAGGAATTTCTAAATCGGGCGGTTTTGTCTTtg                                                                               >  1:310934/1‑150 (MQ=255)
                                                                          gcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTTCGATAAAACGGGTTATTAAATTTTTCGGGTTAAGTAACCGGGCGTGGAAAACAATTCTCAAAATGTGCGTTTTTGTCtt                                                                                 >  1:1015017/1‑147 (MQ=255)
                                                                          gcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTTAGGGCTACCATAAAACGGGTACTGAATTCTTCCGGTTCCAATAACCCGCAGTGGGAAGCAGTTTCTAAATAGAGCGTtttttttttt                                                                                >  2:968151/1‑148 (MQ=255)
                                                                            aTTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGGGGTAAGATAAAACGGTTACTGAATTTTTTCGGGTCAAGTAACCCGCAGTGGAAAGCCGTTCCCAAATTGAGCGGTTTTTTCTTTttgg                                                                            >  1:260524/1‑150 (MQ=255)
                                                                            aTTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGGGCTTCCATAAAACCGGTACTTAATTCTTTCCGGTCAAGTAACCCGCAGTGGGAAAGAGTTTCTAAATCGAGCGTTTTTGTCTTTttgg                                                                            >  2:792716/1‑150 (MQ=255)
                                                                                    aGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGGTAAAATTCAGTGCCAAGGTAAAAGGGGTTCTGAATTTTTCCGGTTTAAGGAACCGGGCGAGGAAAGCAGATTCTAAATCGAGCGTTTTTGTCTTTGTGGTGTTCAc                                                                     >  2:1262737/1‑149 (MQ=255)
                                                                                    aGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCTGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgg                                                                   <  1:115778/151‑1 (MQ=255)
                                                                                     gAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACAGGTAATGAATTTTTTCGGGTCAAGGACCCCGCAGTGGAAAACAGTTCCTTAATTGAGCGTTTTTTTCTTTGTGGTGTTCAAggg                                                                  >  1:1241356/1‑151 (MQ=255)
                                                                                       aaCGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCACGCGCTACCATAAAACGGGTACTGCATTCTTCCGGTTCAAGTACCCGGCAGGGGGAAGCAGTTCCTCAATTCAGCGTTTTTTTCTTTGTGGTGTCCACggggg                                                                >  2:125486/1‑151 (MQ=255)
                                                                                         cGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTTCGATAAAAGGGTTACTGAATTCTTCCGGTTTTAGAACCCCGCCGGGGAAAGCAGGTCCTAAAAAGCGCGCTTTTTTTTTTGTGGTGTTCCCcgggg                                                                >  1:348139/1‑149 (MQ=255)
                                                                                          gCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTACGATAAAACGGTTACTGAAATCTTTCGGGTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTTTGGTGTCCCCGGGGGCAt                                                             >  1:1150389/1‑151 (MQ=255)
                                                                                          gCTCTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCTGGTGCTTCGATAAAACGGGTTCTTAATTCTTCCGGGTTAAGTACCCGGCAGTGGAAAGCAGTTTCTTAATCGAGCGTTTTTTTCTTTGTTGTGTCCACGGGGGCAt                                                             >  2:170140/1‑151 (MQ=255)
                                                                                           cTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGTCTCGATAAAACGGGTACTGAAATTTTTCGGGTCTAGTACCCCGCAGTGGGAAGGAGTTTCTAAAATGAGCGTTTTTTTCTTTTTGGGGTCCACGGGGGCAtt                                                            >  2:628093/1‑151 (MQ=255)
                                                                                                 cttcttcttTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCAttttt                                                         <  2:1018201/148‑1 (MQ=255)
                                                                                                  ttcttcttTTCCGGCGTTTGTACGGCGGTTTACTTCAAGAGCTACGATTAAGCGGTTTCTGCCATTTTCCGGCTTTATTATCCGGGAGCGGGAAACAACTCCTCAATAGAGGGGTTTTGTTTTTGTGGGGTCCACAGGGGGGTTTTTActc                                                     >  1:1228905/1‑149 (MQ=255)
                                                                                                  ttcttcttTTCCGGCGTTTGTACGGCGGTTAAATTTAAGTGCTACGATAAAACGGTTAATGAAATCTTCCGGGTCAAGGACCCGGCAGTGGAAAGCAGTTCCTAAATCGCGCGGTTTTGTGTTTGTGGGGTCCACGGGGGGATTTTTTCGc                                                     >  1:1083730/1‑151 (MQ=255)
                                                                                                    cttcttTTCCGGCGTTTGTACGGCGGTTAAATTTAGGGGCTACGATAAAACGGGTTCTGAATTTTTCCGGGTCAAGTACCCGGCAGTGGGAAACAGTTCCTAAATCGAGCGGTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGcc                                                    >  1:932951/1‑150 (MQ=255)
                                                                                                       cttTTCCGGCGTTTGTACGGCGGTTAAATTCCGGTGCTTCGATAAAAGGGTTAATGAATTCTTCCGGTTTAAGGAACCGGCAGTGGGAAGCCGTTCCTAAATTGTGCGTTTTTGTCTTTTTGGTGTTCACGGGGGGATTTTTACGCcaca                                                 >  1:943143/1‑150 (MQ=255)
                                                                                                       cttTTCCGGCGTTTGTACGGCGGTTAAATTCAGGGGCTACGATAAAACGGGTACTGAATTCTTCCGGGTCAAGGACCCCGCAGTGGGAAGCAGTTCCTAAATCGAGCGGTTTTTTCTTTGTGGGGTCCACGGGGGGATTTTTACGCcac                                                  >  1:499234/1‑149 (MQ=255)
                                                                                                         tttCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCcaca                                                 <  1:39859/148‑1 (MQ=255)
                                                                                                         tttCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACaaa                                               <  1:993530/150‑1 (MQ=255)
                                                                                                         tttCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGGTTCTTAATTCTTTCGGGTCAAGGACCCGGCAGTGGAAAGAAGTTTCCAAATCGAGCGGTTTTTTCTTTTTGGGGTCCACGGGGGCCTTTTTTCGCCCCaa                                                >  2:209603/1‑149 (MQ=255)
                                                                                                          ttCCGGCGTTTGTACGGCGGTTAAATTCATGGGTTACTGTAAAGCGGGTATTTAAATCTTCCGGGTCAAGGACCCGGCAGTGGAAAGCAGTTTCTTAAAACAGCGGTTTTGTTTTTTTGGGGTTTCCGGGGGGATTTTTACGCCACaaaa                                              >  2:1032813/1‑149 (MQ=255)
                                                                                                          ttCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTTTCTTTGTGGTGTCCACGGGGGCCTTTTTACGCCACaaa                                               >  1:944286/1‑149 (MQ=255)
                                                                                                                 gTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGGTTCTGAATTCTTTCGGGTTAAATACCCGGGAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTTTCTTTTTGGTGTCCACGGGGGCATTTTTACGCCCCAAATATAGt                                         >  2:25712/1‑148 (MQ=255)
                                                                                                                       aCGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTa                                <  2:1210480/151‑1 (MQ=255)
                                                                                                                           cggTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATc                              >  1:1281330/1‑149 (MQ=255)
                                                                                                                             ggTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGc                           >  1:414220/1‑151 (MQ=255)
                                                                                                                              gTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCa                          <  1:494006/151‑1 (MQ=255)
                                                                                                                                tAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAAt                        >  2:1271414/1‑151 (MQ=255)
                                                                                                                                 aaaTTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGTGCATTTTTACGCCACAAATATAGTGCCGGAGTATCATCAAt                        <  2:772048/150‑1 (MQ=255)
                                                                                                                                   aTTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTc                     <  2:1180407/151‑1 (MQ=255)
                                                                                                                                      cAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                  <  1:1086500/151‑1 (MQ=255)
                                                                                                                                      cAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                  <  1:317348/151‑1 (MQ=255)
                                                                                                                                      cAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                  <  1:50149/151‑1 (MQ=255)
                                                                                                                                         gttcTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGATTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACACCACAAATATAGAGCCGGAGTATCAGAATTGTCATCGa                >  2:1232438/4‑150 (MQ=255)
                                                                                                                                            cTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTg            >  1:603554/1‑151 (MQ=255)
                                                                                                                                              aCGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTg            <  1:813303/149‑1 (MQ=255)
                                                                                                                                                 aTAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACt       >  2:602928/1‑151 (MQ=255)
                                                                                                                                                   aaaGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACtt      <  1:924024/150‑1 (MQ=255)
                                                                                                                                                   aaaGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACtt      <  1:1030854/150‑1 (MQ=255)
                                                                                                                                                   aaaGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTa     <  1:1200459/151‑1 (MQ=255)
                                                                                                                                                   aaaGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTa     <  1:353090/151‑1 (MQ=255)
                                                                                                                                                      gCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTg  >  1:412260/1‑151 (MQ=255)
                                                                                                                                                      gCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTg  >  2:882882/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTG  >  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: