Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913412,1240AG64.2% 268.0 / inf 313S952P (TCC→CCC) sbcCATP‑dependent dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (62/50);  new base G (119/82);  total (181/132)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 339 total aligned reads displayed.

AAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAG  >  NC_000913/412032‑412215
                                                                                            |                                                                                           
 aGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGaa                                                                                            <  1:5897293/93‑1 (MQ=255)
  ggAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGaaa                                                                                           <  2:564606/93‑1 (MQ=255)
   gAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTCACGAGGAAACTTTGGCCGCCGGAAAg                                                                                          >  1:4698271/1‑93 (MQ=255)
     aCAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAggg                                                                                        >  2:6882236/1‑93 (MQ=255)
     aCAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAggg                                                                                        <  2:3583820/93‑1 (MQ=255)
     aCAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAggg                                                                                        >  1:649618/1‑93 (MQ=255)
      cAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGt                                                                                       >  2:2681281/1‑93 (MQ=255)
       aGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTa                                                                                      >  1:367703/1‑93 (MQ=255)
       aGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTa                                                                                      <  1:3913781/93‑1 (MQ=255)
           aGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACggg                                                                                  >  2:6991660/1‑93 (MQ=255)
           aGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACggg                                                                                  <  2:573980/93‑1 (MQ=255)
           aGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACggg                                                                                  <  2:8289054/93‑1 (MQ=255)
           aGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACggg                                                                                  <  2:6225722/93‑1 (MQ=255)
            gTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGt                                                                                 <  1:8813507/93‑1 (MQ=255)
            gTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGt                                                                                 >  1:7200571/1‑93 (MQ=255)
                aTACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGCAACTTTCGCCGCCGGGAAGGGTACGGGTATcg                                                                             >  2:2989647/1‑93 (MQ=255)
                 tACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATcgc                                                                            <  2:4077724/93‑1 (MQ=255)
                  tcGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGGCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATcgcg                                                                           <  1:5446238/92‑1 (MQ=255)
                  aCGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATcgcg                                                                           <  1:1961744/93‑1 (MQ=255)
                   cGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGt                                                                          >  1:7917814/1‑93 (MQ=255)
                    gtgtTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTa                                                                         >  2:7901347/1‑93 (MQ=255)
                     tgtTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTAc                                                                        >  1:5302170/1‑93 (MQ=255)
                      gtTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTAcc                                                                       >  1:5866608/1‑93 (MQ=255)
                       ttttATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCg                                                                      >  1:2451830/1‑93 (MQ=255)
                       ttttATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCg                                                                      >  2:6449743/1‑93 (MQ=255)
                           aTGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATc                                                                  >  1:2917634/1‑93 (MQ=255)
                               cTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATctgcc                                                              >  1:3845494/1‑93 (MQ=255)
                                tGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATctgcct                                                             >  1:7680196/1‑93 (MQ=255)
                                 gACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACGAACGAGGAAACCACCGCCGCCGGGAAGAGCACGGGAAGCGCGCACCGCAgccgccgg                                                            >  1:6708679/1‑91 (MQ=255)
                                  aCCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATctgcctgc                                                           >  1:2741729/1‑93 (MQ=255)
                                    cAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCa                                                         >  2:5135083/1‑93 (MQ=255)
                                     aGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAg                                                        >  1:2185519/1‑93 (MQ=255)
                                        tCCGAAAGCGCCAGCGCCAGCTCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTa                                                     <  2:381757/93‑1 (MQ=255)
                                        tCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTa                                                     >  2:3261240/1‑93 (MQ=255)
                                        tCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTa                                                     <  2:5915808/93‑1 (MQ=255)
                                          cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATc                                                   >  1:1116204/1‑93 (MQ=255)
                                          cGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGTTATc                                                   >  1:765789/1‑93 (MQ=255)
                                           gAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATca                                                  >  2:261583/1‑93 (MQ=255)
                                             aaGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATcaac                                                <  2:3335947/93‑1 (MQ=255)
                                             aaGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATcaac                                                >  2:4825266/1‑93 (MQ=255)
                                                cgccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCaacaac                                             <  2:2439368/93‑1 (MQ=255)
                                                cgccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCaacaac                                             <  1:754523/93‑1 (MQ=255)
                                                 gccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAAcc                                            <  2:8386057/93‑1 (MQ=255)
                                                  ccagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCt                                           >  1:6212843/1‑93 (MQ=255)
                                                   cagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTc                                          <  1:1507615/93‑1 (MQ=255)
                                                   cagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTc                                          <  1:2017833/93‑1 (MQ=255)
                                                   cagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTc                                          <  2:1409572/93‑1 (MQ=255)
                                                   cagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTc                                          <  1:7001066/93‑1 (MQ=255)
                                                   cagcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTc                                          >  1:1020757/1‑93 (MQ=255)
                                                    agcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCg                                         >  2:4312522/1‑93 (MQ=255)
                                                     gcgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGa                                        >  1:707906/1‑93 (MQ=255)
                                                      cgccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGAc                                       <  2:3448907/93‑1 (MQ=255)
                                                       gccagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATAGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACt                                      <  1:979088/93‑1 (MQ=255)
                                                       gccagcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACt                                      >  1:6301267/1‑93 (MQ=255)
                                                         cagcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGTGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTc                                    >  1:6359833/1‑93 (MQ=255)
                                                         cagcgccaaACTAACGAGGACGCTTTCGCCGCCGGAAAGGGTACGGGTATCGCTTACCGCATCTGCCTGCCATGTATCAACAACCTCGACTTc                                    <  2:1222680/93‑1 (MQ=255)
                                                           gcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCa                                  >  1:686427/1‑93 (MQ=255)
                                                           gcgccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCa                                  <  2:5304812/93‑1 (MQ=255)
                                                           gcgccagACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCa                                  <  1:4514261/93‑1 (MQ=255)
                                                             gccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAgc                                >  1:1779203/1‑93 (MQ=255)
                                                             gccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAgc                                <  2:4141483/93‑1 (MQ=255)
                                                              ccagACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAgcg                               >  2:4920114/1‑93 (MQ=255)
                                                                agACTAACGAGGAAACTTTCGCCGCCGGGAAGGGGACGGGTATCGCGCACCGCATCTTCCTGCCAGGTATCAACAACCTCTACTTCCAGCGcc                             >  2:1892952/1‑93 (MQ=255)
                                                                 gACTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCt                            >  2:3668432/1‑93 (MQ=255)
                                                                  aCTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTc                           >  1:3297616/1‑93 (MQ=255)
                                                                  aCTAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTc                           >  1:1234456/1‑93 (MQ=255)
                                                                  aCTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTc                           <  2:8139136/93‑1 (MQ=255)
                                                                  aCTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTc                           <  2:7210583/93‑1 (MQ=255)
                                                                    tAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGc                         <  1:4958542/93‑1 (MQ=255)
                                                                    tAACGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGc                         <  1:928964/93‑1 (MQ=255)
                                                                    tAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGc                         >  2:172781/1‑93 (MQ=255)
                                                                    tAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGc                         >  1:4252812/1‑93 (MQ=255)
                                                                    tAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGc                         <  1:29533/93‑1 (MQ=255)
                                                                      aCGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTg                       <  1:4096621/93‑1 (MQ=255)
                                                                      aCGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTg                       >  1:508770/1‑93 (MQ=255)
                                                                       cGAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCATGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTgg                      <  1:7319922/93‑1 (MQ=255)
                                                                        gAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGc                     >  1:1587194/1‑93 (MQ=255)
                                                                        gAGGAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGc                     >  1:3444811/1‑93 (MQ=255)
                                                                         aGGAAACTTTCGCCGCCGGAAAGGGTACGGGGATAGCGTACCGCATCTGCCTGCCACGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCt                    >  2:5826649/1‑93 (MQ=255)
                                                                          ggAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCtt                   >  1:147403/1‑93 (MQ=255)
                                                                          ggAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCtt                   >  1:3484781/1‑93 (MQ=255)
                                                                          ggAAACTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCtt                   <  2:1022132/93‑1 (MQ=255)
                                                                          ggAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCtt                   <  1:2780710/93‑1 (MQ=255)
                                                                             aaCTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTgc                <  1:7723593/93‑1 (MQ=255)
                                                                             aaCTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTgc                >  1:8213522/1‑93 (MQ=255)
                                                                              aCTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTgcg               <  1:2124786/93‑1 (MQ=255)
                                                                              aCTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTgcg               >  2:7525298/1‑93 (MQ=255)
                                                                               cTTTCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTgcgc              <  1:5362228/93‑1 (MQ=255)
                                                                                tttCGCCGCCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCt             >  2:8175888/1‑93 (MQ=255)
                                                                                   cgccgcCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCGGCGCCTCGCTGGCTTTGCGCTGTa          >  2:4637276/1‑93 (MQ=255)
                                                                                   cgccgcCGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTa          >  2:7467443/1‑93 (MQ=255)
                                                                                    gccgccGGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTaa         >  2:158851/1‑93 (MQ=255)
                                                                                       gccgGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAg      <  2:7120112/93‑1 (MQ=255)
                                                                                       gccgGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAg      >  2:6624350/1‑93 (MQ=255)
                                                                                       gccgGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAg      >  1:5962542/1‑93 (MQ=255)
                                                                                        ccgGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGa     >  1:4592210/1‑93 (MQ=255)
                                                                                         cgGGAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGAt    <  1:8520024/93‑1 (MQ=255)
                                                                                         cgGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGAt    <  1:261292/93‑1 (MQ=255)
                                                                                          ggAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATa   <  1:1724840/93‑1 (MQ=255)
                                                                                           gAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCCCTGGCTTTGCGCTGTAACAGATAg  >  2:6746929/1‑93 (MQ=255)
                                                                                            |                                                                                           
AAGGAACAGCGAGTCAATACGTGTTTTATGGCTGACCAGATCCGAAAGCGCCAGCGCCAGCGCCAGACTAACGAGGAAACTTTCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAG  >  NC_000913/412032‑412215

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