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

GCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTT  >  NC_000913/3194115‑3194610
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
gCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTgggggggcg                                                                                                                                                                                                                                                        >  2:85519/1‑248 (MQ=255)
      gCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGGACTTTCTGATGAACAAACCAGCCTTAAATGCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGGGCCGATGAAGTCAATCGACGTTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAAGGTGGGGGGGCGGAAAcc                                                                                                                                                                                                                                                  <  2:430315/251‑1 (MQ=255)
           acGTAAAAAGGGTCTGGCTGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAGCACGGGTAATGTGGGTTGGCGGAGACCTGGCg                                                                                                                                                                                                                                             <  2:260776/251‑1 (MQ=255)
             gTAAAAAGGGTCCGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCgg                                                                                                                                                                                                                                            >  2:215460/1‑250 (MQ=255)
                    gggTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGAAGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGcc                                                                                                                                                                                                                                    >  2:111652/1‑250 (MQ=255)
                     ggTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCAt                                                                                                                                                                                                                                   >  1:323829/1‑251 (MQ=255)
                     ggTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCAt                                                                                                                                                                                                                                   >  1:307654/1‑251 (MQ=255)
                      gtcgggCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGGGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCAtt                                                                                                                                                                                                                                  <  2:111628/247‑1 (MQ=255)
                          ggCGGAAGCAGAAGCCCACCGTGCGCTGAACGATGCTATCAACGTACTTGCTGATGAACAAACCAGCCGTAAAGTCAAACAGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAAGCGACGGTATCAAGAGTATTCAGGTCTATGGTCTGAATCGTGGCGGCGCTGCGGGTGAGGCAAACACGGGGAATGGGGGGGGGCGGAAACCTGGCGGAACAAGCATTATc                                                                                                                                                                                                                               <  2:67148/250‑1 (MQ=255)
                                       gcaCAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCtt                                                                                                                                                                                                                 <  1:84967/251‑1 (MQ=255)
                                              gTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGCGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGcac                                                                                                                                                                                                          >  1:413824/1‑251 (MQ=255)
                                                  gcTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGcacac                                                                                                                                                                                                        >  1:308799/1‑249 (MQ=255)
                                                   cTGAACGATTCTATCAACCTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGATAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAGTGTGGGGGGGCGGAAACCTGGTGGAACAAGCATTATCAGCCGCTCTCTCTTACCGcacaca                                                                                                                                                                                                       <  1:25988/249‑1 (MQ=255)
                                                      aaCGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCAc                                                                                                                                                                                                  >  1:283838/1‑251 (MQ=255)
                                                            gCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTAAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTg                                                                                                                                                                                             >  1:372128/1‑250 (MQ=255)
                                                               aTCAACGTACTTTCGGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGGACACAGGCACCGCTGATTg                                                                                                                                                                                         <  2:47820/251‑1 (MQ=255)
                                                                tCAACGTACTTTCTGATGCAAAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCCGGCGCTGCCTGCGGTACTAGAGAAATCCGTTGAGCCGATGCAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACCCGGGTAATGTGGGGGGGCGGAACACTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCt                                                                                                                                                                                              >  2:43495/1‑245 (MQ=255)
                                                                tCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTAGTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAGCCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCt                                                                                                                                                                                              <  1:43495/245‑1 (MQ=255)
                                                                        cTTTCTGATGAACAAACCAGCCTTACATTCAAACTAGCCTTGTTGCCGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGAGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACCGGCAACGCTGATTGACTCCTTGc                                                                                                                                                                                >  2:61137/1‑251 (MQ=255)
                                                                        cTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAg                                                                                                                                                                                                      >  2:55387/1‑229 (MQ=255)
                                                                        cTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAg                                                                                                                                                                                                      <  1:55387/229‑1 (MQ=255)
                                                                                ggaACAAACCAGCCTTAAATTCAAACGGGCCGTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGaa                                                                                                                                                                        <  2:352836/250‑1 (MQ=255)
                                                                                  aaCAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGTTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCTCTCTCTCGTACCGCACACATGCACCGCTGATGGACTCCTTGCTCAAGGAAAt                                                                                                                                                                      <  2:200977/251‑1 (MQ=255)
                                                                                  aaCAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGaact                                                                                                                                                                      >  1:317509/1‑249 (MQ=255)
                                                                                        ccAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGt                                                                                                                                                                >  1:329588/1‑251 (MQ=255)
                                                                                          aGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGGTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCAGTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAACTTGGCGttt                                                                                                                                                              <  2:317509/251‑1 (MQ=255)
                                                                                                  aTTCAAACTGGCCTTGTTGCAGGCGATGCCTGCGGTAATAGAGAAATCCGATGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGAGGCGCTGCGGGTGATGCAAACCCGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAAGGGGCGTTACAggcggc                                                                                                                                                      <  2:46459/251‑1 (MQ=255)
                                                                                                     cAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATGCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGTTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCGCAATGAAATTGGCGTTTCAGGCGGCTCa                                                                                                                                                   <  2:350231/251‑1 (MQ=255)
                                                                                                      aaaCTGGCCTTGTAGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCAc                                                                                                                                                  <  1:308146/251‑1 (MQ=255)
                                                                                                       aaCTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACt                                                                                                                                                 >  1:417629/1‑251 (MQ=255)
                                                                                                              cTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCTGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCtt                                                                                                                                                                                  <  1:361563/211‑1 (MQ=255)
                                                                                                              cTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGGAATGTGGGGGGGCGGCAACCTGGCGGACCAAGAATTATCAGCCGCTCTCTCTTACCGCCCACGGGCCCCGCTGATTGACTCCtt                                                                                                                                                                                  >  2:361563/1‑211 (MQ=255)
                                                                                                               ttgttgCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCAt                                                                                                                                         >  2:91305/1‑251 (MQ=255)
                                                                                                                    gCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGAc                                                                                                                                     <  1:61137/250‑1 (MQ=255)
                                                                                                                          gcTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTCCCGCACACAGGCCCCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCCTTGACTTCAccc                                                                                                                              >  2:260503/1‑251 (MQ=255)
                                                                                                                           cggcctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAAACGACGGTATCAAGATTATGCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTTGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTGCAGGCGGCTCACTGGCGGCAGTGACTTCAccat                                                                                                                             <  2:40026/249‑3 (MQ=255)
                                                                                                                             gcctgcGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGAGGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCt                                                                                                                             <  1:215460/249‑1 (MQ=255)
                                                                                                                                          aagaAATCCGGTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTc                                                                                                              <  2:308799/250‑1 (MQ=255)
                                                                                                                                           agaAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTcc                                                                                                             >  2:136373/1‑251 (MQ=255)
                                                                                                                                             aaaTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCCTTATCAGCCGCTCTCTCTTACCGCACAAAGGCACCGCTGATTGAATCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATGGCCTTCACCCTTATCCTCAACAACTTCCg                                                                                                           >  2:44057/1‑251 (MQ=255)
                                                                                                                                              aaTCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGt                                                                                                          <  1:37394/251‑1 (MQ=255)
                                                                                                                                                      gAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGaaa                                                                                                  <  1:260503/251‑1 (MQ=255)
                                                                                                                                                          cGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGc                                                                                              <  1:85519/251‑1 (MQ=255)
                                                                                                                                                               aaGTCAATCGACGGTATCAAGATTATTCAGTTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGGAAACACGTGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTTCTTTCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGTCGGCATTGACTTCACCCTTAACCGCAACAACTCCCGTCGAAGAAAAAGCAGAGt                                                                                         <  1:257531/251‑1 (MQ=255)
                                                                                                                                                               aaGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGt                                                                                         <  1:4111/251‑1 (MQ=255)
                                                                                                                                                                 gTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTaa                                                                                       <  1:91305/251‑1 (MQ=255)
                                                                                                                                                                 gTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTa                                                                                        >  1:259813/1‑250 (MQ=255)
                                                                                                                                                                   cAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAc                                                                                     <  1:325370/251‑1 (MQ=255)
                                                                                                                                                                     aTCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAccc                                                                                   <  2:417629/251‑1 (MQ=255)
                                                                                                                                                                                   atttattCAGTTCGATGGTCTGAATCGGGGAGGAGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGaataa                                                                     <  2:337667/249‑1 (MQ=255)
                                                                                                                                                                                         tCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                <  1:114947/250‑1 (MQ=255)
                                                                                                                                                                                                gATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAACCCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                        <  1:166470/251‑1 (MQ=255)
                                                                                                                                                                                                 aTGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGAAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGCTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCCTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCCGGAATAATCCCTGGCCCaaa                                                        >  1:35695/1‑250 (MQ=255)
                                                                                                                                                                                                  tGGTCTGCCTCGTGGCGGCGGTGCGGGTGATGCCAACCCGGGTAATGTGGGGGGGCGGACACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGCCGAAGACAAAGCCGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaa                                                       >  2:254034/1‑250 (MQ=255)
                                                                                                                                                                                                   ggTCTGAATCGTGGCGGCTCGGCGGGTGATGCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCGCTTACCGCCCTCAGGCACCGCTGATTGCCTCCTTGCTCAAGGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAACCTCTTCCAGGAATAATCCCTGGCCCAATAAg                                                     <  2:233606/251‑1 (MQ=255)
                                                                                                                                                                                                    gTCTGAATCGTGGAGGCGCTGCGGGAGATGCAAACACGGGTAATGTGGGGGGGCGGACACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCCCTGGCGGCATTGACTTCTCCCTTAACCTCCACAACTCCCGTCGAAGAAAAAGCAGAGTAACCCCTCTTCCCGGAATAATCCCTGGCCACAAAAGc                                                    >  2:380353/1‑251 (MQ=255)
                                                                                                                                                                                                              tGGCGACGCTGCGGGTGATGCCAACCCGGGTAATGTGGGGGGGCGGCAGCCTGGCGGAACCCGCATTATCCGACGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAAAGAAATTGGCGTTTCCGGCGGCTCACTGGCGGCCTTGACTTCACCCTTCACCTCAACAACTCCCGTAGAAGCAAAAGCAGAGTAACCCCTCTTACAGGACTAATCCCTGACCCAACAAGCCGGGGAtttt                                          >  2:411127/1‑251 (MQ=255)
                                                                                                                                                                                                                ggggcgCTGCGGGGGATGCAAACACGGGTAATGTGGGTGGGCGTAAACCTGGCGGAACAAGCATTATCAGCCGCGCTCTCTTACCGCACACAGGCACCGCTGATTGAGTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCGCAACAACGCCCGTCGAAGAAAAATCGGAGGAAACCCTCTTCCATGAATAATGTCTGGCCCAAAAATCCGGGGAtttttt                                        <  2:65791/249‑1 (MQ=255)
                                                                                                                                                                                                                gcggcgCTGCGGGTGATGCAAACACTGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCTCACAGGCACCGCTGATTGACTCCTTGCTCAATGAAAGTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAGCCCGCTGCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                        <  1:296553/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggcgCTGCGGTTGATGCAAACACGGGTAATGTGAGGGGGCGGAAACCTGGCGGATCAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGGAATTGGTGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTGATCCTCAACAACTCCCTTCGTAGAGAAAGCAGAGTAAGCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGACGGGGATTTTTtat                                      <  1:399803/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCAAACACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCt                        <  1:15289/250‑1 (MQ=255)
                                                                                                                                                                                                                                   caaCACGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGCTTGCCTCCTTGCTCAATGCAATTGGCGTTTCCGGCGGCTCACTGGCGGCATTGACTTCACCCTTACCCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAACCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                      >  1:88540/2‑250 (MQ=255)
                                                                                                                                                                                                                                        gggaaaagggggggggCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCCCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGc               <  2:239323/241‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGGTAATGTGGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTg                >  1:136119/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       gggCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt  <  1:136373/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
GCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGTGGGGGGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTT  >  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: