Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,902,178 G→T 12.5% D169Y (GAC→TAC)  PP_2551 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,902,1780GT12.5% 105.2 / 3.1 32D169Y (GAC→TAC) PP_2551LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (19/9);  new base T (2/2);  total (21/11)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GTCGCTGGAGTACCGCACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATCGACCG  >  NC_002947/2902054‑2902308
                                                                                                                            |                                                                                                                                  
gtCGCTGGAGTACCGCACCTGGCCGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGCCCTGGCGACCACCGTGTTAGACGCGTCTCCGCCCACCCTGCACGGACGGCGGTTGGGTGAGGCCCGCGCGGTa                                                                                                                         >  2:146766/1‑136 (MQ=255)
   gCTGGAGTACCGCACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGc                                                                                                                      <  2:14423/136‑1 (MQ=255)
    cTGGAGTACCGCACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGcc                                                                                                                     >  1:301710/1‑136 (MQ=255)
    cTGGAGTACCGCACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGcc                                                                                                                     >  1:587138/1‑136 (MQ=255)
               cACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAg                                                                                                          <  1:204472/136‑1 (MQ=255)
               cACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAg                                                                                                          >  1:56399/1‑136 (MQ=255)
                  cTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGAt                                                                                                       >  1:178586/1‑136 (MQ=255)
                  cTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGAt                                                                                                       >  2:325368/1‑136 (MQ=255)
                     gCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGTACCGCGCGGTATGCCTGATGCGCCAGGATCAc                                                                                                    >  2:562082/1‑136 (MQ=255)
                     gCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGTACCGCGCGGTATGCCTGATGCGCCAGGATCAc                                                                                                    >  2:156363/1‑136 (MQ=255)
                                 gTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCcgcg                                                                                        >  2:577619/1‑136 (MQ=255)
                                 gTAAGCCTAGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTCCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCcgcg                                                                                        >  2:249885/1‑136 (MQ=255)
                                         tGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGATCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCcgc                                                                                <  2:423566/136‑1 (MQ=255)
                                             ggCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTc                                                                            >  2:634398/1‑136 (MQ=255)
                                                      agaAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGcc                                                                   >  1:93335/1‑136 (MQ=255)
                                                          aTCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGTACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACt                                                               <  1:156363/136‑1 (MQ=255)
                                                          aTCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGTACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACt                                                               <  1:562082/136‑1 (MQ=255)
                                                            cGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTAc                                                             <  1:236298/136‑1 (MQ=255)
                                                                   gCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGt                                                      >  1:80448/1‑136 (MQ=255)
                                                                                tCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTAc                                                             >  2:165935/1‑116 (MQ=255)
                                                                                tCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTAc                                                             <  1:165935/116‑1 (MQ=255)
                                                                                     gaGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGGTGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCAttt                                    >  2:587528/1‑136 (MQ=255)
                                                                                                    cTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACa                     >  1:61862/1‑136 (MQ=255)
                                                                                                         cggacggCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGAc                >  2:658882/1‑136 (MQ=255)
                                                                                                           gacggCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACTAGGACAg              >  1:124620/1‑136 (MQ=255)
                                                                                                           gacggCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAg              >  1:227461/1‑136 (MQ=255)
                                                                                                              ggCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCtt           >  2:175803/1‑136 (MQ=255)
                                                                                                                  gTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATc       >  1:374982/1‑136 (MQ=255)
                                                                                                                  gTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATc       <  2:533397/136‑1 (MQ=255)
                                                                                                                  gTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATc       <  2:554536/136‑1 (MQ=255)
                                                                                                                  gTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATc       <  1:181996/136‑1 (MQ=255)
                                                                                                                       gTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATCGACCg  <  1:634398/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                                  
GTCGCTGGAGTACCGCACCTGGCAGGCCGGCCAGTACGCCTTGCTGGCCAGCGGAGAAATCGACCTGGCGACCACCTTGTTCGACGAGTCGCCGCCCAACCTGCACGGACGGCTGTTGGGTGAGGACCGCGCGGTATGCCTGATGCGCCAGGATCACCCGCTGGCCGCGCAGGCCGCGCTCAGCCAGGCCGACTACCTGGCGTACAAGCATGTGCGCATTTCGGGCGGTGGCGACAAGGACAGCTTCATCGACCG  >  NC_002947/2902054‑2902308

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: ‑