Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,194,362 (G)6→7 coding (1499/1662 nt) yqiK → flotillin family inner membrane protein YqiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,194,3621.G95.2% 178.8 / 3.7 62coding (1499/1662 nt)yqiKflotillin family inner membrane protein YqiK
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (26/33);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

GCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GG‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTT  >  NC_000913/3194115‑3194610
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
gCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTg‑gg‑g‑gggcg                                                                                                                                                                                                                                                        >  2:85519/1‑248 (MQ=255)
      gCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGGACTTTCTGATGAACAAACCAGCCTTAAATGCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGGGCCGATGAAGTCAATCGACGTTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAAGGTG‑GG‑G‑GGGCGGAAAcc                                                                                                                                                                                                                                                  <  2:430315/251‑1 (MQ=255)
           acGTAAAAAGGGTCTGGCTGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAGCACGGGTAATGTG‑GG‑TTGG‑CGGAGACCTGGCg                                                                                                                                                                                                                                             <  2:260776/251‑1 (MQ=255)
             gTAAAAAGGGTCCGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCgg                                                                                                                                                                                                                                            >  2:215460/1‑250 (MQ=255)
                    gggTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGAAGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGcc                                                                                                                                                                                                                                    >  2:111652/1‑250 (MQ=255)
                     ggTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCAt                                                                                                                                                                                                                                   >  1:323829/1‑251 (MQ=255)
                     ggTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCAt                                                                                                                                                                                                                                   >  1:307654/1‑251 (MQ=255)
                      gtcgggCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGGGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCAtt                                                                                                                                                                                                                                  <  2:111628/247‑1 (MQ=255)
                          ggCGGAAGCAGAAGCCCACCGTGCGCTGAACGATGCTATCAACGTACTTGCTGATGAACAAACCAGCCGTAAAGTCAAACAGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAAGCGACGGTATCAAGAGTATTCAGGTCTATGGTCTGAATCGTGGCGGCGCTGCGGGTGAGGCAAACACGGGGAATGGG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATc                                                                                                                                                                                                                               <  2:67148/250‑1 (MQ=255)
                                       gcaCAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCtt                                                                                                                                                                                                                 <  1:84967/251‑1 (MQ=255)
                                              gTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGCGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGcac                                                                                                                                                                                                          >  1:413824/1‑251 (MQ=255)
                                                  gcTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGcacac                                                                                                                                                                                                        >  1:308799/1‑249 (MQ=255)
                                                   cTGAACGATTCTATCAACCTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGATAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAGTGTG‑GG‑G‑GGGCGGAAACCTGGTGGAACAAGCATTATCAGCCGCTCTCTCTTACCGcacaca                                                                                                                                                                                                       <  1:25988/249‑1 (MQ=255)
                                                      aaCGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCAc                                                                                                                                                                                                  >  1:283838/1‑251 (MQ=255)
                                                            gCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTAAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTg                                                                                                                                                                                             >  1:372128/1‑250 (MQ=255)
                                                               aTCAACGTACTTTCGGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGGACACAGGCACCGCTGATTg                                                                                                                                                                                         <  2:47820/251‑1 (MQ=255)
                                                                tCAACGTACTTTCTGATGCAAAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCCGGCGCTGCCTGCGGTACTAGAGAAATCCGTTGAGCCGATGCAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACCCGGGTAATGTG‑GG‑G‑GGGCGGAACACTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCt                                                                                                                                                                                              >  2:43495/1‑245 (MQ=255)
                                                                tCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTAGTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAGCCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCt                                                                                                                                                                                              <  1:43495/245‑1 (MQ=255)
                                                                        cTTTCTGATGAACAAACCAGCCTTACATTCAAACTAGCCTTGTTGCCGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGAGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACCGGCAACGCTGATTGACTCCTTGc                                                                                                                                                                                >  2:61137/1‑251 (MQ=255)
                                                                        cTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAg                                                                                                                                                                                                      >  2:55387/1‑229 (MQ=255)
                                                                        cTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAg                                                                                                                                                                                                      <  1:55387/229‑1 (MQ=255)
                                                                                ggaACAAACCAGCCTTAAATTCAAACGGGCCGTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGaa                                                                                                                                                                        <  2:352836/250‑1 (MQ=255)
                                                                                  aaCAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGTTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCTCTCTCTCGTACCGCACACATGCACCGCTGATGGACTCCTTGCTCAAGGAAAt                                                                                                                                                                      <  2:200977/251‑1 (MQ=255)
                                                                                  aaCAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGaact                                                                                                                                                                      >  1:317509/1‑249 (MQ=255)
                                                                                        ccAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGt                                                                                                                                                                >  1:329588/1‑251 (MQ=255)
                                                                                          aGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGGTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCAGTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAACTTGGCGttt                                                                                                                                                              <  2:317509/251‑1 (MQ=255)
                                                                                                  aTTCAAACTGGCCTTGTTGCAGGCGATGCCTGCGGTAATAGAGAAATCCGATGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGAGGCGCTGCGGGTGATGCAAACCCGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAAGGGGCGTTACAggcggc                                                                                                                                                      <  2:46459/251‑1 (MQ=255)
                                                                                                     cAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATGCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGTTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCGCAATGAAATTGGCGTTTCAGGCGGCTCa                                                                                                                                                   <  2:350231/251‑1 (MQ=255)
                                                                                                      aaaCTGGCCTTGTAGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCAc                                                                                                                                                  <  1:308146/251‑1 (MQ=255)
                                                                                                       aaCTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACt                                                                                                                                                 >  1:417629/1‑251 (MQ=255)
                                                                                                              cTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCTGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCtt                                                                                                                                                                                  <  1:361563/211‑1 (MQ=255)
                                                                                                              cTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGGAATGTG‑GG‑G‑GGGCGGCAACCTGGCGGACCAAGAATTATCAGCCGCTCTCTCTTACCGCCCACGGGCCCCGCTGATTGACTCCtt                                                                                                                                                                                  >  2:361563/1‑211 (MQ=255)
                                                                                                               ttgttgCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCAt                                                                                                                                         >  2:91305/1‑251 (MQ=255)
                                                                                                                    gCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGAc                                                                                                                                     <  1:61137/250‑1 (MQ=255)
                                                                                                                          gcTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTCCCGCACACAGGCCCCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCCTTGACTTCAccc                                                                                                                              >  2:260503/1‑251 (MQ=255)
                                                                                                                           cggcctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAAACGACGGTATCAAGATTATGCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTTGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTGCAGGCGGCTCACTGGCGGCAGTGACTTCAccat                                                                                                                             <  2:40026/249‑3 (MQ=255)
                                                                                                                             gcctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGAGGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCt                                                                                                                             <  1:215460/249‑1 (MQ=255)
                                                                                                                                          aagaAATCCGGTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTc                                                                                                              <  2:308799/250‑1 (MQ=255)
                                                                                                                                           agaAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTcc                                                                                                             >  2:136373/1‑251 (MQ=255)
                                                                                                                                             aaaTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCCTTATCAGCCGCTCTCTCTTACCGCACAAAGGCACCGCTGATTGAATCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATGGCCTTCACCCTTATCCTCAACAACTTCCg                                                                                                           >  2:44057/1‑251 (MQ=255)
                                                                                                                                              aaTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGt                                                                                                          <  1:37394/251‑1 (MQ=255)
                                                                                                                                                      gAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGaaa                                                                                                  <  1:260503/251‑1 (MQ=255)
                                                                                                                                                          cGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGc                                                                                              <  1:85519/251‑1 (MQ=255)
                                                                                                                                                               aaGTCAATCGACGGTATCAAGATTATTCAGTTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGGAAACACGTGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTTCTTTCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGTCGGCATTGACTTCACCCTTAACCGCAACAACTCCCGTCGAAGAAAAAGCAGAGt                                                                                         <  1:257531/251‑1 (MQ=255)
                                                                                                                                                               aaGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGt                                                                                         <  1:4111/251‑1 (MQ=255)
                                                                                                                                                                 gTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTaa                                                                                       <  1:91305/251‑1 (MQ=255)
                                                                                                                                                                 gTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTa                                                                                        >  1:259813/1‑250 (MQ=255)
                                                                                                                                                                   cAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAc                                                                                     <  1:325370/251‑1 (MQ=255)
                                                                                                                                                                     aTCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAccc                                                                                   <  2:417629/251‑1 (MQ=255)
                                                                                                                                                                                   atttattCAGTTCGATGGTCTGAATCGGGGAGGAGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGaataa                                                                     <  2:337667/249‑1 (MQ=255)
                                                                                                                                                                                         tCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                <  1:114947/250‑1 (MQ=255)
                                                                                                                                                                                                gATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAACCCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                        <  1:166470/251‑1 (MQ=255)
                                                                                                                                                                                                 aTGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGAAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGCTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCCTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCCGGAATAATCCCTGGCCCaaa                                                        >  1:35695/1‑250 (MQ=255)
                                                                                                                                                                                                  tGGTCTGCCTCGTGGCGGCGGTGCGGGTGATGCCAACCCGGGTAATGTG‑GG‑G‑GGGCGGACACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGCCGAAGACAAAGCCGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaa                                                       >  2:254034/1‑250 (MQ=255)
                                                                                                                                                                                                   ggTCTGAATCGTGGCGGCTCGGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCGCTTACCGCCCTCAGGCACCGCTGATTGCCTCCTTGCTCAAGGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAACCTCTTCCAGGAATAATCCCTGGCCCAATAAg                                                     <  2:233606/251‑1 (MQ=255)
                                                                                                                                                                                                    gTCTGAATCGTGGAGGCGCTGCGGGAGATGCAAACACGGGTAATGTG‑GG‑G‑GGGCGGACACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCCCTGGCGGCATTGACTTCTCCCTTAACCTCCACAACTCCCGTCGAAGAAAAAGCAGAGTAACCCCTCTTCCCGGAATAATCCCTGGCCACAAAAGc                                                    >  2:380353/1‑251 (MQ=255)
                                                                                                                                                                                                              tGGCGACGCTGCGGGTGATGCCAACCCGGGTAATGTG‑GG‑G‑GGGCGGCAGCCTGGCGGAACCCGCATTATCCGACGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAAAGAAATTGGCGTTTCCGGCGGCTCACTGGCGGCCTTGACTTCACCCTTCACCTCAACAACTCCCGTAGAAGCAAAAGCAGAGTAACCCCTCTTACAGGACTAATCCCTGACCCAACAAGCCGGGGAtttt                                          >  2:411127/1‑251 (MQ=255)
                                                                                                                                                                                                                ggggcgCTGCGGGGGATGCAAACACGGGTAATGTG‑GGTG‑GG‑CGTAAACCTGGCGGAACAAGCATTATCAGCCGCGCTCTCTTACCGCACACAGGCACCGCTGATTGAGTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCGCAACAACGCCCGTCGAAGAAAAATCGGAGGAAACCCTCTTCCATGAATAATGTCTGGCCCAAAAATCCGGGGAtttttt                                        <  2:65791/249‑1 (MQ=255)
                                                                                                                                                                                                                gcggcgCTGCGGGTGATGCAAACACTGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCTCACAGGCACCGCTGATTGACTCCTTGCTCAATGAAAGTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAGCCCGCTGCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                        <  1:296553/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggcgCTGCGGTTGATGCAAACACGGGTAATGTGAGG‑G‑GG‑CGGAAACCTGGCGGATCAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGGAATTGGTGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTGATCCTCAACAACTCCCTTCGTAGAGAAAGCAGAGTAAGCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGACGGGGATTTTTtat                                      <  1:399803/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCAAACACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCt                        <  1:15289/250‑1 (MQ=255)
                                                                                                                                                                                                                                   caaCACGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGCTTGCCTCCTTGCTCAATGCAATTGGCGTTTCCGGCGGCTCACTGGCGGCATTGACTTCACCCTTACCCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAACCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                      >  1:88540/2‑250 (MQ=255)
                                                                                                                                                                                                                                        gggaaaagggg‑gg‑g‑gg‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCCCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGc               <  2:239323/241‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGGTAATGTG‑GG‑G‑GGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTg                >  1:136119/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       g‑gg‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt  <  1:136373/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
GCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTG‑GG‑G‑GG‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTT  >  NC_000913/3194115‑3194610

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