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,8300TC92.5% 148.7 / ‑0.5 53A114A (GCT→GCCnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (28/21);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGT  >  NC_000913/17586‑18075
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cAACCAGTGGATGGTAGAAATACTTTCTGGAGACGCTGGTTTAGCTCCGGGTTTGTGTACTTGAAACCAACAAAAACATGTTGGTAGGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTATGCGTTCTGGAAGTGAAACGTTAACTGATGCAAGGATCGCGAGCCAGCTTACTCCAGGCCGCATTTCTAGTTAGCGCCGCTATTGGCGGGATGATTGTTCCGGCATTACTCTATCTGGCctttcac                                                                                                                                                                                                                                                 <  2:477947/251‑4 (MQ=255)
 acccAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGCCGGGATGGAAGATAAACGTGAACTGCTTCAAGGTTCGCTAGCAAGCTTACGCCAGGCCGCATTTCCCGTTATCTCCCCTATTGGGGGGGTGATTGTTCCGGCATTACTCATTCTGGACTTTTact                                                                                                                                                                                                                                                 >  2:386026/3‑249 (MQ=255)
   ccAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTAt                                                                                                                                                                                                                                              <  1:284025/251‑1 (MQ=255)
                  aCGACTTTCTGGAGACGCCGGTTCAGCGCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTAcc                                                                                                                                                                                                                               <  2:100091/251‑1 (MQ=255)
                               gaCGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTggg                                                                                                                                                                                                                  <  1:221835/251‑1 (MQ=255)
                                    cGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGGGGGCGATc                                                                                                                                                                                                             >  2:132201/1‑251 (MQ=255)
                                                    tggttCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGAcc                                                                                                                                                                                             >  1:36916/1‑250 (MQ=255)
                                                    tggttCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGAc                                                                                                                                                                                              <  2:36916/250‑1 (MQ=255)
                                                         cACTCGACATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGGGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGcc                                                                                                                                                                                        >  1:169766/1‑250 (MQ=255)
                                                                        aaaaCATGCTGTTATGGATAAATGACGTGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTAATCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGGAGGGTGGGCGATCCCGGGGGCTACTGACATTTCTTTTTCACTTGgtgt                                                                                                                                                                          >  2:182706/1‑250 (MQ=255)
                                                                           aCATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTg                                                                                                                                                                      <  1:185605/251‑1 (MQ=255)
                                                                                cTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGTTTACCCGCGAAGGGTGGGCGATCCcg                                                                                                                                                                                                          >  2:140657/1‑210 (MQ=255)
                                                                                cTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCcg                                                                                                                                                                                                          <  1:140657/210‑1 (MQ=255)
                                                                                                   cgcTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGGGCTGTTGGGGAGTCGTGTTCCg                                                                                                                                              >  1:322430/1‑251 (MQ=255)
                                                                                                    gcTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGCAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGt                                                                                                                                             >  1:382567/1‑251 (MQ=255)
                                                                                                           gCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCt                                                                                                                                      >  1:290448/1‑251 (MQ=255)
                                                                                                           gCGGGATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGTTCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCGGGCGCTGTTGGGGAGTCGTGTTCCGTTAg                                                                                                                                          >  2:285238/1‑247 (MQ=255)
                                                                                                            cGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTCGCGCTg                                                                                                                                     >  2:185451/1‑251 (MQ=255)
                                                                                                                aTTTTTCCTTTTCGTCGGTCTGGAAGTTAAACGTGAACTGATGCTAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTCTTGGTGGGCTGATTGTTCCTGCATTACTCTATCTGTCCTTTAACTCTGCCGAGCTGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTCATCACATTGCATTTACCATTGGTGTGCTGTGGCGGTTTGGCAGTCGTGTTCCGttgtctctgaggg                                                                                                                                 >  2:361447/1‑240 (MQ=255)
                                                                                                                 tttttCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGa                                                                                                                                    >  1:282238/1‑247 (MQ=255)
                                                                                                                    ttCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCTCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCt                                                                                                                              <  1:383286/250‑1 (MQ=255)
                                                                                                                    ttCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCt                                                                                                                              <  1:132201/250‑1 (MQ=255)
                                                                                                                            aGTCGGTCTGGAAGTTAAACGTGACCTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATg                                                                                                                                                                                                                                             >  2:105019/1‑131 (MQ=255)
                                                                                                                            aGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATg                                                                                                                                                                                                                                             <  1:105019/131‑1 (MQ=255)
                                                                                                                               aggTCTGGCAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCAGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGAGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGctc                                                                                                                  >  2:399097/2‑251 (MQ=255)
                                                                                                                                  tCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGTCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTCCTGGCGCTGTTGGGAAGTCGTTTTCCGTTAGAGCTGACGATCTTTTTGCTGGCTCTgg                                                                                                               >  2:210267/1‑251 (MQ=255)
                                                                                                                                     ggAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGtt                                                                                                                                                              >  1:50448/1‑201 (MQ=255)
                                                                                                                                       aaGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGGGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtatt                                                                                                          >  2:332050/1‑251 (MQ=255)
                                                                                                                                        aGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGgtgtg                                                                                                                                                                         <  2:351114/187‑2 (MQ=255)
                                                                                                                                        aGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGgtgtg                                                                                                                                                                         >  1:351114/1‑186 (MQ=255)
                                                                                                                                                              ggATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGCCATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtattaagatcggaagggcgtcgtggtg                                                                                   >  2:105386/1‑229 (MQ=255)
                                                                                                                                                              ggATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtatta                                                                                                         <  1:105386/229‑1 (MQ=255)
                                                                                                                                                                             ttACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGtt                                                                    >  2:392386/1‑251 (MQ=255)
                                                                                                                                                                              tACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGtt                                                                    <  1:210267/250‑1 (MQ=255)
                                                                                                                                                                                     ggacGCATTTCCAGTTATCGCCGCTATTGGGGGGATGATTGTTCCGGCATTACGCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGAGCTTGGGGCCATCATTATCATCGCATTGTTCTACACTa                                                            <  2:344002/248‑1 (MQ=255)
                                                                                                                                                                                       ccGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATGGTTCCGGCATGACTCTACCTGGCCCTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAAt                                                          <  2:80854/251‑1 (MQ=255)
                                                                                                                                                                                           aTTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGCTCTTGGGGCCATCATTATCATCGCCTTGTTCTACACTAATGACt                                                      >  2:260224/1‑251 (MQ=255)
                                                                                                                                                                                                   ttATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATg                                              >  1:146200/1‑251 (MQ=255)
                                                                                                                                                                                                     aTCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGc                                            <  1:332050/251‑1 (MQ=255)
                                                                                                                                                                                                                tggtggGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGAGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCCTTGTTCTAAACTAATGACTTATCGATGGCCTCTCTTGGCg                                 >  1:480088/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggtggGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCg                                 <  2:412924/250‑1 (MQ=255)
                                                                                                                                                                                                                  gttggATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGAGCCGATTACCCGCGTAGGGTGGGCGAGCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGTAAGTCGTGTTCCGTTAGCGCTGAAGATCATTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCAGTGTTCTACACTAAGGACTTATCGATGGCCTCTCTTTGCGTc                               <  2:478237/248‑1 (MQ=255)
                                                                                                                                                                                                                         gatTGTTCCGGCAGTACTCTATCGGGCCTTTAACGATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCt                         <  2:290448/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCgg            <  1:132697/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGAGAGCGATACCGACGGCTAATGACATTGCCTTTGCACTTGGTGTGCTGGCGATGTTGGGAAGTCGTGTTACGTTAGAGCTGAAGATCATTTTGATGGATCTGGCTATTATCGATGATCTTGGGGCCATCATTATCATCACATTGTTCTACACTAATGCCTTATCGATGACCTCTCTTGGCGTCGCGGCTGTAGCAAATGCgg            >  1:35912/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGAGTGCTGGCGCTGTTGGGCAGTCGTGTTCCGTTAGCGCTGCAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTcgc     >  1:274423/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGCCGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCCATTGCGGTACTcgcg    >  2:6967/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTcgc     <  1:260224/250‑1 (MQ=255)
                                                                                                                                                                                                                                               cTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGAt                                                                                                                       <  1:209320/134‑1 (MQ=255)
                                                                                                                                                                                                                                               cTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGAt                                                                                                                       >  2:209320/1‑134 (MQ=255)
                                                                                                                                                                                                                                               cTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtatta                                                                                                         <  2:442611/148‑1 (MQ=255)
                                                                                                                                                                                                                                               cTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtatta                                                                                                         >  1:442611/1‑148 (MQ=255)
                                                                                                                                                                                                                                                tGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGt  >  2:276746/1‑250 (MQ=255)
                                                                                                                                                                                                                                                tGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGt  >  1:40707/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGT  >  NC_000913/17586‑18075

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