Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5140GA47.9% 364.5 / inf 309R153Q (CGG→CAG) ‡nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (79/82);  new base A (67/81);  total (146/163)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 326 total aligned reads displayed.

ATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCA  >  NC_000913/580438‑580614
                                                                            |                                                                                                    
               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                               >  2:6090144/4‑101 (MQ=0)
               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                               <  1:2445343/98‑1 (MQ=0)
                                  gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                            <  1:12507530/96‑1 (MQ=0)
                                   cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                           <  1:7755946/97‑1 (MQ=0)
                                   cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                           >  1:5901003/5‑101 (MQ=0)
                                          gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGa                                    <  2:2756065/98‑1 (MQ=1)
                                                gAAATCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                              <  1:10898400/101‑1 (MQ=11)
                                                gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                              >  1:5071534/1‑101 (MQ=11)
                                                gAAAACCGACATATTGATTTCCTCAAGAAGGACCTCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                              <  1:5189711/101‑1 (MQ=11)
                                                gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                              >  1:6728672/1‑101 (MQ=11)
                                                gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                              <  2:3109388/101‑1 (MQ=11)
                                                 nnnaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                             <  2:8568709/98‑1 (MQ=14)
                                                 aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTGGCTGAGTGAatatatg                             >  2:4185999/1‑100 (MQ=2)
                                                  nnnccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                            <  2:9181481/98‑1 (MQ=17)
                                                  nnnccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                            <  2:8431523/98‑1 (MQ=11)
                                                  aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                            >  2:5199301/1‑101 (MQ=32)
                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGATGGAATATATCg                            >  2:6769479/1‑101 (MQ=11)
                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                            <  1:10068588/101‑1 (MQ=11)
                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                            <  1:13376500/101‑1 (MQ=11)
                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                            >  1:2120361/1‑101 (MQ=11)
                                                   aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                           <  1:9945328/101‑1 (MQ=33)
                                                   aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                           >  2:94273/1‑101 (MQ=33)
                                                   aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGTTGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                           >  1:3368707/1‑101 (MQ=11)
                                                   aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                           >  1:4052180/1‑101 (MQ=11)
                                                   aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                           <  1:1789826/101‑1 (MQ=11)
                                                    aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                          <  2:1797250/101‑1 (MQ=35)
                                                    aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                          <  1:13424876/101‑1 (MQ=35)
                                                    aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATcnnn                          >  2:12264370/1‑98 (MQ=11)
                                                    aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                          >  1:8469272/1‑101 (MQ=11)
                                                     nnnaCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                         <  2:7334232/98‑1 (MQ=11)
                                                     ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                         <  1:9027107/101‑1 (MQ=35)
                                                     ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                         >  1:13963149/1‑101 (MQ=11)
                                                      nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                        <  2:14780816/98‑1 (MQ=11)
                                                      cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                        <  1:9222576/101‑1 (MQ=38)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                        >  1:5826556/1‑101 (MQ=11)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                        <  1:8108679/101‑1 (MQ=11)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                        >  1:4255952/1‑101 (MQ=11)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                        <  1:12713906/101‑1 (MQ=11)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAAAAATACAGGGGTTGATGAGTGAATATATCGaaaa                        >  1:3488136/1‑99 (MQ=2)
                                                       nnnaTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                       <  2:8932342/98‑1 (MQ=255)
                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                       <  1:10299462/101‑1 (MQ=38)
                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGATTATATCGAACAg                       >  2:4464510/1‑101 (MQ=255)
                                                        nnntattGATTTCCTCAAGAAGGCCATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                      <  2:13795528/96‑1 (MQ=14)
                                                        gcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                      <  2:5743167/100‑1 (MQ=14)
                                                        aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                      >  1:344436/1‑101 (MQ=9)
                                                         nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                     <  2:12747280/98‑1 (MQ=255)
                                                         nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                     <  2:10697805/98‑1 (MQ=255)
                                                         cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                     >  1:722159/1‑101 (MQ=40)
                                                         cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACannn                     >  2:10142447/1‑98 (MQ=2)
                                                         cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                     >  1:204931/1‑101 (MQ=9)
                                                          atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATCCCGGGGTTGCTGAGTGAATATATCGAACAGTCa                    <  1:1374015/98‑1 (MQ=14)
                                                          atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                    >  1:5638108/4‑101 (MQ=14)
                                                          aTGTTGATTCCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                    <  1:1477756/101‑1 (MQ=255)
                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAACGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                   <  1:12590521/99‑1 (MQ=14)
                                                            nnngATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATGCCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  2:12650397/98‑1 (MQ=14)
                                                            gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  1:6277127/101‑1 (MQ=255)
                                                            gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  1:9727896/101‑1 (MQ=255)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  1:10495050/100‑1 (MQ=14)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  2:139798/100‑1 (MQ=14)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  1:14337152/100‑1 (MQ=14)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  1:10969857/100‑1 (MQ=14)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  1:12409723/100‑1 (MQ=14)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                  <  2:6366051/100‑1 (MQ=14)
                                                             ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt                 >  2:1438647/1‑101 (MQ=14)
                                                              tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt                <  2:4656229/101‑1 (MQ=14)
                                                               gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa               <  1:265657/101‑1 (MQ=255)
                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAggnnn               >  2:12085215/1‑98 (MQ=14)
                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa               <  2:2333795/101‑1 (MQ=17)
                                                                aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtnnn              >  2:14484031/1‑98 (MQ=255)
                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc             >  2:3699357/1‑101 (MQ=255)
                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc             >  2:953008/1‑101 (MQ=255)
                                                                  ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa            <  2:5544011/101‑1 (MQ=255)
                                                                  ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGATGCTGGGAGAATATATCGAACAGTCAGGTTAAcg            >  2:6411873/1‑100 (MQ=255)
                                                                  ttCCTGAAACGGGATATCACCAAAGGCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGCACAGTCAGGTTAACa            <  1:4737700/101‑1 (MQ=255)
                                                                  ttCCTGAAACAGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa            <  1:13197517/101‑1 (MQ=255)
                                                                   tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTACTGAGTGAATATATCGAACAGTCAGGTTAACAg           <  1:11944049/101‑1 (MQ=255)
                                                                    cttGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAcnnn          >  2:11414412/3‑98 (MQ=255)
                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg          <  1:14556728/101‑1 (MQ=255)
                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGCGAATATATCGAACAGTCAGGGTAACAgg          >  1:10155597/1‑101 (MQ=255)
                                                                    ccTGAAACGGGATATCAACAAAGCCATGAACACAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg          <  1:5452055/101‑1 (MQ=255)
                                                                     cTGAAACGGGATGTCATCAAAGCCATGAACAAAGCAGCCGCGCTGTATGAAATGATATCGGGGTTGCTGAGTGAATATATCGAAAAGTCAGGTTAACannn         >  2:7720875/1‑98 (MQ=17)
                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc         >  1:1159978/1‑101 (MQ=255)
                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc         >  1:12863638/1‑101 (MQ=255)
                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc         <  1:8842837/101‑1 (MQ=255)
                                                                      tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCt        <  1:14176150/101‑1 (MQ=255)
                                                                       nnnaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg       <  2:13480561/98‑1 (MQ=255)
                                                                        acaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGcnnn      >  2:10907706/3‑98 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTGAGGTTAACAGGCTGc      >  1:13583490/1‑101 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc      <  1:9275113/101‑1 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc      >  1:732498/1‑101 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc      <  1:8211750/101‑1 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc      >  1:11353838/1‑101 (MQ=255)
                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAGTATATCGAACAGTCAGGTTAACAGGCtnnn     >  2:11067904/1‑98 (MQ=255)
                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg     <  1:700469/101‑1 (MQ=255)
                                                                          aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGGGTGAATATATCGAACAGTCAGGTTAACAGGcggcgg    >  1:4500089/1‑101 (MQ=255)
                                                                          aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCgg    <  2:3243150/101‑1 (MQ=255)
                                                                           cGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGc   >  1:13508187/1‑101 (MQ=255)
                                                                           cGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGc   >  1:4857725/1‑101 (MQ=255)
                                                                            gggATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCa  <  1:9876758/101‑1 (MQ=255)
                                                                            gggATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCa  >  1:1310759/1‑101 (MQ=255)
                                                                            |                                                                                                    
ATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCA  >  NC_000913/580438‑580614

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