Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,751,5370AC33.3% 100.6 / 10.0 55V111G (GTT→GGT) yiaTouter membrane protein YiaT
Reads supporting (aligned to +/- strand):  ref base A (28/8);  new base C (0/18);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.61e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGCC  >  NC_000913/3751310‑3751779
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
acttggcaaacccGCTAGGAGGGGCCTCCCTCTCTGAAAGCCCATAATAATAGTCCCGAAAATTCTCGTAATAATAGAGACCGCCCAGCGAGGGTGTCGCCGACAGACGCCCTATCTGCATTTTGTGGAATCCCGATACCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGTTTTCCCGGTCCCCAACCGCTCGGGGTGATACCGGGCACCCCCCGCCATAGCCGTACTATCa                                                                                                                                                                                                                             <  2:106551/240‑1 (MQ=255)
 aCTTGCCAGACCTCAACGACGGGACTCCCCCTCAAAAACGCCATAGAAATAGCCACGGAATTTCCCGACATAATAGAGACCCCCCAGCGGGGGTGTCGGCGACAGACGCCCTATCTGCATTTTGTGGAATCCCGATAGCCCCCCCCCCCAGCCGTGGCTTTTACCCAGACCGCCCGCAGCTGCAGAGGCTTTCCCGGTCCCCCACCGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCAc                                                                                                                                                                                                                            <  2:34230/251‑1 (MQ=255)
           ccGCTACGACGGGACTCACTCTCTGAAATGCCAAAGAAATAGTCACTAAAATTCTCTTCATAATAGAGACCGCCCAGCGCGGGTGTCAGCGAAAGACGACCTATCTGCATTTTGGGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTTTTATCCAGACCGTCCGCAGCTGCAGAGGCTTCCCCGGTTCCCCAACGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAg                                                                                                                                                                                                                  <  1:204629/251‑1 (MQ=255)
                acgacgGGACTCACTCTCTGAAATGCAATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACCCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGACGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATCGCCGTACTATCACGCTTATCAAGCTGtt                                                                                                                                                                                                             >  1:693904/1‑251 (MQ=255)
                acgacgGGACTCACTCTCTAAAATGCCATAGTAATAGTCCCTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCCCGGTTCCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGt                                                                                                                                                                                                              <  1:520091/250‑1 (MQ=255)
                aacgcggcactcatccgatGAAATCCCATAGGAATAGTCACTGAACTTCTCGTCAGAAAAGAGAACGCCCAGCGGAGGTGTCAGCAACAGTCGACTTCTCTGATTTTTGTGAAATCCCGAGAGCCCCCCCCCCCAGCCGTGCTTGTAATCCAGATCGTCCGCAGCTCCAGAGTCTTCCACGTTCCCCCAACGCTTGTGGTGATACCGTGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGt                                                                                                                                                                                                              <  2:398495/232‑1 (MQ=255)
                  gacgGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                                                                                                                           >  2:694502/1‑251 (MQ=255)
                  gacgGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                                                                                                                           >  2:577027/1‑251 (MQ=255)
                  gacgGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                                                                                                                           >  1:282521/1‑251 (MQ=255)
                     gggACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCGACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATg                                                                                                                                                                                                        >  2:746468/1‑251 (MQ=255)
                       gACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCCGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGc                                                                                                                                                                                                      >  2:425428/1‑251 (MQ=255)
                       gACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGc                                                                                                                                                                                                      >  2:11103/1‑251 (MQ=255)
                           ccctctctGAAATGCCATAGTAATAGTCACTGAAATTCTCTTCATAATAGAGACCCCCCAGCGCAGGTGTCAGCGACAGACGACCTTTCTGCATTTTGTGGAATCCCGATATCTCCCCCCCCCAGCCGTTGCTTTTACCCAGAACGCCCGCAGCTGCAGAGGCTTTCCCGTTCCCCCACCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCAt                                                                                                                                                                                                    <  1:673652/247‑1 (MQ=255)
                           cactcTCTAAAATGCCATAGAAATAGTCACTGAAATTCTCGTCATAATAGAGACCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCCCCCACCCGTTGCTTTTATCCAGAACGCCCGCAGCTGCAGAGGCTTTCACGTTTCCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATaa                                                                                                                                                                                                  <  1:297015/251‑1 (MQ=255)
                                 ctaaaaTGCCATAGTAATAGTCACTAAAATTCTCGTCATAATAGAGACCGCCCAGCGCAGGTGCCAGCAACAGACGCCCTATCTTCTTTTTGTGGAATACCAATACCCCCCCCCCCCAGCCTTTGCTTTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGt                                                                                                                                                                                              <  1:103365/246‑1 (MQ=255)
                                 ctGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGAGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                                                                                                                             >  2:766821/1‑250 (MQ=255)
                                   gAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCAtt                                                                                                                                                                                          >  1:747868/1‑251 (MQ=255)
                                    aaaTGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCAtt                                                                                                                                                                                          >  2:393718/1‑250 (MQ=255)
                                          cATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                                                                                                                    >  1:82884/1‑250 (MQ=255)
                                          cATAGTAATAGTCACTGAAATTCTCGGCATAATAGAGAACGCCCAGCGCAGGAGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGCGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                                                                                                                    >  2:470899/1‑250 (MQ=255)
                                               tAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTa                                                                                                                                                                                         >  2:231438/1‑240 (MQ=255)
                                               aaagaggcatggattttCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCTTTTTGTGGAATCCCGATCGCTCCCCCCCCCAGCGTTTGCTGTTACCCAGACCTTCCGCAGCTGCGGAGGCTTCCACGGTTCCCCAACGCTCGTGGTTATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTTTTGCATGGCATAATCGTCATTa                                                                                                                                                                                         <  1:231438/225‑1 (MQ=255)
                                                          gAAATTCTCGTCAAAATAGAGACCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGACCGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATCCCAGGAAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAg                                                                                                                                                                     <  1:740499/249‑1 (MQ=255)
                                                              ttCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCCGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACAGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAACCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGTCAt                                                                                                                                                               >  2:484396/1‑251 (MQ=255)
                                                              ttCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCAt                                                                                                                                                               >  1:250988/1‑251 (MQ=255)
                                                                     cATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGAt                                                                                                                                                        >  1:702220/1‑251 (MQ=255)
                                                                     aaaaagagagAACCCCCAGCGCTGGTGTCAGGAACAGCCGACCTATCTGCTTTTTGTGAAACCCCGATACCCCCCCCCCCCACCCGTTGCTGTTTTCCAGACCGCCCGCACCTGCAGAGGCTTTCCCGTTCCCCCACCGCTCGTGGTAATACCAGACACCCCCCGCCATAGCTGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGAt                                                                                                                                                        <  2:594753/245‑1 (MQ=255)
                                                                       aaaaggagaACGCCCAGCGCGGGTGTCCGCGACGGACCACCTGTTTGCATTTTGGGGAATCCGGAGGGCTCCCCCCCCCAGCCGTTTCTGTTAACCAAACGGTCCGCAGCGGAAGAGGTTTTCCCGGTCCCCCACCGCTCGTGGTGATACCGGGCAGCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGc                                                                                                                                                      <  2:568951/246‑1 (MQ=255)
                                                                         ataGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCgg                                                                                                                                                    >  1:182922/1‑251 (MQ=255)
                                                                                   cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGATGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTAAAATTCCAGCGGCATCCAGGATGAGGTTAGGCTa                                                                                                                                            >  2:102342/1‑249 (MQ=255)
                                                                                      aGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTc                                                                                                                                       >  2:324641/1‑251 (MQ=255)
                                                                                        cgcgggTGTCGGCGCCAGACGACCTATCGGCTTTTTGTGGAATCCCGATAGCTCCCCCCCCCAGCCGTTGTTGTTACCCATACCGCCCGCAGCTGCAGTGGCTTTCACGTTCCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAt                                                                                                                                     <  1:425428/247‑1 (MQ=255)
                                                                                                    cgacaggcgtcccaccTGCTTTTTGTGGAATCCCGAAACCTCCCCCCCCCGGCGGTTGGGGTTGTCGAGATCCTCCGCAGCTGAAGGGGCTTTCACGGTCCCCCAACCCTCGTTGTATACCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCAcgtt                                                                                                                         <  1:387708/236‑4 (MQ=255)
                                                                                                        aGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCg                                                                                                                      >  1:17028/1‑250 (MQ=255)
                                                                                                           cgcccTATCTGCATTTGGTGGAATACCGATAGCTCCCCCCCCCACCCGTTGCTGTTATCCAGACCGTCCGCAGCGGCAGGGGCTTTCCCGGTTCCCCACCGCTCGTGGTAATCCCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAg                                                                                                                  <  1:185229/248‑1 (MQ=255)
                                                                                                                 atctgctgtttGGGGAATACCGAAAGCCCCCCCACCCAGCCGTTGCTGTTATCAAGAACCTCCGCAGCGGCAGAGGCTTTCACGCTCCCCCAACGCTCGGGGTGAAACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATaa                                                                                                             <  1:314530/242‑1 (MQ=255)
                                                                                                                              ggAATACCGATAGCTCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCCCGGTCCCCCAACGCTCGTGGTGATCCCAGCCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTc                                                                                               <  1:304143/251‑1 (MQ=255)
                                                                                                                                     cGATCGCTCCCCCCCCCAGCCGTTCGTGTTATCCAGACCGCCCGCGGCTGCAGAGGTTTTCACGGTCCCCCACCGCTCGGGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGaa                                                                                        <  1:470899/251‑1 (MQ=255)
                                                                                                                                      gATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAAc                                                                                       >  2:735782/1‑251 (MQ=255)
                                                                                                                                      aaTAGCCCCCCCCCCCAGCCGTTCCGGTTAGCCAGAACCTCCGCCGTTGCAGAGTTTTTCGCGTTTCCCCACCGCTCGGGGTGATACCCGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATGTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAAc                                                                                       <  1:324641/250‑1 (MQ=255)
                                                                                                                                               ccccccccgCCGTTGCTCTTATCCAGACCTTCCGCATCTCAAGCGCTTTCCCCGTTCCCCCAACTCGCGTGGTGATCCCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAgn                                                                              <  2:363763/243‑2 (MQ=255)
                                                                                                                                                 cccccaGCCGTTCCTCTTACCAAGAGCGCCCGCCGCTGCAGAGTCTTTCCCGTTCCCCCAACGCTCGGTGTATTACGCGGCCCCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTAt                                                                            <  2:133115/249‑1 (MQ=255)
                                                                                                                                                  cccccgCCGTTGCTTTTATCCAGACCGCCCGCAGCGGCAGAGGCTTCCCCGTTCCCCAACCGCTGGTGTTACTACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATc                                                                           <  2:403807/246‑1 (MQ=255)
                                                                                                                                                  acccaGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATc                                                                           >  1:34296/1‑251 (MQ=255)
                                                                                                                                                    ccaGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATc                                                                           >  1:403807/1‑249 (MQ=255)
                                                                                                                                                          gTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCAt                                                                   >  2:259638/1‑251 (MQ=255)
                                                                                                                                                                                cacgaacatgaaggctttttCACGTTTCCCCAACGATGGTTGAAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCGCCTTCATAACTAATCAGCGGAATTGc                                                <  2:718461/234‑1 (MQ=255)
                                                                                                                                                                                               gCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCAtt                                 >  1:753442/1‑251 (MQ=255)
                                                                                                                                                                                                cTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCg                                                                                                                                                     >  2:561125/1‑134 (MQ=255)
                                                                                                                                                                                                cTTTCACGGTCCCCCAACGCTCGTGGTGATCCCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGTATCCAGGATGCg                                                                                                                                                     <  1:561125/134‑1 (MQ=255)
                                                                                                                                                                                                             ccccaacgcggtgggtgATACCGGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATa                     <  2:250988/238‑1 (MQ=255)
                                                                                                                                                                                                                 accgcgcGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGc                 <  1:577027/245‑1 (MQ=255)
                                                                                                                                                                                                                         gggtgATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCAtt         <  2:753442/250‑1 (MQ=255)
                                                                                                                                                                                                                                aCCAGCCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGcc  >  2:154819/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTAt                                        >  1:61548/1‑207 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
AACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGCC  >  NC_000913/3751310‑3751779

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