Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 18,253 G→A L255L (TTG→TTA nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,2530GA92.8% 175.5 / 3.6 71L255L (TTG→TTAnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (31/33);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

TCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTG  >  NC_000913/18009‑18492
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
tCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCg                                                                                                                                                                                                                                           >  1:547244/1‑251 (MQ=255)
        aaTGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCgtgggtggg                                                                                                                                                                                                                                   >  1:190933/1‑250 (MQ=255)
          tGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGAGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTCTTCCTTTGAAAAAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGt                                                                                                                                                                                                                                 >  1:293561/1‑251 (MQ=255)
                  cGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCCCGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGTGTTCACGCAACCCTGGCGGGGGGAATTGTCGGCTTCTTTATTTCTTTGAAAGAGAAGCCTGGGCGGTCTCCCGCGCAGCGACTGGAGCATGTGTTACACCCGAGGGTGGCGGATCTggtt                                                                                                                                                                                                                         >  2:176916/1‑248 (MQ=255)
                                       cgGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCt                                                                                                                                                                                                    >  1:629589/1‑251 (MQ=255)
                                                      cGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTc                                                                                                                                                                                      >  1:212514/1‑250 (MQ=255)
                                                                 gTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGGGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTGTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAggtgt                                                                                                                                                                          <  1:97844/251‑4 (MQ=255)
                                                                              gtggtgTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATgg                                                                                                                                                              <  1:80326/250‑1 (MQ=255)
                                                                                   gtaCGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGCTGGCTTGAc                                                                                                                                                        >  1:93678/3‑251 (MQ=255)
                                                                                            aCGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCt                                                                                                                                               >  1:429424/1‑251 (MQ=255)
                                                                                             cGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTTCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCCAGGTGTCAAGCTGGATGGCTTGACCTCCATTCTg                                                                                                                                              >  2:198592/1‑251 (MQ=255)
                                                                                               ggCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGTTTTTGCCGCTGTTTCCATTTGCTACTGCTGCCGTTTCACTGCAAGGTTTCCCGCTGGATGGCTTAACCTCCATTCtccc                                                                                                                                            >  1:125051/1‑248 (MQ=255)
                                                                                               ggCGTCTATATTCTGGTTGGCGTAGGGTTGTGGACTTCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGcc                                                                                                                                            <  1:421238/251‑1 (MQ=255)
                                                                                                 cGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGGAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGAATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCGGTTTGCATTTTCTAAGGCTGGCGTTTCACTGCAAGGGGTCACGGTGGATGGCTTGACCTCCATTCTGCCAt                                                                                                                                          >  2:269988/1‑251 (MQ=255)
                                                                                                      atatTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTgggg                                                                                                                                     >  1:652491/1‑251 (MQ=255)
                                                                                                       tattcgggttgGGGTGGTGGTGTGGTCTGCGGTGTTGAAATCGGGGGTGCACGCAACCCTGGCGGGGGTAATTGGCTGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCGGGATGGCTTGACCTCCATTCTGCCATTgggg                                                                                                                                     <  2:333330/243‑1 (MQ=255)
                                                                                                       gagtctgggtgGCGTAGTGTTGTGGACTGCGGTGTTTAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGATGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCGGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTgggg                                                                                                                                     <  2:49909/241‑1 (MQ=255)
                                                                                                         ttCTGGTTGGCGGTGTGTTGGGGACTGCGGTGTTGAAACCGGGGGTTCACGCAACCCTGGCGGGGTTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCGCCAGCGAAGCGACTGTAGCATGTGTTACTCCCGTTGGTGGCGTATCTGATTTTGCCGCTGTTTGCAGTTGCTAATGCTGTCGTTTCTCTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGat                                                                                                                                   <  1:665657/250‑1 (MQ=255)
                                                                                                             gatgggcttatggtggtgGGCTGCGTTGTTGAAATCGGGGGTTTACGCAACCCTGGCGGTGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCg                                                                                                                              <  2:201437/236‑1 (MQ=255)
                                                                                                                    gtagtgtTGTGGACTGCGGTGTTGAAATAGGGGGTTCACGCAAACCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCtt                                                                                                                       >  1:156159/4‑251 (MQ=255)
                                                                                                                    gtaggtgtgcggACTGTGGTGTTGAAATCGGGGGTTCCCGCAACCCTGGCGGGGGTTATTTTCGGGTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGGCTGGAGCATGTGTTACTCCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGGTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCtt                                                                                                                       <  2:135105/241‑1 (MQ=255)
                                                                                                                            gtgGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGAtt                                                                                                                <  1:28521/250‑1 (MQ=255)
                                                                                                                             tgGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTgg                                                                                                              >  1:619220/1‑251 (MQ=255)
                                                                                                                                aCTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCaa                                                                                                           <  1:432280/251‑1 (MQ=255)
                                                                                                                                   gCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAAcc                                                                                                        <  1:269988/251‑1 (MQ=255)
                                                                                                                                   gCGGTGTTGAAAGCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAAcc                                                                                                        <  1:625259/251‑1 (MQ=255)
                                                                                                                                          tGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTgggg                                                                                                                                     <  2:20084/215‑1 (MQ=255)
                                                                                                                                          tGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTgggg                                                                                                                                     >  1:20084/1‑215 (MQ=255)
                                                                                                                                            aaaTCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGAt                                                                                               >  1:439776/1‑251 (MQ=255)
                                                                                                                                                        cgcgcacccctGGCGGGGGTAATTGGCGGCTTCTTTAGTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACTCCCGGGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGTCAAACCGCTGGGGATTAGTCTTTTctg                                                                                   <  2:416352/242‑1 (MQ=255)
                                                                                                                                                        cACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTTAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTct                                                                                    <  2:212514/250‑1 (MQ=255)
                                                                                                                                                            cAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGa                                                                                                                                                                                                                           <  2:174956/111‑1 (MQ=255)
                                                                                                                                                            cAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGa                                                                                                                                                                                                                           >  1:174956/1‑111 (MQ=255)
                                                                                                                                                              accctGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGTAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGtt                                                                             <  1:605916/248‑1 (MQ=255)
                                                                                                                                                                   ggCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGgcg                                                                         <  1:147566/250‑1 (MQ=255)
                                                                                                                                                                   ggCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCAGGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGTTTGgc                                                                          <  2:619220/249‑1 (MQ=255)
                                                                                                                                                                    gCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGTTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCt                                                                       <  1:382493/251‑1 (MQ=255)
                                                                                                                                                                     cGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTg                                                                      <  1:66732/251‑1 (MQ=255)
                                                                                                                                                                       ggggTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCg                                                                    >  2:295524/1‑251 (MQ=255)
                                                                                                                                                                        gggTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCg                                                                                                                              <  1:543834/192‑1 (MQ=255)
                                                                                                                                                                        gggTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCg                                                                                                                              >  2:543834/1‑192 (MQ=255)
                                                                                                                                                                          gTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGttt                                                                 >  1:299148/1‑251 (MQ=255)
                                                                                                                                                                            taTTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGATCTGCCGCGAAGCGACTGGAGCATTTGATACACCCGTGGGCGGCGTATCCGCTTTTGCCGCTGTTTGATTTTGCTAATGCTGGCGTTTAACAGCAAGGTGTCACGCTGGTTGGCTTGACCTCCATTCAGCCATTGGGGATCATCGCTGGCTTGATGATTGGCAAACCGCTGGGGATTTGTCTGTTCTGCTGGTTGGCGCTGcagagac                                                               >  2:253559/2‑245 (MQ=255)
                                                                                                                                                                            aaTTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCCGATTTTGCCGCTGTTTTCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGTATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCg                                                                    <  1:253559/246‑1 (MQ=255)
                                                                                                                                                                                  cGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTg                                                          <  1:176916/250‑1 (MQ=255)
                                                                                                                                                                                   ggCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGgc                                                        >  1:646979/1‑251 (MQ=255)
                                                                                                                                                                                    gCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGgcg                                                       <  1:474213/251‑1 (MQ=255)
                                                                                                                                                                                       tcttTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCAt                                                    >  1:645363/1‑251 (MQ=255)
                                                                                                                                                                                        cttTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATc                                                   >  1:658958/1‑251 (MQ=255)
                                                                                                                                                                                                   tGAAAGAGATGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTTCACCCGTGGGGTGCGTATCTGATTTTGCCGCTGTTTGCAGTTGCTAAGGCTGGCGTCTCACTGCAAGGTGTCACGGTGGATGGCTTGACCTCCATTCTGCCATTGGGGAACATCGCTGGGTTGCTGAGTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGGTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAggg                                         <  2:49227/250‑1 (MQ=255)
                                                                                                                                                                                                   tGAAAGAGAAGCATGAGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTCTCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAACCCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGCggg                                         >  2:272124/1‑250 (MQ=255)
                                                                                                                                                                                                    gAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGaa                                       <  1:206444/251‑1 (MQ=255)
                                                                                                                                                                                                     aaaGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAAc                                      >  1:238671/1‑251 (MQ=255)
                                                                                                                                                                                                       agagAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACg                                     >  1:606845/1‑250 (MQ=255)
                                                                                                                                                                                                         agaAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACt                                  >  1:14061/1‑251 (MQ=255)
                                                                                                                                                                                                           aaGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACtt                                 >  1:472551/1‑250 (MQ=255)
                                                                                                                                                                                                              catgggggtgttcCAGCGCAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCa                             <  2:629589/240‑1 (MQ=255)
                                                                                                                                                                                                               aTGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATc                              >  1:156727/1‑249 (MQ=255)
                                                                                                                                                                                                               aTGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAg                            <  1:373750/251‑1 (MQ=255)
                                                                                                                                                                                                               aTGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGACGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCAATGCAAGGTGTCACGCTGGCTGGCTTGCCCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCACACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGCGGGAACGACTTATc                              >  1:610844/1‑249 (MQ=255)
                                                                                                                                                                                                                   gaggtctcCAGCGAAGCGACTGGGGCATGTGTTACGCCCGTGGGTGGCGTATCTGAATTTGCCTCTGTTTGCTTTGGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCAGTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTGCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCaaa                        <  2:160300/247‑1 (MQ=255)
                                                                                                                                                                                                                   gCGTTCGCCAGGGAAGCGACTGGAGCTTGTGTTGCTCCCGTGGGTGGCGTAGCTGATTTTGCCGCTGGTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCTCGCTGGATGGCTTGACCTCCATTCTGCCATTGGTGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCaaa                        <  2:624202/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTAGCTGATTTTGCCGCTGGTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAAt                       <  2:324926/251‑1 (MQ=255)
                                                                                                                                                                                                                           cAGCGAAGCGACGGGAGCATGTGGTACACCCGGGGGTGGCGTATCTGATTGTGCCGCTGGTTGCATTTGCTAATGCTGGCGTTGCACTGCAAGGTGTCACGCTGGATGGCTTGCCCTCTATTCTGCCATTGGGGATCATCGCTGGCGTGCTGATTGGAAAACCGCTGGGGATTAGTCTGTGCGGCTGGTTGGCGCTGCGTTTGAAACTGGAGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcg                <  2:79408/251‑3 (MQ=255)
                                                                                                                                                                                                                              cgaagcgaCTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcggtg             >  1:565513/1‑251 (MQ=255)
                                                                                                                                                                                                                              cgaagcgaCTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcggtg             >  1:481868/1‑251 (MQ=255)
                                                                                                                                                                                                                                 agcgaCTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTgggg          <  1:138594/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCtg    <  1:198592/250‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAGCATGTGTTACACCCGGGGGGGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCtg    <  2:565513/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gcaGCATGTGTTTCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCAGCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGGTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTAGCAGCAAATTATGGCGGTGGGGATCCt     <  2:429424/247‑1 (MQ=255)
                                                                                                                                                                                                                                         gAGCATGTGTTACACCAGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcggggtgggtggtgtg  >  2:61182/1‑235 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
TCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTG  >  NC_000913/18009‑18492

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