Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,415,7860AC37.8% 123.7 / 18.1 75V534G (GTG→GGG) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (39/7);  new base C (1/27);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.10e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCACTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATC  >  NC_000913/1415641‑1415931
                                                                                                                                                    |                                                                                                                                                  
tgcctcaCTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCGCTTCATCGAATGTTC‑CCTGG‑TGTGTTGCGTCGCCCATCACCg                                                                                                                                                 <  2:531321/145‑1 (MQ=255)
  catcaCTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑cacc                                                                                                                                            <  1:282785/151‑1 (MQ=255)
     caCTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACcaca                                                                                                                                         <  2:686780/151‑1 (MQ=255)
        tGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTc                                                                                                                                      >  1:793118/1‑151 (MQ=255)
           aaTCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAgg                                                                                                                                   >  2:582524/1‑151 (MQ=255)
             tCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCCCAGTCAGGGc                                                                                                                                 >  2:1277913/1‑151 (MQ=255)
              cgcgATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCa                                                                                                                                <  1:491056/151‑1 (MQ=255)
               gcgATGCGGATCGCTGCTAGCATTCTCATTGTGCGAATGTTCAGAGACTCGCA‑TTTCCTCCGGATAATTTTCCGTAGCTTCAACCTGACTCTCTTCATCGAATGCTT‑CCTGG‑TATGTTGCGTCGCCCATCGCCG‑CACCACAgcaagcg                                                                                                                                  >  2:166539/1‑142 (MQ=255)
                 gATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCACCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCCTCCTCGAAAGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACTACAGTCAGGGCa                                                                                                                                <  2:39311/148‑1 (MQ=255)
                  aTGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGtt                                                                                                                             >  1:1444845/1‑150 (MQ=255)
                   tGCGGCTCGCTGCCAGCTTTCTCATTGTGCGGATGTTCAGCGTCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGACTGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACCGTCAGGGCAGTTAt                                                                                                                           >  2:478209/1‑151 (MQ=255)
                     cGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATcc                                                                                                                         >  2:1540727/1‑151 (MQ=255)
                       gATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATcccc                                                                                                                       >  2:694336/1‑151 (MQ=255)
                        aTCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCg                                                                                                                      >  1:782234/1‑151 (MQ=255)
                         tCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCg                                                                                                                      <  2:332326/150‑1 (MQ=255)
                         tCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCg                                                                                                                      <  2:721210/150‑1 (MQ=255)
                            cTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCCCCG‑CACCACAGTCAGGGCAGTTATCCCCGCCa                                                                                                                   <  2:874759/150‑1 (MQ=255)
                            cTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGcc                                                                                                                    >  1:227087/1‑149 (MQ=255)
                                     ttCTCATTGTGCGGCTGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGTATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCGGCCAGTCTGGCCGc                                                                                                         >  2:1090486/1‑151 (MQ=255)
                                      tctcATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGc                                                                                                         >  1:1431966/1‑150 (MQ=255)
                                           ttGTGGGGTTGTTCAGCCCCTTCCT‑TTCCCCCCGGATCTTTTTCTTTACCTCAACCCTGCCTCTCTCCATCAAATTTTT‑CCTGT‑TATTTGGCGCCCCCCCCCCCCG‑CCCCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGc                                                                                                     <  1:59061/149‑1 (MQ=255)
                                               ggggATTTTCAGCGCCTTCCA‑TTCCCTCCGGATCTTTTTCCTTACCTTCACCCTGCCTCTCTTCATCGATTTTTT‑CCTGG‑TTTGTTGCGCCCCCCATCCCCG‑CCCCACAGTCAGGGCATTTATCCCCGCCAGTCTGGCCGCAGGCATTGCa                                                                                               <  1:784779/149‑1 (MQ=255)
                                               gCGGTTTTTCGGCCCCTCCCA‑TTTCCCCCGGATCATTTCCCTTAGCTTCACCCCGCCCCTCTCCA‑CAAATTTTTCCCTGT‑TAGTTGGCGCCGCCCCCCCCCG‑CCCCACGGTCAGGGCATTTACCCCCGCCAGTCTGGCCGCAgg                                                                                                      <  1:805397/144‑1 (MQ=255)
                                                    tGTTCAGCGCCTTCCA‑TTTCCTCGGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCCAAGGCAtccaaggct                                                                                           >  2:1253369/1‑142 (MQ=255)
                                                      ttcccgcccTTCCATTTTCCCCCGGTTTTTTTTCTTTACTTCCACCCTCCCTCCCTCCACCGATTTTTT‑CCGGT‑TATGTTGCGGCCCCCCCCCCCC‑CCCCACAGCCAGGGAATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTa                                                                                          <  2:58394/143‑1 (MQ=255)
                                                       tCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTAttt                                                                                       >  2:282511/1‑151 (MQ=255)
                                                           ccccTCCCT‑TTCCCCCCGAATCTTTTTCTTTATCTTCACCCTGCCTCTCTTCATCAATTTTTT‑CCTGG‑TATGTGGCGCCCCCCACCCCCC‑CCCCACAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCg                                                                                    <  1:1591971/148‑1 (MQ=255)
                                                                     tttCCCCCGAATCTTTTTCCTTACCTCCACCCTGCCTCTCTCCACCGATTGTTT‑CCTGG‑TTTTTGGCGCCCCCCACCCCCC‑CCCCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTc                                                                                <  1:1066997/145‑1 (MQ=255)
                                                                     ttccctccGGTTCTTTTTCCTTACCTTCACCCTGCCTCTCTTCATCGATTGTTT‑CGTGG‑TATGTGGCGCCCCCCACCCCCG‑CACCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTg                                                                             <  1:76952/145‑1 (MQ=255)
                                                                      tccccccGGATCTTTTTCCTTACCAACACCCTCCCTCCCTTCACCGAGTTTTA‑CGTGG‑TGTGTGGCGCCCCCCCCCCCCC‑CCCCACAGGCGGGGCAGTTTTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCt                                                                              <  2:1492299/141‑1 (MQ=255)
                                                                       tcctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGc                                                                         >  2:197562/1‑150 (MQ=255)
                                                                        cctccttttctttttCCTTCGCTCCACCCTGCCTCTCTCCACCTTTTTTTT‑CCTTT‑TATTTTGCGCCACCCCTGCCCC‑CCCCCCAGTCAGGGCATTGCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTg                                                                          <  1:560795/137‑1 (MQ=255)
                                                                        cctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCac                                                                       >  1:447236/1‑151 (MQ=255)
                                                                        cctccGGATCATTTTCCTTACCTTAAACCTGCCTCTCTTCATCGAATGATT‑CCTGG‑TATATTGCGTCGCACATCACCG‑CACCACCGTCAGGGCAGTTATCACTGCCAGTCTGGCTGCAGGCATTGCAGGCTCTTTCCGGTTCCTGTTGc                                                                         >  2:1277917/1‑149 (MQ=255)
                                                                         ctcggtttCATTTTCTTTAGCTTAACCCTGCCTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGCCGCCCATCCCCG‑CCCCACAGTCGGGGCAGTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTg                                                                          <  1:357572/140‑1 (MQ=255)
                                                                         ctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCact                                                                      >  1:572216/1‑151 (MQ=255)
                                                                         ctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCact                                                                      >  1:421162/1‑151 (MQ=255)
                                                                          tcgtgttctttttCTTTAACTCAACCCTGTCTCTCTTTATCGAATTTGC‑CCTGG‑TATTTGGCG‑CCCCCAACCCCGCCCCCAAAGACAGGCCATTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGt                                                                            <  1:1030562/135‑1 (MQ=255)
                                                                               aTCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGc                                                                >  2:1473659/1‑151 (MQ=255)
                                                                               aTCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGc                                                                >  2:34895/1‑151 (MQ=255)
                                                                                tCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCt                                                               >  2:1306112/1‑151 (MQ=255)
                                                                                   ttttCTTTAGCTCCACCCTGCCTCTCTTCATCGATTTTTT‑CCTGG‑TTTGTTGCGCCCCCCACCCCCC‑CCCCACAGGCAGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTc                                                              <  1:1048981/149‑1 (MQ=255)
                                                                                   ttttCCTTAGCTCCACCCTGCCCCTCTCCATCGATTTTTT‑CTTGG‑TTTGTTCCGTCCCCCCCCCCCG‑CCCCACGGTCAGGGCATTTCTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCa                                                             <  2:456460/150‑1 (MQ=255)
                                                                                   ttttCCTTACCTCCACCCTGCCTCTCTTCATCGAATTTTT‑CCTGT‑TATTTTGCGCCCCCCACCCCCC‑CCCCACAGTCAGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCa                                                             <  1:1390446/150‑1 (MQ=255)
                                                                                   tttcCCTTACCTACACCCTGCCTCTCTTCATCGAATTTTT‑CTGGT‑TTTGTGGCGCCCCCCACCCCCC‑CCCCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCt                                                               <  1:680207/148‑1 (MQ=255)
                                                                                    tttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAgg                                                           >  1:730543/1‑151 (MQ=255)
                                                                                     ttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAgg                                                           >  2:1584692/1‑150 (MQ=255)
                                                                                       ccTTATCTTCAACCTGACTCTCTTCACCGATTGTTT‑CTTGG‑TATGTTGCGTCGCCCACCCCCG‑CCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGtt                                                         <  1:1129615/150‑1 (MQ=255)
                                                                                       ccTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTg                                                        >  1:1061968/1‑151 (MQ=255)
                                                                                         ttACCTTCAACTTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCCCCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCCGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTCCCCTACTGGCCCAGGtt                                                         >  2:203370/1‑148 (MQ=255)
                                                                                           aGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTc                                                    >  2:779450/1‑151 (MQ=255)
                                                                                             cTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCat                                                  >  1:761706/1‑151 (MQ=255)
                                                                                                 aaCCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCt                                              >  1:399245/1‑151 (MQ=255)
                                                                                                   ccTGTCTCCCTTCACCAAAT‑TTT‑CCGGGTTTTGTGGCGCCCCCCATCCCCG‑CCCCAGAGCCGGGGCATTAATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTg                                             <  1:492345/150‑1 (MQ=255)
                                                                                                   ccTGCCTCTTTTCATCAAATGTTT‑CCTGT‑TGTGTTGCTTCGCCCCCCCCCG‑CACCACAGTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTgg                                            <  1:914452/151‑1 (MQ=255)
                                                                                                               cttCGATTTTTC‑CCTGG‑TTTGTGGCGTCCCCCTCCCCCC‑CACCACGGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGtttt                                   <  2:1207496/146‑1 (MQ=255)
                                                                                                               catagatttttt‑CTTGG‑TATGTGGCGCCGCCCATCCCCG‑CCCCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGtttttt                                 <  1:745300/141‑1 (MQ=255)
                                                                                                               cATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTc                                >  1:1050210/1‑151 (MQ=255)
                                                                                                                accGAATTTTT‑CCTGG‑TATGTTGCGCCCCCCACCCCCC‑CCCCACAGTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTc                                <  1:250499/148‑1 (MQ=255)
                                                                                                                  caaaTGTTT‑CCTGT‑TAGGTGGCCTCCCCCTCCCCCG‑CACCACAGTCGGTGCATTCTCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGttttt                                  <  2:118991/144‑1 (MQ=255)
                                                                                                                  cGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGt                             >  2:830870/1‑151 (MQ=255)
                                                                                                                   gAATGATT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCACGGCAGTTATCCCCGCCAGTCTGGCCGCACGCACTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGACTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGtt                            >  2:1314767/1‑151 (MQ=255)
                                                                                                                    caTGTTT‑CCTGG‑TAGGTTGCGCCGCCAGCCCCCG‑CCCAACATTCAGCGCGTTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTCTTCCg                              <  2:168601/147‑1 (MQ=255)
                                                                                                                        ttt‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctgg                       >  1:1422779/1‑151 (MQ=255)
                                                                                                                         tt‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggc                      >  2:1453324/1‑151 (MQ=255)
                                                                                                                         tt‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggc                      >  1:1334341/1‑151 (MQ=255)
                                                                                                                          t‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTATTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggc                      >  2:397558/1‑150 (MQ=255)
                                                                                                                            ccTGG‑TATGTTGCGCCCCCCACCCCCC‑CCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggct                     <  1:792570/150‑1 (MQ=255)
                                                                                                                              tGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGt                  >  2:1304169/1‑151 (MQ=255)
                                                                                                                               tgtttttttGCGCCCCCCACCCCCC‑CCCCACAGTGAGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggctg                    <  1:367821/142‑1 (MQ=255)
                                                                                                                                  tgggTTGCGCCGCCCCCCCCCG‑CCCCAAGGCCGGGGGATTCATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGtt                 <  1:121498/146‑1 (MQ=255)
                                                                                                                                   aTGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCCGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCt               >  2:701189/1‑150 (MQ=255)
                                                                                                                                      ttGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCACGCCCGTCTGGCCGCAGGCATTGCAGGCTATTTCAGGTTCCTGTTGCCCTACTGGCTCAGGGTGATTCATCTCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGt            >  2:1339509/1‑150 (MQ=255)
                                                                                                                                               ccATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATc  >  1:1542011/1‑151 (MQ=255)
                                                                                                                                               ccATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATc  >  2:166592/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCATCACTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCA‑TTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTT‑CCTGG‑TATGTTGCGTCGCCCATCACCG‑CACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATC  >  NC_000913/1415641‑1415931

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