Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,922,713 A→C T376P (ACC→CCC)  sdaC → HAAAP family serine/threonine permease SdaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,922,7130AC94.1% 47.7 / ‑3.2 17T376P (ACC→CCC) sdaCHAAAP family serine/threonine permease SdaC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (9/7);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

AACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTA  >  NZ_CP009273/2922632‑2922795
                                                                                 |                                                                                  
aaCCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGc                                                                            <  1:375076/90‑1 (MQ=255)
 aCCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCg                                                                           <  2:452373/90‑1 (MQ=255)
         aCTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAatc                                                                   >  2:114336/1‑90 (MQ=255)
              gcTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATcgc                                                              >  1:285376/1‑90 (MQ=255)
                           gTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTc                                                 >  1:416441/1‑90 (MQ=255)
                                    tGGAGTGTTGCCACCCTGAACCCGAGCATCCTGGCTATAATTGAACCCCTTCGCGGGCCAATCATCGCGATGACACTGTTCCTCACGCCg                                        >  2:175731/1‑90 (MQ=255)
                                       aTTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATg                                     <  1:339772/90‑1 (MQ=255)
                                                    tGAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGa                        <  1:364601/90‑1 (MQ=255)
                                                     gAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGaa                       >  2:224960/1‑90 (MQ=255)
                                                     gAACCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGaa                       <  1:201485/90‑1 (MQ=255)
                                                      aaCCCGAGCATCCTGGGTATGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGaaa                      <  2:304681/90‑1 (MQ=255)
                                                                        aTGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAg    >  2:451486/1‑90 (MQ=255)
                                                                         tGATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGt   >  1:163125/1‑90 (MQ=255)
                                                                          gATTGCACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGAGGCAGATGAACGCAATTCAGAAAGTCCCGGCAAGGCGTAAGTa  >  2:232549/1‑90 (MQ=255)
                                                                          gATTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTa  >  2:44963/1‑90 (MQ=255)
                                                                           aTTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGaaa                      >  1:447504/1‑69 (MQ=255)
                                                                           aTTGAACCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGaaa                      <  2:447504/69‑1 (MQ=255)
                                                                                 |                                                                                  
AACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTA  >  NZ_CP009273/2922632‑2922795

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