Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 18,093 C→T A202V (GCA→GTA)  nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,0930CT93.9% 130.3 / 0.1 50A202V (GCA→GTA) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (19/27);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

GATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGG  >  NC_000913/17848‑18336
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
gatTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTACCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTAcgcc                                                                                                                                                                                                                                                >  1:487514/1‑251 (MQ=255)
    aCCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTGGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCa                                                                                                                                                                                                                                             <  1:63229/250‑1 (MQ=255)
           aaGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTc                                                                                                                                                                                                                                     <  1:264562/251‑1 (MQ=255)
                ggggCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTatat                                                                                                                                                                                                                                <  1:344266/250‑1 (MQ=255)
                  ggCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTc                                                                                                                                                                                                                              >  1:37044/1‑251 (MQ=255)
                         ccGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTgg                                                                                                                                                                                                                       >  1:204841/1‑251 (MQ=255)
                                  attGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAgtgt                                                                                                                                                                                                               <  1:85715/248‑1 (MQ=255)
                                    tgacatgtcctttggcatgggTGTGCTTGCGCGGTTGGGATGACGTGTTCCGTTATCGCTGAAGATTTTTTTGATTGCTCTGGCTATTATCGACGATTTTGGGGCCATCATTTTCATCGCATTGGTCTACGCTAATGACTTATCGATGGCCTCTGTTGGCGGCGCGTCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGGTGGCgtagtgctg                                                                                                                                                                                                             <  2:36399/231‑4 (MQ=255)
                                     gACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGtgttg                                                                                                                                                                                                             <  1:207181/249‑1 (MQ=255)
                                      aCATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTgg                                                                                                                                                                                                          >  1:322053/1‑251 (MQ=255)
                                      aCATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGGACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTgg                                                                                                                                                                                                          <  2:52879/251‑1 (MQ=255)
                                          tgcctttGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGCAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGCCTg                                                                                                                                                                                                      >  1:310817/5‑251 (MQ=255)
                                                  cACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTg                                                                                                                                                                                              <  1:406467/251‑1 (MQ=255)
                                                      tGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGaaa                                                                                                                                                                                           <  2:167830/250‑1 (MQ=255)
                                                        gtgtGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCGTTTTGATGGCGCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTGATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCg                                                                                                                                                                                        <  2:469814/251‑1 (MQ=255)
                                                          gtACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAAATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCTCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCttgt                                                                                                                                                                                                                          <  2:45379/215‑1 (MQ=255)
                                                          gtACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAAATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCTCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCttgt                                                                                                                                                                                                                          >  1:45379/1‑215 (MQ=255)
                                                                   cTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGc                                                                                                                                                                             >  1:237485/1‑251 (MQ=255)
                                                                               cGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGTCGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTAGCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTGGGGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCgg                                                                                                                                                                  <  2:310817/250‑1 (MQ=255)
                                                                                     ccGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCGTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTGGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTaa                                                                                                                                                            <  2:471954/250‑1 (MQ=255)
                                                                                         tAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGGGTTGCAATCGGGGGTTCTCGCAACCCTGGCGGGGGTGATTGTc                                                                                                                                                       >  2:198793/1‑251 (MQ=255)
                                                                                               tGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGGTGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGGATTGAATCTGTGTGGGGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGGAATTGGCgggttt                                                                                                                                                 >  2:54739/1‑247 (MQ=255)
                                                                                                gAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGcttc                                                                                                                                                 <  1:414017/250‑1 (MQ=255)
                                                                                                gAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGcttc                                                                                                                                                 <  1:812/250‑1 (MQ=255)
                                                                                                   gATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCtttt                                                                                                                                             >  2:373634/1‑250 (MQ=255)
                                                                                                        ttttGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTc                                                                                                                                          <  1:462238/249‑1 (MQ=255)
                                                                                                         tttGATGGCTCTGGCTATTATCGACGATCTTGGTGCCAGCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGACTCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTGGTGGTGTGGACTGCGGTGTTGGAATCGGGGGTTCACGCAACCCTGGCggggggtatttgagggttctttattactt                                                                                                                                       >  2:91719/1‑227 (MQ=255)
                                                                                                                               gacgaTCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGCAAGAGAAGCATGGGCGttt                                                                                                                 >  1:49336/1‑250 (MQ=255)
                                                                                                                                        ggggCCAACATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCagcga                                                                                                         <  2:151093/250‑1 (MQ=255)
                                                                                                                                          ggtcatcatTATCATTGCATTGTTCTACACTAATGACTGATCGATGGCCTTTCTGGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGCCTGCGGTGTTGAAATCGGGGGTTCACGCATCCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGTCGTTCTCCagcgaagc                                                                                                      <  1:417213/248‑1 (MQ=255)
                                                                                                                                                         cGCATTGTTCTGCACTAATGTCTTATCGAGGGCCTCTCTTGGCGTCGCGGCTGGAGCAAGTGCGGTACTCGTGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGGTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGTCTGGAGCAtgtg                                                                                        <  2:367871/250‑1 (MQ=255)
                                                                                                                                                            aTTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTAGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGGGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGGAATTGACAGATTCTTTATTCCTTTGGAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGAATGGAGCATGGGtta                                                                                     >  2:449121/1‑249 (MQ=255)
                                                                                                                                                                  cTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCGCCAGCGAAGCGACTGGGGCATGGGTTACACCCgg                                                                              >  1:104020/1‑250 (MQ=255)
                                                                                                                                                                     cacTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCggggg                                                                           >  1:77920/1‑251 (MQ=255)
                                                                                                                                                                       cTAATGACTTATCGATGGCCTCTATTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGCCTGCGGTGTTGACATCGGGGGATCACGCATCCCTGGCGGGGGGAATTGTCGGCTTATTTTTTCCTTTGAAAGAGAAGCATGGGCGATCTACAGCGAAACGACTGGCACATGTGttaccccatgggggg                                                                         >  2:18921/1‑238 (MQ=255)
                                                                                                                                                                          aTGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCCACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCg                                                                      >  1:282436/1‑251 (MQ=255)
                                                                                                                                                                            gattttgcGATTGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTGCTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACACGTGGGTGGCGTa                                                                    <  2:282183/244‑1 (MQ=255)
                                                                                                                                                                             ccttttcggtggcccttctTGGCGCCGCGGCTGTAGCAATTGCGGTACTCGCTGTAATGAATCTTTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACTCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGGGTTACTCCCGTGGGTGGCGTAt                                                                   <  2:439963/235‑1 (MQ=255)
                                                                                                                                                                                    gATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGAttt                                                            <  1:185467/251‑1 (MQ=255)
                                                                                                                                                                                       ggCCTCTCTTGGCGTCGCGGCTGTAGCCATTGCGGTACTCGCGGTATTGAATCTGGGTGGTGTACGCCGCACGGGCGTCTATATTCTGGGTGGCGTAGGGTTGTGGACTGCGGTGTTGAGATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGATTCTTTTTTTCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGGAACGACTGGGGCGTGTGTTGCACACGTGGGGGGCGGATCTGCTTTTGc                                                         >  2:423143/1‑251 (MQ=255)
                                                                                                                                                                                                    gTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCAtt                                            <  1:287907/251‑1 (MQ=255)
                                                                                                                                                                                                    gTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAAGCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTATTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTTGGTGGCGTATCTGATTTTGCCGCTGTTTGCAtt                                            <  1:15718/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaTTGCGGTACTCGCGTTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGTGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTGGCTAATGCTGGCg                              <  2:98053/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCa                         <  1:401710/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTACTCGCGGTATTGAAGCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCCGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCa                         <  1:115283/251‑1 (MQ=255)
                                                                                                                                                                                                                            cTCGCGGTATTGAATCTGTGTGGGGTACGCCGCACGGGCGTCTATATTCTGGGTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGCGAAGCATGGGCGGTCTCCAGCGAAGTGTCTGGTGCATGTGGTGCACCCCGGGGGGGCGTATCTGGTTTTGCCGCTATTTGCCTTTGCTAATCCTGGCGTTTCAATGc                     >  2:12370/1‑250 (MQ=255)
                                                                                                                                                                                                                               gcgGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTCTATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGAAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCAAGGGCGTTCTCCAGCGAAGCGGCTGGAGCATGTGTTACACCCGTGGGTGGCGTCTCTGATTTTGCCGCTGTTTGCGTTTGCTAATGCTGGCGTTTCACTGCAAgg                 >  2:353803/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGAGTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGTGGGGGGTAATTGTCGGCTTCTGTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAggtgt              <  1:458593/251‑4 (MQ=255)
                                                                                                                                                                                                                                      tGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACg          >  1:205500/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gtggtgTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGAGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATgg  <  1:300010/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
GATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGG  >  NC_000913/17848‑18336

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: ‑