Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 18,322 C→T G278G (GGC→GGT nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,3220CT94.4% 191.9 / 2.1 73G278G (GGC→GGTnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (26/41);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

GTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGG  >  NC_000913/18087‑18568
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
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                                                                                                                                                                                                                  >  1:20084/1‑215 (MQ=255)
                                                            tGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTgggg                                                                                                                                                                                                                  <  2:20084/215‑1 (MQ=255)
                                                              aaaTCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGAt                                                                                                                                                                            >  1:439776/1‑251 (MQ=255)
                                                                          cgcgcacccctGGCGGGGGTAATTGGCGGCTTCTTTAGTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACTCCCGGGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGTCAAACCGCTGGGGATTAGTCTTTTctg                                                                                                                                                                <  2:416352/242‑1 (MQ=255)
                                                                          cACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTTAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTct                                                                                                                                                                 <  2:212514/250‑1 (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                                                                                                                                                                                                           >  2:543834/1‑192 (MQ=255)
                                                                                          gggTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCg                                                                                                                                                                                                           <  1:543834/192‑1 (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:481868/1‑251 (MQ=255)
                                                                                                                                                cgaagcgaCTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcggtg                                                                                          >  1:565513/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)
                                                                                                                                                                   gttacACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTCGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGAAAATTATGGCGGCggggggacctctgggcggg                                                                        >  2:384153/5‑235 (MQ=255)
                                                                                                                                                                                gtggCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCCAGGTGTCACGCTGGATGGCTTGACCTTCAGTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcg                                                                                             <  2:50218/216‑3 (MQ=255)
                                                                                                                                                                                gtggCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcg                                                                                             >  1:50218/1‑214 (MQ=255)
                                                                                                                                                                                            tttttGCCGCTGTTGGCTATTGCTAATGCTGTCGTTTCAGTGCAAGGGGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCGGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGTTGGGGATCCTGGGCGGTATCGGTTTTACTATGTCTATCTTTATTg                                              <  2:93678/249‑1 (MQ=255)
                                                                                                                                                                                              tttGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGcc                                            <  2:547244/250‑1 (MQ=255)
                                                                                                                                                                                                  ccGCTGTTTGCATTTGCTAATGCTGGCGTTTCGCTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCt                                       <  2:156159/251‑1 (MQ=255)
                                                                                                                                                                                                           gCATGTGCTGATGCTGGAGTTTCACTGCAAGGTGTCACGCTGGATGGCTTTACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATGGGCAAACCGCTGGGGATTAGTCTGTTCTGCGGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTCTgg                              <  2:162933/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCa                 <  1:106932/249‑1 (MQ=255)
                                                                                                                                                                                                                           gCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGaa              <  1:488408/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGGTCCCGTGCGGTGTCGGTTTTTCTTATGTTTTCTTTTTTGTCATCCTGGGCTTTTGTAGGGTAGAttcaagacagatta       >  1:531547/1‑237 (MQ=255)
                                                                                                                                                                                                                                       aaGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTgg  <  1:196185/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
GTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGG  >  NC_000913/18087‑18568

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: