Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,041,632 T→A 100% V107E (GTG→GAG)  ygfZ → folate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,6320TA98.3% 73.3 / ‑6.6 33V107E (GTG→GAG) ygfZfolate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  major base A (14/15);  minor base C (1/1);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTT  >  NC_000913/3041493‑3041773
                                                                                                                                           |                                                                                                                                             
gccCAGTGCGACGCCAAAGGTAAGAGGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGAc                                                                                                                                            <  2:566014/143‑1 (MQ=255)
  ccATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCa                                                                                                                                          >  2:341270/1‑143 (MQ=255)
       gCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATcggg                                                                                                                                     >  1:324677/1‑141 (MQ=255)
                         tgtgGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGtgtg                                                                                                                   <  1:396681/143‑1 (MQ=255)
                              aGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                                              >  1:80180/1‑143 (MQ=255)
                                         tccggtgcgcgACGGCGAGGGCTTTGCATGGATTGAACGGCGCGGCGTGCGTGAACCGCAGCTGACTGAACTGAAGAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                                   <  2:271928/136‑1 (MQ=255)
                                          cTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGCCTGAACTGAACAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGCGCGTGTGCCGCTTGgtgctccccg                                                                                                  >  2:239392/1‑136 (MQ=255)
                                            gTTCCGCGACGGCGATGGCTTTGCATGGAGTGAACGGCGCAGCGTGCGGGAACCGCAGCTGACTGAACTGAAAAAATATGCGGGAATCTATAGAGCGAGCATCGCGACAGCACACGAGGACGTGCTGCCTGGCGTTGACggag                                                                                                >  2:413506/1‑141 (MQ=255)
                                               ccGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAGCCGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTCGGCGCTGCCGGtcctc                                                                                             >  2:8856/1‑139 (MQ=255)
                                                cgcgACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTACAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCCTGGCGTTGCCCGGTTTCa                                                                                            >  2:83877/1‑143 (MQ=255)
                                                cgcgACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGAGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                                                                            >  1:175022/1‑143 (MQ=255)
                                                         gATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGGGGGTGCTGCTTGGGGTGCGCGGTTTTGAGgcgcgcgc                                                                                   >  2:364366/1‑143 (MQ=255)
                                                          acggCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATCTTCTTTGTTTTTGTTTTTGTTCTTCTCTCCACACCACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgcc                                                                                  <  1:306102/141‑1 (MQ=255)
                                                          aTGGCTTTGCATGGATCGAACGGCGCCGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGCGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgcc                                                                                  <  1:491320/143‑1 (MQ=255)
                                                               tttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                                                             <  2:20270/143‑1 (MQ=255)
                                                                ttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCTGGCGcgcgccgcgccg                                                                            >  1:278413/1‑141 (MQ=255)
                                                                         ttGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGGGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGCCTAAtctc                                                                   >  1:328773/1‑143 (MQ=255)
                                                                          tGAACGGCGCAGCGTGCGTGAGCCGCAGCTGACGGAACTGCAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATctct                                                                  <  2:463665/143‑1 (MQ=255)
                                                                           gAACGCCCCAGCGTCCGTGACCCGCAGCTGACTGAACTGAAAACATATGCGGTAGCCCCTAAAGAGACCACCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCtt                                                                 <  2:467844/143‑1 (MQ=255)
                                                                                           gtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCg                                                 <  1:292221/143‑1 (MQ=255)
                                                                                                 cGCAGCTGTCTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATGGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaa                                           <  1:320063/143‑1 (MQ=255)
                                                                                                           tgaactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGCACTGCCTTCGAAAGAAAAACAGGTAg                                 >  2:17720/1‑143 (MQ=255)
                                                                                                                tgaaGAAGTATCCGGTATTCTCTAAAGAGACCATAGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCGGGTTTTCAGGCGCGCGCCGCGCTGGGAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCaaa                            <  2:297767/143‑1 (MQ=255)
                                                                                                                  aaaaaaTATGCGGTATGCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTGTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaaga                          <  1:69093/143‑1 (MQ=255)
                                                                                                                  aaaaaaGATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaaga                          <  1:635834/143‑1 (MQ=255)
                                                                                                                   aacaggatGCGGGATTCGCGAAAGAGGCCAGCGCAGCAGACGACGAGCGTGTGCGGTTTGGTGTTGCCGGGTGTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaa                         <  2:599149/137‑1 (MQ=255)
                                                                                                                    aaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGACGGGTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaag                        >  1:415175/1‑143 (MQ=255)
                                                                                                                                aTTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGCACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACtt            <  1:339807/143‑1 (MQ=255)
                                                                                                                                 ttCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACctt           >  2:503303/1‑143 (MQ=255)
                                                                                                                                     ctAAAGAGACCATCCCGCCAGTCTACTATCTTGTTCTGCTTGTTGTTTCCGTTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACATGTAGTCAAAGAAGGCGCGACCACTTTGCTa       <  1:612646/143‑1 (MQ=255)
                                                                                                                                       aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATg     <  2:478182/143‑1 (MQ=255)
                                                                                                                                       aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGACGCGACCACTTTGCTATg     >  1:324760/1‑143 (MQ=255)
                                                                                                                                          gTGACCATCGGGCCAGACGACGAGCGCGTTTTGCTTGGTGTTGCCGGTTTTCAGGCTCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGTTATGGtt  <  2:581870/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTT  >  NC_000913/3041493‑3041773

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