Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5240CT66.5% 130.4 / 206.9 219I156I (ATC→ATTnohBbacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base T (68/77);  total (105/114)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 304 total aligned reads displayed.

CCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATT  >  NC_000913/580452‑580616
                                                                        |                                                                                            
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                                          <  1:524023/89‑1 (MQ=0)
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                                          >  1:4534568/5‑93 (MQ=0)
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                                          >  2:580076/5‑93 (MQ=0)
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                                          >  2:2737074/5‑93 (MQ=0)
          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                                <  2:1367907/90‑1 (MQ=0)
          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                                <  2:1521116/90‑1 (MQ=0)
          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                                >  2:2287180/4‑93 (MQ=0)
          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                                <  2:4924481/90‑1 (MQ=0)
                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                                     <  2:2552831/89‑1 (MQ=0)
                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                                     >  1:819345/5‑93 (MQ=0)
                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                                     <  1:5162767/89‑1 (MQ=0)
                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                                     <  1:2916303/89‑1 (MQ=0)
                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGATATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                                    <  1:1550967/90‑1 (MQ=0)
                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                                    >  2:4629390/4‑93 (MQ=0)
                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                                    <  2:586894/90‑1 (MQ=0)
                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                                    <  1:2041824/90‑1 (MQ=0)
                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAACTGCGCTGGATGAAATAATACCg                                                    >  1:3483239/4‑93 (MQ=0)
                            gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                              >  1:1000087/4‑93 (MQ=0)
                            gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                              >  1:1378002/4‑93 (MQ=0)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                <  2:1710688/93‑1 (MQ=11)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                <  2:5382923/93‑1 (MQ=11)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                >  2:1187199/1‑93 (MQ=11)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                >  2:3842000/1‑93 (MQ=2)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                >  1:92400/1‑93 (MQ=2)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                >  1:5806432/1‑93 (MQ=2)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                >  1:602687/1‑93 (MQ=2)
                                          aCATATTGATTGCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                <  2:3368210/93‑1 (MQ=11)
                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                               >  1:1337043/1‑93 (MQ=11)
                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                               <  1:5619277/93‑1 (MQ=11)
                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                               >  1:2158758/1‑93 (MQ=11)
                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                               >  2:4312260/1‑93 (MQ=2)
                                            atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATca                              >  1:2770265/4‑92 (MQ=2)
                                            atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                              <  1:651362/90‑1 (MQ=11)
                                            aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                              <  2:163070/93‑1 (MQ=33)
                                            aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                              >  2:632370/1‑93 (MQ=33)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                             >  1:6087248/3‑93 (MQ=11)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                             >  2:4267047/3‑93 (MQ=11)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                             >  1:4735511/3‑93 (MQ=2)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGCGAATATATCGa                             >  1:2589016/3‑93 (MQ=11)
                                              gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                            <  2:1915356/93‑1 (MQ=35)
                                               ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                           <  1:4246980/93‑1 (MQ=37)
                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                           <  2:5543956/93‑1 (MQ=11)
                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                           <  2:4452808/93‑1 (MQ=11)
                                                tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                          <  1:4878588/93‑1 (MQ=38)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                          <  2:2987389/93‑1 (MQ=14)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                          <  2:2390409/93‑1 (MQ=14)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                          <  1:4595290/93‑1 (MQ=14)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAATGAATATATCGAACa                          <  1:4136790/93‑1 (MQ=14)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGAAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGCGAATATGTCGAACa                          >  1:387439/1‑93 (MQ=2)
                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                         <  1:3222314/93‑1 (MQ=39)
                                                 gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                         >  1:5132580/1‑93 (MQ=14)
                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGt                        <  1:1053421/93‑1 (MQ=40)
                                                   tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                       <  2:4782948/93‑1 (MQ=255)
                                                   tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                       >  1:4380265/1‑93 (MQ=255)
                                                    ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                      <  2:5327633/93‑1 (MQ=255)
                                                     tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                     <  2:4189965/93‑1 (MQ=255)
                                                      tgtGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCACTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                    <  1:1804990/91‑1 (MQ=255)
                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt                  >  1:3396879/1‑93 (MQ=255)
                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt                  >  1:3449895/1‑93 (MQ=255)
                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt                  >  1:1510486/1‑93 (MQ=255)
                                                          aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa                <  1:4703775/93‑1 (MQ=255)
                                                          aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGATTGCTGAGTGAATATATCGAACAGTCAGGTTaa                <  2:2877187/93‑1 (MQ=255)
                                                           aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc               >  1:2014941/1‑93 (MQ=255)
                                                           aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc               <  2:5177873/93‑1 (MQ=255)
                                                           aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc               >  1:1673271/1‑93 (MQ=255)
                                                            aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa              >  1:3131321/1‑93 (MQ=255)
                                                            aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa              >  1:549637/1‑93 (MQ=255)
                                                             cGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg             >  2:3172730/1‑93 (MQ=255)
                                                                gATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCt          >  1:2319939/1‑93 (MQ=255)
                                                                gATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCt          >  2:4113838/1‑93 (MQ=255)
                                                                 acatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTg         >  1:3988105/3‑93 (MQ=17)
                                                                 acatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTg         >  1:2710291/3‑93 (MQ=17)
                                                                  catcattAAAGCCATGAACAGAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGc        <  1:2272603/92‑1 (MQ=17)
                                                                  catcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGc        <  2:1004016/92‑1 (MQ=17)
                                                                  catcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGc        <  2:4513005/92‑1 (MQ=17)
                                                                  catcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGc        >  2:6068628/2‑93 (MQ=17)
                                                                  catcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGc        >  1:3580957/2‑93 (MQ=17)
                                                                   atcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGca       >  2:3299427/1‑92 (MQ=17)
                                                                   atcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGca       <  2:5658260/93‑2 (MQ=17)
                                                                   atcatcAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg       >  1:309620/1‑93 (MQ=255)
                                                                   atcatcAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg       <  2:4075876/93‑1 (MQ=255)
                                                                   atcatcAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg       >  1:5454729/1‑93 (MQ=255)
                                                                   atcatcAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGCTAACAGGCTGCg       >  1:4316617/1‑93 (MQ=255)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTggag      >  2:737309/1‑90 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  2:3659637/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      <  1:1946554/93‑3 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  1:1460662/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  2:600490/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  1:2071994/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  1:732156/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  2:3564029/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  2:2847640/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  2:2767738/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  2:2590247/1‑91 (MQ=17)
                                                                    tcattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag      >  1:370143/1‑91 (MQ=17)
                                                                    tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCgg      >  1:5391519/1‑93 (MQ=255)
                                                                     cattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGc     <  2:3626283/93‑1 (MQ=17)
                                                                     cattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGc     >  2:2038680/1‑93 (MQ=17)
                                                                     cattaAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGc     >  1:4395706/1‑93 (MQ=17)
                                                                       tcaAAGCCATGAACAACGCAGCCGCGCTGGATGAACTGATACCCAGGCTGATGAGTGAGTATATAGAACAGTCAGGGTAACATGCTGCGGCAt   >  2:83538/1‑93 (MQ=18)
                                                                        |                                                                                            
CCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATT  >  NC_000913/580452‑580616

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: