Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,4980GA61.4% 67.8 / 197.3 158V148I (GTT→ATT) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base A (51/46);  total (85/73)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 171 total aligned reads displayed.

TTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGG  >  NC_000913/580402‑580598
                                                                                                |                                                                                                    
         tcccgGATCTCAGGTGAAATTGCCAGTATTCTCGCCGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGaa                                                                                         <  1:4013856/98‑1 (MQ=0)
                                 aggaTTCTCGACAGGATCCCCCTGTCGGTGCAGGGGCATTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGaa                                                                 <  1:3448462/98‑1 (MQ=1)
                                                 ttcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTg                                                 >  2:1092628/6‑101 (MQ=0)
                                                  tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                <  1:9171237/97‑1 (MQ=0)
                                                   cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgta                                               <  2:3461525/98‑3 (MQ=0)
                                                   cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                               >  1:14720517/4‑101 (MQ=0)
                                                                      gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                            <  1:4220264/96‑1 (MQ=0)
                                                                        ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCt                          <  1:11778593/98‑1 (MQ=1)
                                                                              gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGa                    <  1:7831689/98‑1 (MQ=1)
                                                                                    nnnaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat              <  2:9534114/98‑1 (MQ=11)
                                                                                    nnnaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat              <  2:14418297/98‑1 (MQ=11)
                                                                                    gAAAACGGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat              >  1:14072888/1‑101 (MQ=11)
                                                                                    gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat              >  1:447903/1‑101 (MQ=11)
                                                                                    gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat              >  1:10657069/1‑101 (MQ=11)
                                                                                    gAAAACCGACATATTGATTTCCTCAAGAAGGACCTCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat              <  1:5189711/101‑1 (MQ=11)
                                                                                    gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat              <  1:12203651/101‑1 (MQ=11)
                                                                                    gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat              <  2:3109388/101‑1 (MQ=11)
                                                                                    gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat              >  1:6728672/1‑101 (MQ=11)
                                                                                     nnnaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc             <  2:8568709/98‑1 (MQ=14)
                                                                                     nnnaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCCGAGTGAATATATc             <  2:14295655/98‑1 (MQ=9)
                                                                                     aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc             <  1:14350887/101‑1 (MQ=31)
                                                                                     aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc             >  2:6159687/1‑101 (MQ=31)
                                                                                     aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc             >  1:9965410/1‑101 (MQ=31)
                                                                                     aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc             >  1:13932318/1‑101 (MQ=11)
                                                                                     aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc             <  1:9228391/101‑1 (MQ=11)
                                                                                     aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc             <  2:5780982/101‑1 (MQ=11)
                                                                                     aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc             <  1:2703871/101‑1 (MQ=11)
                                                                                     aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc             <  1:3591959/101‑1 (MQ=11)
                                                                                      aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg            <  1:5657184/101‑1 (MQ=32)
                                                                                      aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg            >  1:1797250/1‑101 (MQ=32)
                                                                                      aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg            <  1:9243754/101‑1 (MQ=32)
                                                                                      aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCg            >  1:3000928/1‑101 (MQ=17)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGATGGAATATATCg            >  2:6769479/1‑101 (MQ=11)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg            <  1:10068588/101‑1 (MQ=11)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg            >  1:4887818/1‑101 (MQ=11)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg            >  1:2120361/1‑101 (MQ=11)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg            <  2:6406781/101‑1 (MQ=11)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg            <  1:13376500/101‑1 (MQ=11)
                                                                                      aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg            <  1:6061366/101‑1 (MQ=11)
                                                                                       nnncGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa           <  2:12147090/98‑1 (MQ=11)
                                                                                       aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa           <  1:9945328/101‑1 (MQ=33)
                                                                                       aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGTTGAAATAATACCGGGGTTGCTGAGTGAATATATCGa           >  1:3368707/1‑101 (MQ=11)
                                                                                       aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa           <  1:1789826/101‑1 (MQ=11)
                                                                                       aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa           <  1:13176356/101‑1 (MQ=11)
                                                                                        aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATcnnn          >  2:13478024/1‑98 (MQ=14)
                                                                                        aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa          <  2:1797250/101‑1 (MQ=35)
                                                                                        aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa          >  1:5586525/1‑101 (MQ=35)
                                                                                        aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa          <  1:13424876/101‑1 (MQ=35)
                                                                                        aCCGACATATTGATTTCCTCAAGAAGGATATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa          <  1:13682067/101‑1 (MQ=11)
                                                                                        aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATcnnn          >  2:12264370/1‑98 (MQ=11)
                                                                                        aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa          >  1:8469272/1‑101 (MQ=11)
                                                                                        aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa          <  2:5724343/101‑1 (MQ=11)
                                                                                        aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa          <  1:722215/101‑1 (MQ=11)
                                                                                        aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa          >  1:13692666/1‑101 (MQ=11)
                                                                                         nnnaCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc         <  2:7334232/98‑1 (MQ=11)
                                                                                         ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc         >  2:6380918/1‑101 (MQ=35)
                                                                                         ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc         <  1:4096760/101‑1 (MQ=11)
                                                                                         ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc         >  1:83141/1‑101 (MQ=11)
                                                                                          nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaaaa        <  2:14720517/98‑3 (MQ=11)
                                                                                          nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa        <  2:14780816/98‑1 (MQ=11)
                                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAATCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa        >  1:10604540/1‑101 (MQ=38)
                                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa        >  2:297457/1‑101 (MQ=38)
                                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa        >  1:1244784/1‑101 (MQ=38)
                                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGTAATAATACCGGGGTTGCTGAGTGAATATATCGAACa        >  1:4789765/1‑101 (MQ=11)
                                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGannn        >  2:11176917/1‑98 (MQ=11)
                                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa        <  1:8108679/101‑1 (MQ=11)
                                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa        >  1:9908078/1‑101 (MQ=11)
                                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAAAAATACAGGGGTTGATGAGTGAATATATCGaaaa        >  1:3488136/1‑99 (MQ=2)
                                                                                           nnnaTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg       <  2:8932342/98‑1 (MQ=255)
                                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaannn       >  2:9043237/1‑98 (MQ=25)
                                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGCACAg       >  2:4937974/1‑101 (MQ=25)
                                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg       <  2:5416161/101‑1 (MQ=38)
                                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg       >  1:10930731/1‑101 (MQ=38)
                                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg       >  1:7449198/1‑101 (MQ=14)
                                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg       <  1:10275343/101‑1 (MQ=14)
                                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg       >  1:1732193/1‑101 (MQ=14)
                                                                                           gACATATTGATTTCCTCAAGAAGGACATCAGTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg       <  1:12780147/101‑1 (MQ=14)
                                                                                           ctcATGTTGATTTCATGAAACGGGAAATCATCAACGCCATCAACAAAGCAGCCGAGCTGGATGAACTGATCCCGGGGTTGCTGAGTGAATATATCGAACAg       <  2:4238117/99‑1 (MQ=14)
                                                                                            gcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt      <  2:5743167/100‑1 (MQ=14)
                                                                                            aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATAGATCGAACAGt      >  1:13112047/1‑101 (MQ=255)
                                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAcnnn      >  2:9195353/1‑98 (MQ=2)
                                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAgc      >  1:5833639/1‑100 (MQ=14)
                                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt      <  1:8951572/101‑1 (MQ=14)
                                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt      >  1:1255696/1‑101 (MQ=14)
                                                                                             nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc     <  2:12747280/98‑1 (MQ=255)
                                                                                             nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc     <  2:10697805/98‑1 (MQ=255)
                                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc     >  1:10121400/1‑101 (MQ=40)
                                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc     >  1:722159/1‑101 (MQ=40)
                                                                                             cATATTGATTTCGTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACannn     >  2:7616095/1‑98 (MQ=11)
                                                                                             cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc     <  2:758686/101‑1 (MQ=14)
                                                                                             cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc     >  1:11544639/1‑101 (MQ=14)
                                                                                             cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGGTGCTGAGTGAATATATCGAACannn     >  2:12579737/1‑98 (MQ=11)
                                                                                              aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa    >  2:508321/1‑101 (MQ=255)
                                                                                              aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa    >  2:3894571/1‑101 (MQ=255)
                                                                                              aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGCCa    >  1:8994500/1‑101 (MQ=255)
                                                                                               tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg   <  1:12550944/101‑1 (MQ=255)
                                                                                                gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg  <  1:6277127/101‑1 (MQ=255)
                                                                                                gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg  <  1:5451696/101‑1 (MQ=255)
                                                                                                gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg  <  1:6894158/101‑1 (MQ=255)
                                                                                                gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg  <  1:9727896/101‑1 (MQ=255)
                                                                                                |                                                                                                    
TTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGG  >  NC_000913/580402‑580598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑