Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5090GC54.4% 226.3 / inf 270L151L (CTG→CTCnohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (59/64);  new base C (67/80);  total (126/144)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 287 total aligned reads displayed.

TCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTG  >  NC_000913/580411‑580609
                                                                                                  |                                                                                                    
tcccgGATCTCAGGTGAAATTGCCAGTATTCTCGCCGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGaa                                                                                                    <  1:4013856/98‑1 (MQ=0)
                                        ttcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTg                                                            >  1:3122214/6‑101 (MQ=0)
                                        ttcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTg                                                            >  2:1092628/6‑101 (MQ=0)
                                         tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                           <  1:9171237/97‑1 (MQ=0)
                                          cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgta                                                          <  2:3461525/98‑3 (MQ=0)
                                                             gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                       <  1:4220264/96‑1 (MQ=0)
                                                              cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                      >  1:5901003/5‑101 (MQ=0)
                                                               ttcccGGAGCTGGAAAACCGACATATTGATTTCATCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCt                                     <  2:3122214/98‑1 (MQ=0)
                                                                           gAAAACGGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                         >  1:14072888/1‑101 (MQ=11)
                                                                           gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                         <  1:4652625/101‑1 (MQ=11)
                                                                           gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                         >  1:10657069/1‑101 (MQ=11)
                                                                           gAAAACCGACATATTGATTTCCTCAAGAAGGACCTCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                         <  1:5189711/101‑1 (MQ=11)
                                                                           gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                         >  1:6728672/1‑101 (MQ=11)
                                                                            nnnaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCCGAGTGAATATATc                        <  2:14295655/98‑1 (MQ=9)
                                                                            aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                        <  2:5475083/101‑1 (MQ=31)
                                                                            aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                        >  1:7309166/1‑101 (MQ=11)
                                                                            aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                        <  1:2703871/101‑1 (MQ=11)
                                                                            aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                        >  1:13932318/1‑101 (MQ=11)
                                                                            aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTGGCTGAGTGAatatatg                        >  2:4185999/1‑100 (MQ=2)
                                                                             nnnccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                       <  2:9181481/98‑1 (MQ=17)
                                                                             aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                       <  1:9243754/101‑1 (MQ=32)
                                                                             aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                       >  1:1797250/1‑101 (MQ=32)
                                                                             aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                       >  2:5199301/1‑101 (MQ=32)
                                                                             aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                       <  1:13376500/101‑1 (MQ=11)
                                                                             aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                       <  1:10068588/101‑1 (MQ=11)
                                                                             aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATCAAATAATACCGGGGTTGCTGAGTGAATATATCg                       <  1:10124699/101‑1 (MQ=11)
                                                                              aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                      <  1:9945328/101‑1 (MQ=33)
                                                                              aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                      >  2:94273/1‑101 (MQ=33)
                                                                               aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATcnnn                     >  2:13478024/1‑98 (MQ=14)
                                                                               aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                     >  1:8687565/1‑101 (MQ=35)
                                                                               aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                     <  1:9690722/101‑1 (MQ=35)
                                                                               aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                     <  2:6176156/101‑1 (MQ=11)
                                                                               aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                     >  1:13692666/1‑101 (MQ=11)
                                                                                ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                    >  2:6380918/1‑101 (MQ=35)
                                                                                ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                    >  1:83141/1‑101 (MQ=11)
                                                                                ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                    >  1:12422827/1‑101 (MQ=11)
                                                                                 nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaaaa                   <  2:14720517/98‑3 (MQ=11)
                                                                                 cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                   <  2:5742046/101‑1 (MQ=38)
                                                                                 cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                   >  2:297457/1‑101 (MQ=38)
                                                                                 cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                   >  1:5826556/1‑101 (MQ=11)
                                                                                 cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                   <  1:12713906/101‑1 (MQ=11)
                                                                                 cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                   <  1:8108679/101‑1 (MQ=11)
                                                                                 cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAAAAATACAGGGGTTGATGAGTGAATATATCGaaaa                   >  1:3488136/1‑99 (MQ=2)
                                                                                  ggcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGCTGCTGAGTGAATATATCGaannn                  >  2:12240249/3‑98 (MQ=2)
                                                                                  gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaannn                  >  2:9043237/1‑98 (MQ=25)
                                                                                  gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                  <  1:10299462/101‑1 (MQ=38)
                                                                                  gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGATTATATCGAACAg                  >  2:4464510/1‑101 (MQ=255)
                                                                                  gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                  >  1:1732193/1‑101 (MQ=14)
                                                                                   aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATAGATCGAACAGt                 >  1:13112047/1‑101 (MQ=255)
                                                                                   aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                 >  1:344436/1‑101 (MQ=9)
                                                                                    nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                <  2:12747280/98‑1 (MQ=255)
                                                                                    cATATTGATTTCGTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACannn                >  2:7616095/1‑98 (MQ=11)
                                                                                    cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                >  1:11544639/1‑101 (MQ=14)
                                                                                    cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                >  1:204931/1‑101 (MQ=9)
                                                                                    cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGGTGCTGAGTGAATATATCGAACannn                >  2:12579737/1‑98 (MQ=11)
                                                                                     atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa               >  1:5638108/4‑101 (MQ=14)
                                                                                     atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa               <  2:6846006/98‑1 (MQ=14)
                                                                                     aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGTcg               >  1:14696985/1‑100 (MQ=255)
                                                                                     aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGCCa               >  1:8994500/1‑101 (MQ=255)
                                                                                      tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg              <  1:13539510/99‑1 (MQ=14)
                                                                                      tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg              <  2:6437919/99‑1 (MQ=14)
                                                                                      tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGtnnn              >  2:7390025/1‑98 (MQ=255)
                                                                                      tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg              <  1:12550944/101‑1 (MQ=255)
                                                                                      nnntGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg              <  2:14461727/98‑1 (MQ=14)
                                                                                      nnntGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg              <  2:13935957/98‑1 (MQ=14)
                                                                                       gatGATTTCCTCAAACGGTATATCAACAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg             <  2:872337/99‑1 (MQ=14)
                                                                                       gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTcnnn             >  2:11359233/1‑98 (MQ=255)
                                                                                       gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg             <  1:6894158/101‑1 (MQ=255)
                                                                                       gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg             <  1:6277127/101‑1 (MQ=255)
                                                                                       attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg             <  1:10495050/100‑1 (MQ=14)
                                                                                       attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg             <  1:14337152/100‑1 (MQ=14)
                                                                                         tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCCGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCCAACAGTCAGGtt           <  1:12603315/101‑1 (MQ=14)
                                                                                         tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgnnn           >  2:13594439/1‑98 (MQ=14)
                                                                                         tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgggt           >  2:1996240/1‑99 (MQ=14)
                                                                                         tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt           >  1:7140454/1‑101 (MQ=14)
                                                                                         tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt           <  1:10948311/101‑1 (MQ=14)
                                                                                         tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt           <  2:4656229/101‑1 (MQ=14)
                                                                                          gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa          >  1:7549804/1‑101 (MQ=17)
                                                                                          gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa          >  1:5003194/1‑101 (MQ=17)
                                                                                           aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtnnn         >  2:11991787/1‑98 (MQ=255)
                                                                                           aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa         <  1:2161757/101‑1 (MQ=255)
                                                                                           aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa         <  1:2502013/101‑1 (MQ=255)
                                                                                           aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATAGCGAACAGTCAGGTTaa         >  2:5976721/1‑101 (MQ=255)
                                                                                           aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGTTTGCTGAGTGAATATATCGAACAGTCAGGTTaa         <  1:14464767/101‑1 (MQ=17)
                                                                                           aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtnnn         >  2:14848106/1‑98 (MQ=14)
                                                                                           aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa         >  1:9830018/1‑101 (MQ=17)
                                                                                            tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc        >  2:953008/1‑101 (MQ=255)
                                                                                            tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc        <  1:8516342/101‑1 (MQ=255)
                                                                                            nnnccTGAAACGGGATATCATCAAAGCCCTGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc        <  2:10657069/98‑1 (MQ=255)
                                                                                             nnncTGAAACGGGATATCATCAAAGCCACGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa       <  2:10390509/98‑1 (MQ=255)
                                                                                              tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCTCTAGATGAACTCATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg      <  1:13231470/101‑1 (MQ=255)
                                                                                              tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg      >  1:11833669/1‑101 (MQ=255)
                                                                                              tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg      <  2:2372361/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)
                                                                                               ccTGAAACGGGATATCAACAAAGCCATGAACACAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg     <  1:5452055/101‑1 (MQ=255)
                                                                                                cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc    >  1:12863638/1‑101 (MQ=255)
                                                                                                  gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg  <  1:4848505/101‑1 (MQ=255)
                                                                                                  gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg  <  1:13447/101‑1 (MQ=255)
                                                                                                  |                                                                                                    
TCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTG  >  NC_000913/580411‑580609

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