Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,492,8430CT28.9% 60.2 / 10.6 45F92F (TTC→TTTtsgAputative transporter TsgA
Reads supporting (aligned to +/- strand):  ref base C (5/27);  new base T (13/0);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCT  >  NC_000913/3492599‑3493092
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                   
ggATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCGGTTTTTTGTTCAGCcac                                                                                                                                                                                                                                                               >  2:411827/1‑251 (MQ=255)
               ccTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTCCGGTTTGTTGTTTAGCCCCAGCCTGGCGCTGTtc                                                                                                                                                                                                                                               >  1:419671/1‑251 (MQ=255)
                 tACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGCAATATCGCCGATTATTTCACTCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGGTGTTTAGCCCCACCCCGGCGCTGTtctt                                                                                                                                                                                                                                             >  1:200586/1‑250 (MQ=255)
                  aCGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATCGTGATGGGCAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGTCGGGTTTATGTTTAGCCCCAGCCTGgcggcgttcttc                                                                                                                                                                                                                                             >  1:230145/1‑242 (MQ=255)
                      actgactgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGTTGTTTAGCCCCAACCTGGCGCTGGTCTcggcgg                                                                                                                                                                                                                                        >  1:136434/1‑251 (MQ=255)
                                                ccGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGc                                                                                                                                                                                                           <  1:128407/251‑1 (MQ=255)
                                                     aTGGTGATGGGAAATATCGCCGATTATTGCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTGGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCAt                                                                                                                                                                                                      <  2:225120/251‑1 (MQ=255)
                                                       ggTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGGCACAACCTGGCGCTGGTCTCGGGGGGGGTGTTCCTTTTCGGGGGGGGCAGCGGccac                                                                                                                                                                                                     >  1:416034/1‑247 (MQ=255)
                                                        gTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCAc                                                                                                                                                                                                   <  1:149251/251‑1 (MQ=255)
                                                         tGATGGGACACATCGCCGATTATTTCAATCTGCCTGTTTCCAGAATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGCTGGTGCTGGCGGTTGCCGGGTTTATTTTTAGCCACACCCCGGGCTTCTTCTCGGGGGCGATGTTATTTCTCCGGGGGGGTCggcggcgtccc                                                                                                                                                                                                   >  2:182421/1‑246 (MQ=255)
                                                                  atatCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTCGCCCCAGCCCGGGGGTGTTTTCGGGGGCGGCTGTGCTTCTCGGGGTGGGCAGCAGCATCAccctctagctt                                                                                                                                                                                         >  2:113135/1‑242 (MQ=255)
                                                                       gCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCCGGTTGGTGGTCTGGCACAGCCCTGCGGGGTTTTCCGCGGCGGATGTTCTTTCTCGGGGGGGTAGCGGcctccacctcgcggattggtc                                                                                                                                                                                    >  1:4875/1‑231 (MQ=255)
                                                                           attattTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTc                                                                                                                                                                               <  2:362360/251‑1 (MQ=255)
                                                                             tattTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTcc                                                                                                                                                                              <  1:44551/250‑1 (MQ=255)
                                                                                ttCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGGTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTgg                                                                                                                                                                           <  2:175497/250‑1 (MQ=255)
                                                                                  cAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTAGCCCCAGCCTCGCGCCGTTCTCGGCGGCGATGTTCATTTTCGGGGTGGGGGGCGGCATCCACATGTCGATTGGTACCTTTCTggga                                                                                                                                                                         >  1:350429/1‑248 (MQ=255)
                                                                                        gCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAacac                                                                                                                                                                     <  1:17094/248‑1 (MQ=255)
                                                                                                ccAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTTAGCCAAAGCCAGGCGCTGGTTTCGGCGGGGGTGGTTCTTTTCGGGGGGGGTAGCGGCCTCACCAAGTCGGTTTGTACATTCCTGGGAAAACAAATGTAtta                                                                                                                                                           >  1:205826/1‑248 (MQ=255)
                                                                                                   gTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACTTTCCTGGTAACACAAAGGTATGAAggg                                                                                                                                                       <  1:366924/251‑1 (MQ=255)
                                                                                                       gAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTc                                                                                                                                                  <  1:411827/251‑1 (MQ=255)
                                                                                                           aaTACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGGTGTTTCGCCCCACCCCGGCGCTGTTTTCGGCGGCGGTGTTCCTTTTCGGGGGGGGCAGCGGGATCACCCTGTCGGTTTGTTCATTCCTGGTAACACAAATATATGAAGGGCGTCAgcg                                                                                                                                              >  1:339247/1‑251 (MQ=255)
                                                                                                                              aCGCCGGCATTTTAATCTCTAGCTGCCTCAACGCCTGGCTTAGGGAAAGCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCGGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTGCATTCTCGGGGTGGTCGGCGGCATCACCATGTCGATGGGTACATTCCTGGGAACGCAAATGTAGGAAGGGCGTCAGCGCGTTTCCCGCCTGTTAttt                                                                                                                           <  2:61318/251‑1 (MQ=255)
                                                                                                                                              ctctATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTCGCCGTTTTTTGTTTAGCCCCAACCTCGGGCTGTTTTCCGCGGCGGTGTTTATTTTCGGGGGGGGCACCGGCGTCACCCTTTCCCTTTGTACATTTCCTGTGACACAAACATATTAAGGGGGGCAGCGCGGGTCCCGCCTGTTAtttttccaatccctctttc                                                                                                           >  1:280325/1‑235 (MQ=255)
                                                                                                                                                      tccccAACGCCTGGCTGAGGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTGCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGTTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATg                                                                                                    <  1:305231/244‑1 (MQ=255)
                                                                                                                                                       ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTACTCCTGATGGTGCTGGCGGTTGCCGGGTTTATGGTCAGCCCCAGCCCCGCGCGGTTCTCCGCCGGGGTGTTTATTCTCCGGGGGGGCGGGGGGCGTCCCCAGGCGGTTGGTAGCTTTCTTGTAACACAAACGTATTAAGGGGGTCAGAGCGGGTCCCGCCCGTTATTTATACACCCCTTCTTCATCATGGCTg                                                                                                >  2:98587/1‑251 (MQ=255)
                                                                                                                                                       ccTCAACGCCTGGCTGAGGGAAATCGTCCCGTGGAAAACGCATTGAGGTTTTGTCTTTCTCCTGATGGTGCTGGGGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCGCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTg                                                                                                <  2:137734/251‑1 (MQ=255)
                                                                                                                                                            ccGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGgatg                                                                                           <  1:101089/250‑1 (MQ=255)
                                                                                                                                                              gccggtcTGATGGAAACCGTCCCGTTGAAACGGTGGTTAGGTTTTGTTTTTTTCTTGATGGTGCTGGGGGTTGCCGTTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATGGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGgatgat                                                                                        <  2:264055/245‑1 (MQ=255)
                                                                                                                                                               ccTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAcc                                                                                                                        <  1:129040/221‑1 (MQ=255)
                                                                                                                                                                 ttgctgtgtgAAATCGTCCAGCTGAAAACGCAGTAGCGTTTATGCTTTCTCCTGTGGGTGCTGGCGGGTGCCGGTTTGAGGTGCAGCCACACCCTGGCGCGGTTCTCGGCGGCGATGTTCTGTCTCGGGGTTGTCAGCGGCATCACCATGTCGATTGGTACATTCCGGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGAtt                                                                                       <  2:313816/240‑1 (MQ=255)
                                                                                                                                                                  ggCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTGCGGTTTGATGTTCAGCCCCAGCCCGGCGCTGTTCTCGGCGGCGATGTTCATTTTTGGGGGGGGCAGCGGCATCACCCTGTTGATTGGTACATTCCTGGGGACACAAAAGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGGTATTTTCCGACTCCTTCTTTAGTATGGCTGGGGTGGttttt                                                                                    >  1:372887/1‑250 (MQ=255)
                                                                                                                                                                      gATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTAGCCACAGCCCGGCGCTGTTTTCGGCGGGGGTGGTCATTTTCGGGGGGGGCAGCGGCCTCACCCTGTTGATTTGTACCTTCCCGGGAACACAAAAGTATGAAAGGCGGTCACGCGGTTTCCCCCCGTTTTTTTTCCACTCCTTCTTCCGTAtgggtggggttgtttttccc                                                                                  >  1:127414/1‑238 (MQ=255)
                                                                                                                                                                       aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCCGGTTGGTTTTTAGCCCCCACCTCGCGCTGTTTTCGGGGGGGGTGTTTATTTTCGGGGGGGGCAGCGGGATCCCCATGGCGAGTGGGAAATTTCTTGTAACACAAATGTATGAAGGGGGGCCGCGCGGGTCCCGCCTGGTATTTACCGACGCCTTCTTCAGTAATGGTGGGATGAGTTTTccac                                                                                 >  1:12426/1‑248 (MQ=255)
                                                                                                                                                                                   ccGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTAACCACCACCCGGGGCCGTTCTCCGGGGCGGTGTTCATTCTCGGGGCGGTCAGCGGCACCACCATGTCGGTTTGTACAATTCTGGGCACCCAAATGAATGAAGGGGGTCAACGCGGTTCCCGCCTGTTTTTTACCGACTCCTTTCTTAGTTTGGGTGGGGTGGTTTTTCCAATggaccccgcggt                                                                   >  1:299794/1‑240 (MQ=255)
                                                                                                                                                                                      gtGAAAACGCAGTTACGTGTTGGCTGTCTCCTGATGGGTCTGGCGTTTGCCGTTTTGAGGTTCAGCCACAGCCGGGCGCTGTTCTCGGCGGCGAGTTTCATCCTCGGGGGGGTCAGCGGCATCACCATTTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCt                                                               <  2:419671/250‑1 (MQ=255)
                                                                                                                                                                                         aaaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACt                                                            <  1:159452/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTgg                                         <  1:160672/251‑1 (MQ=255)
                                                                                                                                                                                                                           gtcTGGTGGTTGCCGGTTTGTTGTTCGGCCACGGCCTGGCGCTGTTCTGGGCTGCGATTTTCGTTCTCGGGGTGGTCGGCGGCATCACCAGGGCGATTGGGTGTTTCCTGGTAAGACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGGTATTTACCTACGCTTTCTTCAGTATGGCTGGGATGAGTTTCCCAAGGATCGCCGCGTGTCGACGTTCGCGCAGCAGTGAGTGGTAGTGTGTTTATGcc                          <  2:220726/249‑1 (MQ=255)
                                                                                                                                                                                                                           ggCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCGCGGCGGCGATGTTCATTCTCGGGGTGGGCAGCGGCATCACCATGGCGATTGGTACAGTCCTGGTAACACAAATGTATGAAGGGCGTCGGCGCGGTTCCCGCCTGTGATTTACCGACTCCTTCTTCGGTATGGCTGGGATGAGTTTCCCAATGACCGCCGGTATTCGACTGGCGCGCAGCATTGAGTGGGACTGGGTTTATGcc                          <  1:422089/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gggTTGCGGGTTTGATGTTCAGCCACGGCCTGGCGCTTTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATc                    <  2:111803/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGGTCATTCCCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACGCAAATGTATGAAGGGCGTCAGCGCGGTGCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATc                    <  1:447597/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCg                   <  2:136434/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCTTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGgt            <  1:98587/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gTTTGAGGTTCAGCCAACGCCGGTCGCTGTTCGCGGCTGCGATGTCCATGCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTAGGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGgtt           <  2:16660/251‑2 (MQ=255)
                                                                                                                                                                                                                                                ggATGTCCAGCCACTGCCTGGGGGTGTTCCCGGTGGGGGTGTTCAGTCCCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACAGTCCTGGTAACACAAAGGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAATTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTGCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATTGCTGCATCGGGCTGGTGTAt        <  2:280325/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTAAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGTGCTGGTGTATGTCGCt  <  1:144625/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                   
GGATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCT  >  NC_000913/3492599‑3493092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: