Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 18,127 G→A V213V (GTG→GTA)  nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,1270GA94.4% 89.1 / ‑2.1 37V213V (GTG→GTA) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (20/14);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

ACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGC  >  NC_000913/17886‑18369
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
aCATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTgg                                                                                                                                                                                                                                            >  1:268768/1‑251 (MQ=255)
       ctttGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCgg                                                                                                                                                                                                                                     >  1:252939/2‑251 (MQ=255)
                 ggTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATc                                                                                                                                                                                                                           <  1:231511/251‑1 (MQ=255)
                                 tGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAAcc                                                                                                                                                                                                           >  1:371468/1‑250 (MQ=255)
                                                                     tGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCttt                                                                                                                                                                        <  1:367534/250‑1 (MQ=255)
                                                                              tGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAgaga                                                                                                                                                                >  1:141307/1‑249 (MQ=255)
                                                                                     tatCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATgggg                                                                                                                                                       >  1:160465/1‑250 (MQ=255)
                                                                                     tatCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGGCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGc                                                                                                                                                       >  1:7653/1‑251 (MQ=255)
                                                                                       tCGACGATCTTGGGGCCATCATTATCATAGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTCTTGAATCTGTGTGGTGTACGCCGCACGGGAGTCTATATTCTGGGTGGCGTAGTGTTGTGGACGGCGGTGTTGAAATAGGGGGTTCACGCAACCCTGGCGGGGGGAATTGTCGGCTTCTGTTGTCCTTTTCAAGAGTAGCATGGGCGt                                                                                                                                                     >  2:131864/1‑251 (MQ=255)
                                                                                             aTCTTGGGGCCATCATTATCATCGCACTGTTCTACACTAATGTCTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCACCCCTGGCGGGGGGAATTGTCGGCTTCTTTATTCCTTTGAAAGAAGAGCATGGGCGTTTTcc                                                                                                                                                >  1:332392/1‑250 (MQ=255)
                                                                                               cTTGGGGCCATCATTCTCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCACTATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGGAATCGGGGGTTCACGCCACACTGGCGGGGGATATCTTCGGCTTCTTTATTCCTTTTGAATAGAACCATGGGCGTTCTCCatg                                                                                                                                             >  2:147949/1‑249 (MQ=255)
                                                                                                 tGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCag                                                                                                                                              <  1:162672/248‑1 (MQ=255)
                                                                                                   gggCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTCTTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCagcga                                                                                                                                           >  1:32725/1‑249 (MQ=255)
                                                                                                     gCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCagcgaagcg                                                                                                                                       <  1:192929/251‑1 (MQ=255)
                                                                                                      ccATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGACGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAAATGTCGGCTTTTTTATTCCTTTTAAAGAAAAGACGGGGCGTTCTCCagcgaaacg                                                                                                                                       >  2:31829/1‑247 (MQ=255)
                                                                                                          catTAGCATCGCCTTGTTCTACACTAATGACTTATCGATGGCCTCGCTTGGCGTCGCGGCGTTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAATAGAAGCATGGGTGTTCTCCAGCGAAGCTACTgg                                                                                                                                  <  2:160465/251‑1 (MQ=255)
                                                                                                                catcGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCAtg                                                                                                                            >  1:299289/1‑251 (MQ=255)
                                                                                                                 cacGCATTGTTCGTCACTATTGACTTAGCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCTGTATTGAATCTGGTTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGTCTTCTTTATTCCTTTGAAAGAGAAGAATGGGCGTTCTCCAGCGAAGCGACTGGAGCAtgt                                                                                                                           <  2:33474/249‑1 (MQ=255)
                                                                                                                   cGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATgtgt                                                                                                                         >  1:274786/1‑251 (MQ=255)
                                                                                                                             tACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCgtg                                                                                                               >  1:129482/1‑251 (MQ=255)
                                                                                                                                      gACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTCTGAAAGAGAAGCATGGGCGTTCTCCAGTGACGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGt                                                                                                       >  1:183828/1‑250 (MQ=255)
                                                                                                                                        cTTATCGATGGCCTCTCTTGGCGTCGCGGCTGAAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACAGGCGACTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTACTTTGAAAAAGAAACCTGGGAGATCTCCAGCAAAACGACTGGCGCATGTGTTACACCCATGGGTGGCgattc                                                                                                    >  1:165331/1‑247 (MQ=255)
                                                                                                                                            ttggttgCCTCTCTTGGCGTCGCGTCTGTAGGAATTTCGGTACGTGCGGTATTGAACCTGTTTGGTGTACGCCGCACGGGCGTCTATCTTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGa                                                                                                 <  2:311070/244‑1 (MQ=255)
                                                                                                                                               aTGGCCTCTCTTGGCGTCGCGGCTGTAGCAATGGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACTCCCGTGGGTGGCGTATCTGAtttt                                                                                             <  1:31829/251‑1 (MQ=255)
                                                                                                                                                 ggCCTTTCTTGGCGTCTCGTCTGTATCAATTGTGGTAATAGCGGTATTGAATCTGTGTGTTGTACGCCGCACGGGTGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCGTCTTTATTCCTGTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTg                                                                                            <  2:371468/250‑1 (MQ=255)
                                                                                                                                                   ccTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCg                                                                                         <  1:203872/251‑1 (MQ=255)
                                                                                                                                                           ggCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTg                                                                                  <  1:258969/250‑1 (MQ=255)
                                                                                                                                                              gTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTCACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCAtt                                                                              >  1:148865/1‑251 (MQ=255)
                                                                                                                                                                    gCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGCCTGTAGCATGTGTTACACCCGTGGGGGGCGTATCTTGTTTTGCCGCTGTTTgactttgcga                                                                         >  2:137060/1‑241 (MQ=255)
                                                                                                                                                                         aGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTg                                                                   >  1:226657/1‑251 (MQ=255)
                                                                                                                                                                                 cGGTACTCGTGGTATTGAATCTGTTTGGTGTACGCCGCCCGGGCGGCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCa                                                          <  2:226657/251‑1 (MQ=255)
                                                                                                                                                                                     aCTCGCGGTATTGACTCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCCGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGGGGCGGATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCCCTGc                                                      >  1:86047/1‑251 (MQ=255)
                                                                                                                                                                                                           ggggtaCGCCGCGCGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGGTCACGCAACCCTGGCGTGGGTAATTGTCTGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCACTGCAAGGTGTCACGCTGGATGGCtt                                <  2:274786/246‑1 (MQ=255)
                                                                                                                                                                                                                  gccgcACGAGCGTCTATATTCTGGTTGGCGTCGTGTTGTGGACTGCGGTGTTGCAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCCTGGGCGTTCTCCAGCGAAGCGCCTGGGGCCTGTGTTAAACCCGTTGTGGGCGTCTCTGATTTTGCCGCTTTTTGCATTTGCTAATG‑TGGACGTTTCACTTAAAGGTGTCACGCTGGATGGCTTTACCTcc                         >  2:124633/1‑251 (MQ=255)
                                                                                                                                                                                                                         gggCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGc                  <  1:99714/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggttgGCGTAGTGTTGTGGACTGTGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTTTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGatcat     <  1:262844/248‑1 (MQ=255)
                                                                                                                                                                                                                                         gttgGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGc  >  1:65068/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGG‑CGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGC  >  NC_000913/17886‑18369

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