Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,893 T→C A135A (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,8930TC98.3% 199.5 / NA 61A135A (GCT→GCCnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (37/22);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGAC  >  NC_000913/17647‑18138
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
cGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTg                                                                                                                                                                                                                                                   >  1:522892/1‑251 (MQ=255)
               cATGCTGTTATGGCTAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCCAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTCACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGCCATTGCCCTTGCACTTGGTGTGCTgg                                                                                                                                                                                                                                    >  2:35596/1‑251 (MQ=255)
                    tGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAGGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTg                                                                                                                                                                                                                               >  1:110231/1‑251 (MQ=255)
                       tATGGATAAATTGCGCCTTGAAGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTATCATGCTTACGCCAGGCCGCATTTCCAGTTATCGCCTCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTCACTATGCCGATCCGATTATCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGTTGTGCTGGCGCTGtt                                                                                                                                                                                                                             <  2:342360/250‑1 (MQ=255)
                        aTGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTg                                                                                                                                                                                                                            <  1:160993/250‑1 (MQ=255)
                          ggATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTggg                                                                                                                                                                                                                          >  1:225054/1‑250 (MQ=255)
                          ggATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTggg                                                                                                                                                                                                                          >  1:481605/1‑250 (MQ=255)
                            aTAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGa                                                                                                                                                                                                                         >  1:193066/1‑249 (MQ=255)
                              aaaTGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGGGCTGGCGCTGTTGGGAAGTc                                                                                                                                                                                                                     >  1:189769/1‑251 (MQ=255)
                              aaaTGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGACAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGATATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGt                                                                                                                                                                                                                      <  2:141432/250‑1 (MQ=255)
                                                         ccTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCtttt                                                                                                                                                                                          >  1:617564/1‑251 (MQ=255)
                                                            gTTAGTCGGTCTGGAAGTTAAACGTTAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTg                                                                                                                                                                                        <  1:346816/250‑1 (MQ=255)
                                                                  cGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGTCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGct                                                                                                                                                                                  <  2:76147/250‑1 (MQ=255)
                                                                  cGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGTCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGct                                                                                                                                                                                  >  1:76147/1‑250 (MQ=255)
                                                                    gTCTGGAAGTTAAACGTGCACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATGGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACACGCGAAGGGTGGGCGATCCCGGCGGCTCCTGACATTGCATTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGATGAAGATCTTTTTGATGGCTCTg                                                                                                                                                                               >  1:464891/1‑251 (MQ=255)
                                                                        ggAAGATAAACGTTAACTGATGCAAGGATCGCTAGCCAGCGTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTGTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGc                                                                                                                                                                             <  1:63286/249‑1 (MQ=255)
                                                                          aaGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGGAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtctt                                                                                                                                                                         >  2:225729/1‑248 (MQ=255)
                                                                              tAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATcg                                                                                                                                                                     >  1:605389/1‑251 (MQ=255)
                                                                                   gTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCt                                                                                                                                                                                                     >  2:505384/1‑214 (MQ=255)
                                                                                   gTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCt                                                                                                                                                                                                     <  1:505384/214‑1 (MQ=255)
                                                                                     gAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATc                                                                                                                                                               <  2:470273/250‑1 (MQ=255)
                                                                                            tGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGGGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGCTCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGcc                                                                                                                                                       >  1:331329/1‑251 (MQ=255)
                                                                                               aaGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCcatc                                                                                                                                                    >  2:243645/1‑251 (MQ=255)
                                                                                                aGGATCGCTAGCCAGCTTACGCCAGGTCGCATTTCCAGTTATCGCCTCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGTCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGAGCGATCCCGGCGTCTACTGACATTGCCTTTGCACTTGGTGTGCTGTCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGACCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCcatca                                                                                                                                                   <  1:291562/251‑1 (MQ=255)
                                                                                                 ggATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGCCATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGACTATTATCGACGATCTTGGGGAcatcat                                                                                                                                                  >  1:332984/1‑251 (MQ=255)
                                                                                                         aGCCTGCTTACGCCAGGCCGCATTTCCAGTTTTCGCCGCTATTGGTGGGATGATTGTTCCGGCATTTCGCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTATTGACATTGCCTTTGCTCTTGGTGTGCTGGCGCTGTTGGGAAGTCTTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCTACGATCTTGGGGCCATCATTatcatc                                                                                                                                           <  2:443913/250‑1 (MQ=255)
                                                                                                          gCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTGTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATTGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGc                                                                                                                                         <  1:312838/251‑1 (MQ=255)
                                                                                                              gCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGACGGGTGGGCGATCCCGGCGGCTCCTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTg                                                                                                                                     >  1:83457/1‑251 (MQ=255)
                                                                                                                ttACGCCAGGCCGCATTTCCAGATATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTACCTTTGCACTTGGTGTGCTGGCTCTGTTGGGAAGTCGTGTTCCGTTCGCGCTGAATCTCTTTTTGATTGCTCTGGCTATTATCGACGATCTTGGGGCCCTCCTTCTCCTCGCATTGtt                                                                                                                                   >  2:462794/1‑251 (MQ=255)
                                                                                                                           cGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGCTTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAAtt                                                                                                                        >  1:452484/1‑250 (MQ=255)
                                                                                                                             cATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGAc                                                                                                                      <  1:555463/251‑1 (MQ=255)
                                                                                                                               tttCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTg                                                                                                                                                                                                                                     >  1:281938/1‑138 (MQ=255)
                                                                                                                               tttCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTg                                                                                                                                                                                                                                     <  2:281938/138‑1 (MQ=255)
                                                                                                                                ttCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACtt                                                                                                                    <  2:158800/250‑1 (MQ=255)
                                                                                                                                  ccAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATc                                                                                                                 >  2:257721/1‑251 (MQ=255)
                                                                                                                                   cAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATACCGGCGGATACTGACATTGCCTTTGCACTTGGTGTGCTGGAGCTGTTGGGAAGTCGTGTTCCGTTAGCGATGAAGATCTTTTTGATGGCTCTGGCTATTATAGACGATCTTGGGGCAATCATAATCATCGCATTGTTCTACCCTAATGACTTAt                                                                                                                  >  1:582769/1‑249 (MQ=255)
                                                                                                                                     gTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGAt                                                                                                              <  1:438784/251‑1 (MQ=255)
                                                                                                                                       tATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATgg                                                                                                            >  2:218369/1‑251 (MQ=255)
                                                                                                                                         tcgacgcTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCAGATCAGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGCTGGCTCTGGCTATTATCGACGATCTTGGGGCCATAATTATCATCGAATTGTTCTACACTAATGACTTATCGATGGcc                                                                                                          >  2:1493/5‑251 (MQ=255)
                                                                                                                                           gccgcTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCg                                                                                                                                          >  2:439524/1‑217 (MQ=255)
                                                                                                                                           gccgcTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCg                                                                                                                                          <  1:439524/217‑1 (MQ=255)
                                                                                                                                           gccgcTATTGGTGGGATGATTGATCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCctc                                                                                                        >  1:176829/1‑251 (MQ=255)
                                                                                                                                                  ttGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGc                                                                                                 >  2:269065/1‑251 (MQ=255)
                                                                                                                                                     gttggATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTc                                                                                              >  1:57963/4‑251 (MQ=255)
                                                                                                                                                                   ccGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAAt                                                                                <  1:497519/251‑1 (MQ=255)
                                                                                                                                                                      gCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGCTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGc                                                                             <  1:35596/251‑1 (MQ=255)
                                                                                                                                                                                  cTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTTTCGACGATCTTGGGGCCATCATTTTCATCGCATTGTTCAACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGt                                                                 >  1:351464/1‑251 (MQ=255)
                                                                                                                                                                                      cctttAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTg                                                             <  1:257721/249‑1 (MQ=255)
                                                                                                                                                                                            aCTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCtg                                                       <  1:218369/251‑1 (MQ=255)
                                                                                                                                                                                               aTGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCtgtgt                                                    >  1:527294/1‑251 (MQ=255)
                                                                                                                                                                                                        cgatTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTAcgc                                           >  1:192032/1‑251 (MQ=255)
                                                                                                                                                                                                               ccGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAAGTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACgg                                    <  2:282333/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcgAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACggg                                   >  1:300563/1‑251 (MQ=255)
                                                                                                                                                                                                                             cGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCtg                      <  2:189769/251‑2 (MQ=255)
                                                                                                                                                                                                                              gATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCtgg                     >  1:177964/1‑249 (MQ=255)
                                                                                                                                                                                                                                tCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGCTAGCGCTCCCGATCTTTTTGATGGCTCTGGCTATTATCGCCGACCTTGGGGCCCTCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCCCGGGCGTCTATATTCtggct                   >  1:560970/1‑247 (MQ=255)
                                                                                                                                                                                                                                      aggcTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGCAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCgt              >  2:60647/2‑250 (MQ=255)
                                                                                                                                                                                                                                       ggcTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCgt              >  1:525793/1‑249 (MQ=255)
                                                                                                                                                                                                                                        gcTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCgtagt           >  1:82276/1‑248 (MQ=255)
                                                                                                                                                                                                                                           aCTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAgtgt         <  1:258118/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTGCCTTTGCACTTGGTGTGCATGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGACTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTggcc  >  2:366499/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGAC  >  NC_000913/17647‑18138

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