Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,830 T→C A114A (GCT→GCC nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,8300TC97.8% 142.2 / ‑4.4 46A114A (GCT→GCCnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/24);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

CCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACT  >  NC_000913/17589‑18069
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
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)
            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)
                                                                                                                                                                                                                                gCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCTCTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGCTCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGTCGGCCTCTCTTTGCGTCCCGGCTGTAGCCATTg         >  2:66732/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tACTCTATCTGGCCTTTAACTATGCCGATCCGATTGCCCGCGAAGGGTGGGCGAGCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTg         <  2:198983/246‑1 (MQ=255)
                                                                                                                                                                                                                                    tACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTg         >  1:198983/1‑246 (MQ=255)
                                                                                                                                                                                                                                      ctctATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACt  >  1:324926/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
CCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACT  >  NC_000913/17589‑18069

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