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

TTCATCACTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATC  >  NC_000913/1415641‑1415931
                                                                                                                                                    |                                                                                                                                                  
tgcctcaCTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCGCTTCATCGAATGTTCCCTGGTGTGTTGCGTCGCCCATCACCg                                                                                                                                                 <  2:531321/145‑1 (MQ=255)
  catcaCTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGcacc                                                                                                                                            <  1:282785/151‑1 (MQ=255)
     caCTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACcaca                                                                                                                                         <  2:686780/151‑1 (MQ=255)
        tGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTc                                                                                                                                      >  1:793118/1‑151 (MQ=255)
           aaTCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAgg                                                                                                                                   >  2:582524/1‑151 (MQ=255)
             tCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCCCAGTCAGGGc                                                                                                                                 >  2:1277913/1‑151 (MQ=255)
              cgcgATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCa                                                                                                                                <  1:491056/151‑1 (MQ=255)
               gcgATGCGGATCGCTGCTAGCATTCTCATTGTGCGAATGTTCAGAGACTCGCATTTCCTCCGGATAATTTTCCGTAGCTTCAACCTGACTCTCTTCATCGAATGCTTCCTGGTATGTTGCGTCGCCCATCGCCGCACCACAgcaagcg                                                                                                                                  >  2:166539/1‑142 (MQ=255)
                 gATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCACCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCCTCCTCGAAAGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACTACAGTCAGGGCa                                                                                                                                <  2:39311/148‑1 (MQ=255)
                  aTGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGtt                                                                                                                             >  1:1444845/1‑150 (MQ=255)
                   tGCGGCTCGCTGCCAGCTTTCTCATTGTGCGGATGTTCAGCGTCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGACTGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACCGTCAGGGCAGTTAt                                                                                                                           >  2:478209/1‑151 (MQ=255)
                     cGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATcc                                                                                                                         >  2:1540727/1‑151 (MQ=255)
                       gATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATcccc                                                                                                                       >  2:694336/1‑151 (MQ=255)
                        aTCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCg                                                                                                                      >  1:782234/1‑151 (MQ=255)
                         tCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCg                                                                                                                      <  2:332326/150‑1 (MQ=255)
                         tCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCg                                                                                                                      <  2:721210/150‑1 (MQ=255)
                            cTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCCCCGCACCACAGTCAGGGCAGTTATCCCCGCCa                                                                                                                   <  2:874759/150‑1 (MQ=255)
                            cTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGcc                                                                                                                    >  1:227087/1‑149 (MQ=255)
                                     ttCTCATTGTGCGGCTGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGTATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCGGCCAGTCTGGCCGc                                                                                                         >  2:1090486/1‑151 (MQ=255)
                                      tctcATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGc                                                                                                         >  1:1431966/1‑150 (MQ=255)
                                           ttGTGGGGTTGTTCAGCCCCTTCCTTTCCCCCCGGATCTTTTTCTTTACCTCAACCCTGCCTCTCTCCATCAAATTTTTCCTGTTATTTGGCGCCCCCCCCCCCCGCCCCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGc                                                                                                     <  1:59061/149‑1 (MQ=255)
                                               ggggATTTTCAGCGCCTTCCATTCCCTCCGGATCTTTTTCCTTACCTTCACCCTGCCTCTCTTCATCGATTTTTTCCTGGTTTGTTGCGCCCCCCATCCCCGCCCCACAGTCAGGGCATTTATCCCCGCCAGTCTGGCCGCAGGCATTGCa                                                                                               <  1:784779/149‑1 (MQ=255)
                                               gCGGTTTTTCGGCCCCTCCCATTTCCCCCGGATCATTTCCCTTAGCTTCACCCCGCCCCTCTCCACAAATTTTTCCCTGTTAGTTGGCGCCGCCCCCCCCCGCCCCACGGTCAGGGCATTTACCCCCGCCAGTCTGGCCGCAgg                                                                                                      <  1:805397/144‑1 (MQ=255)
                                                    tGTTCAGCGCCTTCCATTTCCTCGGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCCAAGGCAtccaaggct                                                                                           >  2:1253369/1‑142 (MQ=255)
                                                      ttcccgcccTTCCATTTTCCCCCGGTTTTTTTTCTTTACTTCCACCCTCCCTCCCTCCACCGATTTTTTCCGGTTATGTTGCGGCCCCCCCCCCCCCCCCACAGCCAGGGAATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTa                                                                                          <  2:58394/143‑1 (MQ=255)
                                                       tCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTAttt                                                                                       >  2:282511/1‑151 (MQ=255)
                                                           ccccTCCCTTTCCCCCCGAATCTTTTTCTTTATCTTCACCCTGCCTCTCTTCATCAATTTTTTCCTGGTATGTGGCGCCCCCCACCCCCCCCCCACAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCg                                                                                    <  1:1591971/148‑1 (MQ=255)
                                                                     tttCCCCCGAATCTTTTTCCTTACCTCCACCCTGCCTCTCTCCACCGATTGTTTCCTGGTTTTTGGCGCCCCCCACCCCCCCCCCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTc                                                                                <  1:1066997/145‑1 (MQ=255)
                                                                     ttccctccGGTTCTTTTTCCTTACCTTCACCCTGCCTCTCTTCATCGATTGTTTCGTGGTATGTGGCGCCCCCCACCCCCGCACCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTg                                                                             <  1:76952/145‑1 (MQ=255)
                                                                      tccccccGGATCTTTTTCCTTACCAACACCCTCCCTCCCTTCACCGAGTTTTACGTGGTGTGTGGCGCCCCCCCCCCCCCCCCCACAGGCGGGGCAGTTTTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCt                                                                              <  2:1492299/141‑1 (MQ=255)
                                                                       tcctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGc                                                                         >  2:197562/1‑150 (MQ=255)
                                                                        cctccttttctttttCCTTCGCTCCACCCTGCCTCTCTCCACCTTTTTTTTCCTTTTATTTTGCGCCACCCCTGCCCCCCCCCCAGTCAGGGCATTGCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTg                                                                          <  1:560795/137‑1 (MQ=255)
                                                                        cctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCac                                                                       >  1:447236/1‑151 (MQ=255)
                                                                        cctccGGATCATTTTCCTTACCTTAAACCTGCCTCTCTTCATCGAATGATTCCTGGTATATTGCGTCGCACATCACCGCACCACCGTCAGGGCAGTTATCACTGCCAGTCTGGCTGCAGGCATTGCAGGCTCTTTCCGGTTCCTGTTGc                                                                         >  2:1277917/1‑149 (MQ=255)
                                                                         ctcggtttCATTTTCTTTAGCTTAACCCTGCCTCTCTTCATCGAATGTTTCCTGGTATGTTGCGCCGCCCATCCCCGCCCCACAGTCGGGGCAGTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTg                                                                          <  1:357572/140‑1 (MQ=255)
                                                                         ctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCact                                                                      >  1:572216/1‑151 (MQ=255)
                                                                         ctccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCact                                                                      >  1:421162/1‑151 (MQ=255)
                                                                          tcgtgttctttttCTTTAACTCAACCCTGTCTCTCTTTATCGAATTTGCCCTGGTATTTGGCGCCCCCAACCCCGCCCCCAAAGACAGGCCATTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGt                                                                            <  1:1030562/135‑1 (MQ=255)
                                                                               aTCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGc                                                                >  2:1473659/1‑151 (MQ=255)
                                                                               aTCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGc                                                                >  2:34895/1‑151 (MQ=255)
                                                                                tCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCt                                                               >  2:1306112/1‑151 (MQ=255)
                                                                                   ttttCTTTAGCTCCACCCTGCCTCTCTTCATCGATTTTTTCCTGGTTTGTTGCGCCCCCCACCCCCCCCCCACAGGCAGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTc                                                              <  1:1048981/149‑1 (MQ=255)
                                                                                   ttttCCTTAGCTCCACCCTGCCCCTCTCCATCGATTTTTTCTTGGTTTGTTCCGTCCCCCCCCCCCGCCCCACGGTCAGGGCATTTCTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCa                                                             <  2:456460/150‑1 (MQ=255)
                                                                                   ttttCCTTACCTCCACCCTGCCTCTCTTCATCGAATTTTTCCTGTTATTTTGCGCCCCCCACCCCCCCCCCACAGTCAGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCa                                                             <  1:1390446/150‑1 (MQ=255)
                                                                                   tttcCCTTACCTACACCCTGCCTCTCTTCATCGAATTTTTCTGGTTTTGTGGCGCCCCCCACCCCCCCCCCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCt                                                               <  1:680207/148‑1 (MQ=255)
                                                                                    tttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAgg                                                           >  1:730543/1‑151 (MQ=255)
                                                                                     ttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAgg                                                           >  2:1584692/1‑150 (MQ=255)
                                                                                       ccTTATCTTCAACCTGACTCTCTTCACCGATTGTTTCTTGGTATGTTGCGTCGCCCACCCCCGCCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGtt                                                         <  1:1129615/150‑1 (MQ=255)
                                                                                       ccTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTg                                                        >  1:1061968/1‑151 (MQ=255)
                                                                                         ttACCTTCAACTTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCCCCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCCGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTCCCCTACTGGCCCAGGtt                                                         >  2:203370/1‑148 (MQ=255)
                                                                                           aGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTc                                                    >  2:779450/1‑151 (MQ=255)
                                                                                             cTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCat                                                  >  1:761706/1‑151 (MQ=255)
                                                                                                 aaCCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCt                                              >  1:399245/1‑151 (MQ=255)
                                                                                                   ccTGTCTCCCTTCACCAAATTTTCCGGGTTTTGTGGCGCCCCCCATCCCCGCCCCAGAGCCGGGGCATTAATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTg                                             <  1:492345/150‑1 (MQ=255)
                                                                                                   ccTGCCTCTTTTCATCAAATGTTTCCTGTTGTGTTGCTTCGCCCCCCCCCGCACCACAGTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTgg                                            <  1:914452/151‑1 (MQ=255)
                                                                                                               cttCGATTTTTCCCTGGTTTGTGGCGTCCCCCTCCCCCCCACCACGGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGtttt                                   <  2:1207496/146‑1 (MQ=255)
                                                                                                               catagattttttCTTGGTATGTGGCGCCGCCCATCCCCGCCCCACAGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGtttttt                                 <  1:745300/141‑1 (MQ=255)
                                                                                                               cATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTc                                >  1:1050210/1‑151 (MQ=255)
                                                                                                                accGAATTTTTCCTGGTATGTTGCGCCCCCCACCCCCCCCCCACAGTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTc                                <  1:250499/148‑1 (MQ=255)
                                                                                                                  caaaTGTTTCCTGTTAGGTGGCCTCCCCCTCCCCCGCACCACAGTCGGTGCATTCTCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGttttt                                  <  2:118991/144‑1 (MQ=255)
                                                                                                                  cGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGt                             >  2:830870/1‑151 (MQ=255)
                                                                                                                   gAATGATTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCACGGCAGTTATCCCCGCCAGTCTGGCCGCACGCACTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGACTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGtt                            >  2:1314767/1‑151 (MQ=255)
                                                                                                                    caTGTTTCCTGGTAGGTTGCGCCGCCAGCCCCCGCCCAACATTCAGCGCGTTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTCTTCCg                              <  2:168601/147‑1 (MQ=255)
                                                                                                                        tttCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctgg                       >  1:1422779/1‑151 (MQ=255)
                                                                                                                         ttCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggc                      >  2:1453324/1‑151 (MQ=255)
                                                                                                                         ttCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggc                      >  1:1334341/1‑151 (MQ=255)
                                                                                                                          tCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTATTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggc                      >  2:397558/1‑150 (MQ=255)
                                                                                                                            ccTGGTATGTTGCGCCCCCCACCCCCCCCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggct                     <  1:792570/150‑1 (MQ=255)
                                                                                                                              tGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGt                  >  2:1304169/1‑151 (MQ=255)
                                                                                                                               tgtttttttGCGCCCCCCACCCCCCCCCCACAGTGAGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggctg                    <  1:367821/142‑1 (MQ=255)
                                                                                                                                  tgggTTGCGCCGCCCCCCCCCGCCCCAAGGCCGGGGGATTCATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGtt                 <  1:121498/146‑1 (MQ=255)
                                                                                                                                   aTGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCCGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCt               >  2:701189/1‑150 (MQ=255)
                                                                                                                                      ttGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCACGCCCGTCTGGCCGCAGGCATTGCAGGCTATTTCAGGTTCCTGTTGCCCTACTGGCTCAGGGTGATTCATCTCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGt            >  2:1339509/1‑150 (MQ=255)
                                                                                                                                               ccATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATc  >  1:1542011/1‑151 (MQ=255)
                                                                                                                                               ccATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATc  >  2:166592/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCATCACTGCAATCGCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATC  >  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: