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,6350AT96.8% 98.1 / NA 31L36Q (CTG→CAG) folDbifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (12/18);  minor base G (0/1);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTG  >  NC_000913/557391‑557853
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          
gTCCGGATGAATACGTTCCAGCACTTAGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCCGAATGCCCTCGATGGCGTTGTCGGCCTTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTAATAAGAGCGGGAGACGAACCCGAATTCTTCACAAGCCTTGCGTTTGCGTGCGACATAAATTTGCGATGCAGAGTTACTACCCACCAGCCCAACGCCCTGTCCTg                                                                                                                                                                                                                      >  2:114101/1‑251 (MQ=255)
        gAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGc                                                                                                                                                                                                              >  1:92332/1‑251 (MQ=255)
            aCGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAggc                                                                                                                                                                                                          >  2:165249/1‑249 (MQ=255)
                gccagctttttgatgttacaattcgccgccgGTAACGGCATTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGTGAGGCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGGGACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCgg                                                                                                                                                                                                       <  1:140662/225‑1 (MQ=255)
                            aCGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACgcgc                                                                                                                                                                                          >  1:230939/1‑251 (MQ=255)
                                          cgccgGTAACGGCAGTTGAACCAGAATGCCATCGATGGGTTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCGTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCTGCTGCAATACTCGCCTGACCTTTTTGAg                                                                                                                                                                            <  2:251406/251‑1 (MQ=255)
                                                  aCGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTc                                                                                                                                                                     >  1:275757/1‑250 (MQ=255)
                                                        gTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAgcgc                                                                                                                                                              >  1:169116/1‑251 (MQ=255)
                                                           gAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCAc                                                                                                                                                            <  1:58354/250‑1 (MQ=255)
                                                               cAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACctgct                                                                                                                                                       >  2:107951/1‑251 (MQ=255)
                                                               cAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACctgct                                                                                                                                                       <  1:164882/251‑1 (MQ=255)
                                                                aGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGACCCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCCATACGCGCCTGAAATTTTTGAGCAACTTCAGAGCGCACctgctg                                                                                                                                                      >  1:134788/1‑251 (MQ=255)
                                                                aGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACctgctg                                                                                                                                                      <  1:303997/251‑1 (MQ=255)
                                                                aGAATGCCATCGATGGGTTTGTCGTCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACctgctg                                                                                                                                                      <  2:79178/251‑1 (MQ=255)
                                                                 gAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAATCCTTGCGTTTGCTTGCGACATAAATTTGCGATTCAGGTTTACTACGCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCTGTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCtgctgc                                                                                                                                                     <  1:141474/251‑1 (MQ=255)
                                                                        aTCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCg                                                                                                                                              <  1:257702/251‑1 (MQ=255)
                                                                         gcGATGGGGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGt                                                                                                                                             <  2:169116/250‑1 (MQ=255)
                                                                          cGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAAGACGCGCCTGAACTTTGTGAGCAACTTCAGAGCGCAcc                                                                                                                                                           <  1:248607/236‑1 (MQ=255)
                                                                          cGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCAGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCCACTTCAGCGCGCAcc                                                                                                                                                           >  2:248607/1‑236 (MQ=255)
                                                                                             cAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCTCCAGCACAACGGCCTGTCCTGGTGCCCGCAGGCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCttt                                                                                                                          <  1:140748/250‑1 (MQ=255)
                                                                                                      gATAAGCTCCAGCAGCACCGATTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAACTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATc                                                                                                                 >  1:221229/1‑250 (MQ=255)
                                                                                                         aaGCTTCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCagag                                                                                                             <  1:69502/251‑1 (MQ=255)
                                                                                                                  cagcTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGCTGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGCTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTc                                                                                                    >  1:119117/1‑251 (MQ=255)
                                                                                                                                tggaggtTTCCGGGAGGGCATCAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATACATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                                                        >  2:64701/5‑249 (MQ=255)
                                                                                                                                                  cATAAGTGCGGGAGACGAATCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCAGCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATTCGCCt                                                                    <  2:234422/251‑1 (MQ=255)
                                                                                                                                                     aaGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGGCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAAGCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAt                                                                  <  1:195597/250‑1 (MQ=255)
                                                                                                                                                                               tcACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTc                                         <  1:110363/249‑1 (MQ=255)
                                                                                                                                                                                        gtGCGGTTGCTTGCGACATAAATGTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTGTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCGGTCAGgtttcgttt                              <  2:119117/250‑1 (MQ=255)
                                                                                                                                                                                                         aTAAAGTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCTCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGATAGTATTCCATCTGTTACGTTGATCGAAGGGTATGCGCCTAGTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGaa             <  1:130706/251‑1 (MQ=255)
                                                                                                                                                                                                           aaaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCGGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACATCAGAGCGCACCTGCTGCGCAAGCGTTTTACCGGCAAGAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTGACGTAGATCGAAGGGGATGCGCCTATTGTGTCATAAGCGGGGCGCGCTGTCAGGTTTCTTGTCAGATGTATCGCTTGAAGc           <  2:134788/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCTCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTGCAGATTTATCGCGTGAAGCGACCTCTTg  <  2:248619/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          
GTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTG  >  NC_000913/557391‑557853

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