Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 557,635 A→T L36Q (CTG→CAG)  folD ← bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913557,6350AT93.8% 39.8 / NA 16L36Q (CTG→CAG) folDbifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (10/5);  minor base G (0/1);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCG  >  NC_000913/557453‑557836
                                                                                                                                                                                      |                                                                                                                                                                                                         
ccAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGAGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAAc                                                                                                                                                       >  1:271607/1‑235 (MQ=255)
            cGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTTCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtt                                                                                                                           <  1:307002/251‑1 (MQ=255)
                        cGGCATTCAGAGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCATGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAAt                                                                                                                               >  1:260390/1‑235 (MQ=255)
                                        gATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTACGGGAGGTCATACGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATc                                                                                                >  1:118132/1‑250 (MQ=255)
                                            aGCTCCAGCAGCTCCGCTTCACTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAAGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCag                                                                                              >  1:273504/1‑248 (MQ=255)
                                                                 cTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGTTGCAATACGCGCCTGAACTTTGTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGttn                                                                                                                          <  1:60978/199‑2 (MQ=255)
                                                                                         gagCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAACTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTCGTGCCCGCAGTCCGGCTGCAATACGCGCCTGCACTTTTTGAGCACCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGCGGATTCCATCTGTTACGt                                                                       >  1:226507/1‑226 (MQ=255)
                                                                                             gggAGNNGANNNCGNNTTCTTCACAAGCCTTGCGTTTGCNTGCGACATAAATTTGCGATGCAGGGTTACNANCNACCANNACANCGGCCTGTCCNGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTNCNGCGCGCACCTGCTGCGCAATCGTTTTACCGTNAATAATCTTtgct                                                                                                      >  2:243366/1‑191 (MQ=255)
                                                                                             gggAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAGCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgct                                                                                                      <  1:243366/191‑1 (MQ=255)
                                                                                                  aCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCAGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGcc                                                    >  1:263311/1‑236 (MQ=255)
                                                                                                                                      gcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATcn  <  1:198811/250‑2 (MQ=255)
                                                                                                                                                    cGATGNCGNGNTANNANGCACNAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGcn                                                                                                                                                                                 >  2:146937/1‑60 (MQ=255)
                                                                                                                                                    cGATGCAGTGTTACTACCCACCAGCACAACGGCCGGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGATCAACTTCAGAGCGCACCTGCTGGGCa                                                                                                                                 <  1:146937/109‑1 (MQ=255)
                                                                                                                                                      aTGCAGTGTTACTACCCACCAGCACAACGGCCTGTCCTGGTCCCCGCAGTTCGGCTGCAATGCGCGCCTGACCTTTTTGAGCAACTTCagag                                                                                                                                                <  1:296239/92‑1 (MQ=255)
                                                                                                                                                            ggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGACTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGAAATCg                                                                                                                             >  1:182200/1‑105 (MQ=255)
                                                                                                                                                             gTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGAAATAAGCGCCTGCACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCCTCTGTTACGTAGATCGAAggg                                                            >  1:127240/1‑169 (MQ=255)
                                                                                                                                                                                      |                                                                                                                                                                                                         
CCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCG  >  NC_000913/557453‑557836

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