Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,174,739 G→A 10.8% A14T (GCA→ACA)  PP_2787 → transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,174,7390GA10.8% 96.6 / 3.0 37A14T (GCA→ACA) PP_2787transporter
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base A (2/2);  total (17/20)
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.74e-01

CTGTTCACGCAGACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCATTGATCTGCAGCGCGAGTTCACCG  >  NC_002947/3174592‑3174880
                                                                                                                                                   |                                                                                                                                             
cTGTTCACGCAGACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGc                                                                                                                                              <  6:172373/149‑1 (MQ=255)
  gTTCACGCAGACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCaa                                                                                                                                            >  6:237710/1‑149 (MQ=255)
        gCAGACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTgtgt                                                                                                                                      <  5:237710/149‑1 (MQ=255)
        gCAGACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGACAATTgtgt                                                                                                                                      >  4:192230/1‑149 (MQ=255)
           gACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCg                                                                                                                                   <  3:137563/149‑1 (MQ=255)
                       gcgcGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCATCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGa                                                                                                                       >  3:222343/1‑149 (MQ=255)
                                               gCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACAcc                                                                                                <  7:279403/148‑1 (MQ=255)
                                               gCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACAcc                                                                                                <  6:266684/148‑1 (MQ=255)
                                               gCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGACAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACAcc                                                                                                >  3:291353/1‑148 (MQ=255)
                                                   gAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCgg                                                                                           <  4:19538/149‑1 (MQ=255)
                                                    aaTACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTgtgt                                                                                                                                      <  4:31971/105‑1 (MQ=255)
                                                    aaTACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTgtgt                                                                                                                                      >  3:31971/1‑105 (MQ=255)
                                                                  gCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCa                                                                                                        <  2:154208/121‑1 (MQ=255)
                                                                  gCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCa                                                                                                        >  1:154208/1‑121 (MQ=255)
                                                                   cccGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCg                                                                              <  6:26020/146‑1 (MQ=255)
                                                                   cccGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCg                                                                              <  5:22108/146‑1 (MQ=255)
                                                                           gcgcCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCTGCGCGCTGATGTCTGacac                                                                         <  6:18865/143‑1 (MQ=255)
                                                                           gcgcCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGacac                                                                         >  5:18865/1‑143 (MQ=255)
                                                                                  tctcGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGaca                                                                          >  4:254976/1‑135 (MQ=255)
                                                                                  tctcGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGaca                                                                          <  3:254976/135‑1 (MQ=255)
                                                                                   ctcGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGc                                                           >  2:209922/1‑149 (MQ=255)
                                                                                       aTTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTc                                                         >  8:168612/1‑147 (MQ=255)
                                                                                            tCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGaca                                                  >  2:35358/1‑149 (MQ=255)
                                                                                                gAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTTGCGATCACCGCCTACTTCACGTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGacac                                                 <  4:222343/146‑1 (MQ=255)
                                                                                                gAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGaca                                                  <  3:240975/145‑1 (MQ=255)
                                                                                                gAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGACAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGaca                                                  <  3:192230/145‑1 (MQ=255)
                                                                                                      tcattcatGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGt                                                                 >  5:78524/1‑124 (MQ=255)
                                                                                                      tcattcatGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGt                                                                 <  6:78524/124‑1 (MQ=255)
                                                                                                              ggAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGa                                >  5:243458/1‑149 (MQ=255)
                                                                                                                         tGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGAAGATCATTGAt                     >  5:305160/1‑149 (MQ=255)
                                                                                                                          gAACGTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCATTGATc                    >  5:37052/1‑149 (MQ=255)
                                                                                                                           aaCTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCATTGAt                     <  2:62700/147‑1 (MQ=255)
                                                                                                                            aCTTCGTCGAGCGCTATGCGCGGACAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCATTGATCTg                  <  4:291353/149‑1 (MQ=255)
                                                                                                                                    gaccgcTATGCGCGGGCAATTGTGTTCGTCTTGGTGGCGATAACCGCCTACTTCACCTTCACCCTCGGCGCGATGATGTCGGACACAAATCCGTATCTGCTCAAGGACAATCACCCGTCGCTCAAGACGAACATTCAACTGCAgcga            >  4:220662/4‑146 (MQ=255)
                                                                                                                                          tATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCAt                         <  2:137537/128‑1 (MQ=255)
                                                                                                                                          tATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCAt                         >  1:137537/1‑128 (MQ=255)
                                                                                                                                            tGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCATTGATCTGCAGCGCGAGTTCACCg  <  6:243458/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTGTTCACGCAGACCGCGATCAAGCGCGGGCTGAAAGCCGGCGACTTGCCGGAATACGCCGTGTCTGCCCGCTAAGCGCCGCTCTCGATTCATCTGGAACCCTCATTCATGGAAAGATACCTGAACTTCGTCGAGCGCTATGCGCGGGCAATTGTGTTCGTGTTGGTGGCGATCACCGCCTACTTCACCTACACCCTCGGCGCGCTGATGTCGGACACCAATCCGTATCTGCTCAAGGACACTCACCCGGCGCGCAAGACGATCATTGATCTGCAGCGCGAGTTCACCG  >  NC_002947/3174592‑3174880

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