Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,809 G→T V107V (GTG→GTT nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,8090GT96.0% 140.1 / ‑5.2 50V107V (GTG→GTTnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (22/26);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCG  >  NC_000913/17566‑18047
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
gatTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAt                                                                                                                                                                                                                                         <  1:302690/251‑1 (MQ=255)
gatTATGGCCAATAGCGGCGCAACCAGAGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGCATTTTTCCTGTTAGTCGGTATGGAAGTTAAACGTGAACTGATGCAAGGGTCGCAAGCCCGCTTACGCCCGGCCGCCTTTCCCATTATCTCCGCTTTTGGGGGGGTGGTTGTTTCGGCAt                                                                                                                                                                                                                                         >  2:489911/1‑251 (MQ=255)
 atTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGGAGTTAAACGTGAACTGATGCAAGGTTCGCTAGCCAGCTTTCGCCAGGCCGCCTTTCCAGTTATCGCCGCTATTTGGGGGATGTGTGTTCCGGCAtt                                                                                                                                                                                                                                        >  2:8056/1‑251 (MQ=255)
  ttATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAtt                                                                                                                                                                                                                                        >  1:85633/1‑250 (MQ=255)
          aataGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGGCGGTCTGGAAGGTAAACGTGAACTGATGCAAGGATCGCTATCCAGCTTACGCCAGGCCGCATTTCCATTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCt                                                                                                                                                                                                                               <  2:293227/248‑1 (MQ=255)
                ggcgCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGcctt                                                                                                                                                                                                                         <  1:489911/251‑4 (MQ=255)
                       ccAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTAt                                                                                                                                                                                                                  <  1:507176/251‑1 (MQ=255)
                         aGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATTATTGTTCCGGCATTACTCTATCGGGCCTTTATCTATGc                                                                                                                                                                                                                <  2:85633/251‑1 (MQ=255)
                               tGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGccgatcc                                                                                                                                                                                                          >  1:385949/1‑251 (MQ=255)
                                ggtattatgacttgttggcgtcgctggTTCAGCTCCGGGTTGGTTCCATCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTGCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCATCTTACGCCAGGCCGCATTTCCAGTTATCGCCTCTATTGGGGGGAGGATTGTGCCGGCATTACTCTATCTGGGCTTTAACTATGccgatccg                                                                                                                                                                                                         <  2:86043/226‑1 (MQ=255)
                                 gttccACGACTTTCTGGAGACGCCTGTTCATCTCCGGGTTGGTTCTCTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATTCAAGGATCGCTAGCCAGCTTACTCCAGGCCGCATTTTCAGTTAGCGTCGCTATTGGTGGGATGATTGTTTCGTCATTACTCTATCTGGCCTTTAACTCTGccgatccg                                                                                                                                                                                                         <  2:572770/246‑1 (MQ=255)
                                  tATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAtt                                                                                                                                                                                                                                        <  2:88214/218‑1 (MQ=255)
                                  tATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAtt                                                                                                                                                                                                                                        >  1:88214/1‑218 (MQ=255)
                                   atcactcctttatgtagaCGCCCGTCCAGATACGGGTTGTTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTTATGGCGGTATTTTTCCTGTTATTCCGTCTGGAAGTTAACCGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCTTTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGTCCTTTACCTTTGccgatccgat                                                                                                                                                                                                       <  2:367840/235‑1 (MQ=255)
                                                       ccGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATATGGCCTTTAACTTTGCCGGTCCGATTACCCGCGAAGGGTGGGCGAt                                                                                                                                                                                  >  1:121784/1‑251 (MQ=255)
                                                                    ggggtggatCACTAGAAAGCAACAAAAACTTGCTTTTATGTATAAATGACGCGCTGATTGCGGTATGTCTCCTGTTAGTCGGTCTGGAAGGTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCGTACGCCAGGCCGCATTTCCAGTTATCGCCGCTTTTGGTGGGATGATTGTTCCGGCATTCCTCTATCTGGCCTGTACCTAGTCCGATCCGATTACTCGCTAAGGTTGGGCTAGCCCGGCGGTTACt                                                                                                                                                                     <  2:326576/243‑1 (MQ=255)
                                                                                  gAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTATATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGc                                                                                                                                                       >  1:73513/1‑251 (MQ=255)
                                                                                        aacaaAACCATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGCTAGCAGGAAGGGAAGTTAAACGTGAACTGGTGCAAGGGTCGCCAGCCAGCTTACGCCAGGCCGGATTTCCAGTTATGCCCGCCATTGGTGGGATGATTGTTCCGGCATTTCTCTTTCCGGCCTTTAACTATGCCGCTCCGATTACCCGCGCAGGGTGGGCGGTCCCggcggcgcagaccattgcatttgcccttgg                                                                                                                                                 >  2:611336/1‑227 (MQ=255)
                                                                                                        tATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGtt                                                                                                                                  <  1:205449/250‑1 (MQ=255)
                                                                                                             aTAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTATGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCCGTTATCGCCGCTATTGGTGGGATGATTGTTCAGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGCTTACCCGCGAAAGGTGGGGGATCCCGGCGGCTACTGCCATTGCCTTTGCACTTTGTGTGCTGGCGCTGTTgggga                                                                                                                             >  2:423580/1‑248 (MQ=255)
                                                                                                                aaTGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACCTTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCg                                                                                                                         >  1:626040/1‑251 (MQ=255)
                                                                                                                aaTGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCg                                                                                                                         >  1:258334/1‑251 (MQ=255)
                                                                                                                                cGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTg                                                                                                         >  1:135105/1‑251 (MQ=255)
                                                                                                                                         tCCTGTTAGTCGGTCGTGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCtt                                                                                                 <  1:25090/250‑1 (MQ=255)
                                                                                                                                             gTTAGGCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTg                                                                                             <  1:605847/250‑1 (MQ=255)
                                                                                                                                              ttAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACACGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCCCTTGGGGGGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAACGCTGTAGCTCTTTTTGAt                                                                                           >  2:225756/1‑251 (MQ=255)
                                                                                                                                                agtacgTCTGAAAGTGAAACTTGAATTCATGCAAGGAGTGCTAGCCTGCTTACTCCAGGCCGCATTTCCAGTTAGCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATgg                                                                                         <  2:472212/246‑1 (MQ=255)
                                                                                                                                                   cGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGctc                                                                                      <  1:8056/251‑1 (MQ=255)
                                                                                                                                                         ggAAGTTAAACGTGAACTGATGCAAGGTTCGCTAGCCAGCTTGCGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGACATTTCTCTATCTGGGCTTTAACTATGCCGCGCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGTCTTTGCACTTGGTGTCCTGGCGCTGGTGGGAAGGCGTGTTTCGTTAGCGCTGAAGATCTTTTTAATGGCCCTGGCta                                                                                <  2:238315/251‑1 (MQ=255)
                                                                                                                                                            aGTTAAACGTGAACTGTTGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTTTCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGTGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTTTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtatta                                                                             <  1:562080/251‑1 (MQ=255)
                                                                                                                                                             gTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtattat                                                                            >  1:9131/1‑251 (MQ=255)
                                                                                                                                                              ttAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATc                                                                           >  1:416352/1‑251 (MQ=255)
                                                                                                                                                               gtaccgtgaaatgtatgaaggttcccTAGCCAGCTTACTCCAGGCCACATTTCCATTGATCGTCGCTATTGGTGTGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGAGTCGATTACCCGCGTAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCCCGTGGTGTGCTGTCGCTTTTGTTAAGTCGTGTTCCGTTAGCTCTGAAGATCTTTTTGATGGCGCTGGCTAGTATcg                                                                          <  2:631843/226‑1 (MQ=255)
                                                                                                                                                                     tGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATc                                                                    >  1:427146/1‑251 (MQ=255)
                                                                                                                                                                              gCAAGGATCGCTAGCCAGCTTACGCCTAGCCGCATTTCCAGTTTCCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTGGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGGGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCca                                                           <  2:469232/251‑1 (MQ=255)
                                                                                                                                                                               cAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCcat                                                          <  1:647145/251‑1 (MQ=255)
                                                                                                                                                                                     tCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGCTCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGGGTTCCGTTAGCGCTGAAGGTCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCCTTat                                                    >  2:292247/1‑251 (MQ=255)
                                                                                                                                                                                        cTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCCTTatcat                                                 >  1:201437/1‑251 (MQ=255)
                                                                                                                                                                                            ccAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGc                                              >  1:215764/1‑250 (MQ=255)
                                                                                                                                                                                                cTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGt                                         >  1:333330/1‑251 (MQ=255)
                                                                                                                                                                                                   aCGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCt                                      <  1:423580/251‑1 (MQ=255)
                                                                                                                                                                                                          gCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTaa                               <  1:85317/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGTCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGAt                   <  2:385949/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTTATCGCCGCTAGTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACGTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGAt                   <  1:611336/251‑1 (MQ=255)
                                                                                                                                                                                                                        tATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATgg                 >  1:93461/1‑251 (MQ=255)
                                                                                                                                                                                                                         aTCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATgg                 <  1:603276/250‑1 (MQ=255)
                                                                                                                                                                                                                              cgcTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTTATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCctctc           <  1:290283/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTg        >  1:624202/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggttggATGATTGTTCCTGCATTTCTCTGTCCGGCCTTTAAGTATGCCGATCCGATTACCCGTTGAGGGTGGGCTATCCCGGCTGCTACTGAGATTGCCTTTGCACTTGGTCTGCTGGCTCTGTTGGGAAGTCGTGTTCCGTTATCGCTGAAGATCTTTTTGATGGCTCTGGCTATTACCTACGATCTTGGGGCCATCATGATCATCGAATTGTTCTACAATAATGACTTATCGATGGCCTCTCTTGGCg     <  2:176102/246‑1 (MQ=255)
                                                                                                                                                                                                                                       tggGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTcg  >  1:199864/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
GATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCG  >  NC_000913/17566‑18047

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: