Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5120AG52.9% 232.7 / inf 278K152K (AAA→AAG) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base A (62/69);  new base G (67/80);  total (129/149)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 315 total aligned reads displayed.

AGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGG  >  NC_000913/580435‑580612
                                                                             |                                                                                                    
aggaTTCTCGACAGGATCCCCCTGTCGGTGCAGGGGCATTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGaa                                                                               <  1:3448462/98‑1 (MQ=1)
                ttcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTg                                                               >  2:1092628/6‑101 (MQ=0)
                  cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgta                                                             <  2:3461525/98‑3 (MQ=0)
                  cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                             >  1:14720517/4‑101 (MQ=0)
                                     gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                          >  2:3030520/6‑101 (MQ=0)
                                     gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                          <  1:12507530/96‑1 (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:14418297/98‑1 (MQ=11)
                                                   gAAATCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                            <  1:10898400/101‑1 (MQ=11)
                                                   gAAAACGGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                            >  1:14072888/1‑101 (MQ=11)
                                                    aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           >  1:13932318/1‑101 (MQ=11)
                                                    aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           <  1:3591959/101‑1 (MQ=11)
                                                    aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           >  1:7309166/1‑101 (MQ=11)
                                                    aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           <  1:9228391/101‑1 (MQ=11)
                                                     nnnccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                          <  2:9181481/98‑1 (MQ=17)
                                                     aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCg                          >  1:3000928/1‑101 (MQ=17)
                                                     aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                          <  1:10068588/101‑1 (MQ=11)
                                                     aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                          <  1:13376500/101‑1 (MQ=11)
                                                      nnncGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                         <  2:12147090/98‑1 (MQ=11)
                                                      aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatatnnn                         >  2:14200387/1‑98 (MQ=11)
                                                      aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                         >  2:94273/1‑101 (MQ=33)
                                                       aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                        >  1:5586525/1‑101 (MQ=35)
                                                       aCCGACATATTGATTTCCTCAAGAAGGATATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                        <  1:13682067/101‑1 (MQ=11)
                                                       aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                        <  2:6176156/101‑1 (MQ=11)
                                                       aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                        <  2:5724343/101‑1 (MQ=11)
                                                        ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                       >  2:6380918/1‑101 (MQ=35)
                                                        ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                       >  1:83141/1‑101 (MQ=11)
                                                         nnncATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                      <  2:12681847/98‑1 (MQ=38)
                                                         nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaaaa                      <  2:14720517/98‑3 (MQ=11)
                                                         cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGTAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                      >  1:4789765/1‑101 (MQ=11)
                                                         cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGannn                      >  2:11176917/1‑98 (MQ=11)
                                                         cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                      >  1:9908078/1‑101 (MQ=11)
                                                          gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaannn                     >  2:9043237/1‑98 (MQ=25)
                                                          gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGCACAg                     >  2:4937974/1‑101 (MQ=25)
                                                          gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                     <  2:5416161/101‑1 (MQ=38)
                                                          gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                     >  1:7449198/1‑101 (MQ=14)
                                                          gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAAGAATACCGGGGGTGCTGAGTGAATATATCGaannn                     >  2:10717468/1‑98 (MQ=2)
                                                          gACATATTGATTTCCTCAAGAAGGACATCAGTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                     <  1:12780147/101‑1 (MQ=14)
                                                           aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATAGATCGAACAGt                    >  1:13112047/1‑101 (MQ=255)
                                                           aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                    <  1:8951572/101‑1 (MQ=14)
                                                            cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                   <  1:1744323/101‑1 (MQ=40)
                                                            cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCGGAGTGAATATATCGAACannn                   >  2:13182408/1‑98 (MQ=255)
                                                            cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGGTGCTGAGTGAATATATCGAACannn                   >  2:12579737/1‑98 (MQ=11)
                                                             atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                  >  1:5118609/4‑101 (MQ=14)
                                                             atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                  <  2:6846006/98‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAACGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:12590521/99‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  2:4301774/99‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:10350385/99‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  2:2433323/99‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:13539510/99‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:1819050/99‑1 (MQ=14)
                                                              tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  2:6437919/99‑1 (MQ=14)
                                                              tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGtnnn                 >  2:7390025/1‑98 (MQ=255)
                                                              nnntGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  2:11852643/98‑1 (MQ=255)
                                                               gatGATTTCCTCAAACGGTATATCAACAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  2:872337/99‑1 (MQ=14)
                                                               gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:5451696/101‑1 (MQ=255)
                                                               gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:9727896/101‑1 (MQ=255)
                                                               gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:6894158/101‑1 (MQ=255)
                                                               attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:12409723/100‑1 (MQ=14)
                                                               attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:4504266/100‑1 (MQ=14)
                                                               attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:10969857/100‑1 (MQ=14)
                                                               attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:14337152/100‑1 (MQ=14)
                                                                ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt               <  1:6192252/101‑1 (MQ=255)
                                                                ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt               <  1:10560909/101‑1 (MQ=14)
                                                                nnnaTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt               <  2:13540582/98‑1 (MQ=255)
                                                                 tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGta              >  1:5836247/1‑100 (MQ=255)
                                                                 tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgggt              >  2:1996240/1‑99 (MQ=14)
                                                                 tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              >  1:7140454/1‑101 (MQ=14)
                                                                  gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa             <  2:1388016/101‑1 (MQ=17)
                                                                   aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGACCAGTCAGGTTaa            <  2:4648692/101‑1 (MQ=255)
                                                                   aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  1:12663679/1‑101 (MQ=255)
                                                                   aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATAGCGAACAGTCAGGTTaa            >  2:5976721/1‑101 (MQ=255)
                                                                   aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGTTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            <  1:14464767/101‑1 (MQ=17)
                                                                   aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtnnn            >  2:14848106/1‑98 (MQ=14)
                                                                    tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc           >  2:953008/1‑101 (MQ=255)
                                                                    tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc           <  1:8516342/101‑1 (MQ=255)
                                                                    nttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc           <  1:2267351/100‑1 (MQ=255)
                                                                      tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCTCTAGATGAACTCATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg         <  1:13231470/101‑1 (MQ=255)
                                                                      tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg         <  1:9375790/101‑1 (MQ=255)
                                                                      tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg         <  1:6653218/101‑1 (MQ=255)
                                                                      tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTACTGAGTGAATATATCGAACAGTCAGGTTAACAg         <  1:11944049/101‑1 (MQ=255)
                                                                       nnngAAACGGGATATCATCAAAGACATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg        <  2:8445636/98‑1 (MQ=255)
                                                                       ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg        <  1:4413040/101‑1 (MQ=255)
                                                                       ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGCGAATATATCGAACAGTCAGGGTAACAgg        >  1:10155597/1‑101 (MQ=255)
                                                                        cTGAAACGGGATGTCATCAAAGCCATGAACAAAGCAGCCGCGCTGTATGAAATGATATCGGGGTTGCTGAGTGAATATATCGAAAAGTCAGGTTAACannn       >  2:7720875/1‑98 (MQ=17)
                                                                        cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc       <  1:5660218/101‑1 (MQ=255)
                                                                         tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGGATATATCGAACAGTCAGGTTAACAgnnn      >  2:10413439/1‑98 (MQ=255)
                                                                          nnnaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg     <  2:10097770/98‑1 (MQ=255)
                                                                          gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg     <  1:8138477/101‑1 (MQ=255)
                                                                          gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg     >  1:2504775/1‑101 (MQ=255)
                                                                          gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg     >  1:9135215/1‑101 (MQ=255)
                                                                          gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg     <  1:13447/101‑1 (MQ=255)
                                                                           aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc    <  1:496726/101‑1 (MQ=255)
                                                                           aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc    <  1:3388808/101‑1 (MQ=255)
                                                                           aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc    >  1:11353838/1‑101 (MQ=255)
                                                                           aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGTCAGGTTAACAGGCTGc    >  2:4104079/1‑101 (MQ=255)
                                                                           aaaCGGGATATCATCAAAGCCATGAACAAAACAGCCGCGCTGGATGAACTGAGACCGGGGTTGCTTAGTGAAAATATCAAACAGTCAGGATAACAGGCTgg    >  2:3126275/1‑100 (MQ=17)
                                                                            aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg   <  1:2072358/101‑1 (MQ=255)
                                                                             |                                                                                                    
AGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGG  >  NC_000913/580435‑580612

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