Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,485,547 (G)6→7 coding (429/1098 nt) yhdY → amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,485,5411.G92.9% 91.0 / 1.0 28coding (423/1098 nt)yhdYamino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (15/11);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

TGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCGGCGTTC  >  NC_002947/1485407‑1485672
                                                                                                                                       |                                                                                                                                   
tGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACg                                                                                                                                    >  1:81241/1‑136 (MQ=255)
tGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACg                                                                                                                                    >  2:323068/1‑136 (MQ=255)
tGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACg                                                                                                                                    >  1:362261/1‑136 (MQ=255)
          aCCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACgggggggtatc                                                                                                                          >  1:1281282/1‑131 (MQ=255)
          aCCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACgggggggtata                                                                                                                          >  2:1110802/1‑131 (MQ=255)
                 cGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCt                                                                                                                    >  1:443807/1‑136 (MQ=255)
                      ggCTTGCCGTGCTCGGCGCCGCCCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGaaa                                                                                                               >  2:110277/1‑136 (MQ=255)
                       gCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAg                                                                                                              >  1:908348/1‑136 (MQ=255)
                       gCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAg                                                                                                              <  1:260351/136‑1 (MQ=255)
                         ttGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCg                                                                                                            <  2:1329025/136‑1 (MQ=255)
                            ccGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTAc                                                                                                         >  1:522886/1‑136 (MQ=255)
                              gTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCg                                                                                                       <  1:323068/136‑1 (MQ=255)
                                              cgcTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGgggcg                                                                                       <  1:810962/136‑1 (MQ=255)
                                              cgcTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGgggcg                                                                                       <  1:910018/136‑1 (MQ=255)
                                                                        aaaGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGcc                                                             <  1:706831/136‑1 (MQ=255)
                                                                                            cTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGgcgc                                         >  2:48601/1‑136 (MQ=255)
                                                                                                       tGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTgg                              >  1:1353885/1‑136 (MQ=255)
                                                                                                        gTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTggg                             <  2:294506/136‑1 (MQ=255)
                                                                                                        gTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTggg                             <  2:522886/136‑1 (MQ=255)
                                                                                                            ccGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTg                                                                                            <  2:1025502/69‑1 (MQ=255)
                                                                                                            ccGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTg                                                                                            >  1:1025502/1‑69 (MQ=255)
                                                                                                            ccGGTGCTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATc                         >  2:454750/1‑136 (MQ=255)
                                                                                                                  cTGGCCTACACCCTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCctgctg                   <  1:888576/136‑1 (MQ=255)
                                                                                                                        tACACCCTGCTGCACGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTgg            >  1:305676/1‑136 (MQ=255)
                                                                                                                             ccTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCg        >  2:79446/1‑136 (MQ=255)
                                                                                                                             ccTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCg        >  1:1040490/1‑136 (MQ=255)
                                                                                                                             ccTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCg        >  1:618751/1‑136 (MQ=255)
                                                                                                                             ccTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCg        >  1:572485/1‑136 (MQ=255)
                                                                                                                             ccTGCTGCACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCg        <  1:28428/136‑1 (MQ=255)
                                                                                                                                  tgcACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCGGCGtt   >  2:1047368/1‑136 (MQ=255)
                                                                                                                                  tgcACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCGGCGtt   >  1:787428/1‑136 (MQ=255)
                                                                                                                                  tgcACGGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCGGCGtt   >  1:1231799/1‑136 (MQ=255)
                                                                                                                                  tgcACGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCGGCGTTc  >  2:853729/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                   
TGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGCCGCGCCGCTGTTCATCAAGCGTTTCCCGCGCAAAGCCTTCTACGGCCTGGGCTTCCTGGTGCTGTACCCGGTGCTGGCCTACACCCTGCTGCACGGGGGGTATCTGGGCCTGGAAAGCGTACCGACCAGCCAGTGGGGCGGGCTGATGCTGACCCTGGTGATCGCCACCGTGGGCATCGTCGGCGCCTTGCCGCTGGGCATCCTGCTGGCACTGGGCCGGCGTTC  >  NC_002947/1485407‑1485672

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