Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,424,386 G→A *229* (TAG→TAA ulaF → L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,3860GA98.2% 177.2 / NA 57*229* (TAG→TAAulaFL‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/28);  minor base C (1/0);  total (29/28)
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

ATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGC  >  NC_000913/4424255‑4424520
                                                                                                                                   |                                                                                                                                      
aTGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcgc                                                                                                                                 >  1:1224824/1‑139 (MQ=255)
     agaagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGc                                                                                                                            >  2:869224/1‑139 (MQ=255)
      gaagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCg                                                                                                                           <  1:1240351/139‑1 (MQ=255)
         gTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGaa                                                                                                                        >  1:563773/1‑139 (MQ=255)
             cGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAAttt                                                                                                                    >  2:538450/1‑139 (MQ=255)
                     gCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTc                                                                                                            <  1:869224/139‑1 (MQ=255)
                       gTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCgg                                                                                                          <  2:563773/139‑1 (MQ=255)
                               cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAGCACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAGATTTCCTTCCTCGGGAGATAAc                                                                                                  <  2:709230/139‑1 (MQ=255)
                               cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                  <  2:1210164/139‑1 (MQ=255)
                               cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                  <  2:1224824/139‑1 (MQ=255)
                                    ggCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTc                                                                                             >  1:783976/1‑139 (MQ=255)
                                         tAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTACACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAAtt                                                                                        >  1:458677/1‑139 (MQ=255)
                                                      aaTCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGtttt                                                                           >  2:570907/1‑139 (MQ=255)
                                                      aaTCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGAGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTCGTCTAATTCCGCAGCCGtttt                                                                           >  2:754880/1‑139 (MQ=255)
                                                                      tCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGccccc                                                           >  1:113110/1‑139 (MQ=255)
                                                                          gatgaATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCg                                                       >  1:1220163/1‑139 (MQ=255)
                                                                              aaTAAACACTTCATGTGTAAACACGGTCCTAACGCTTATTACGGACAGAAGTACAACACGCGCTGCGGAAATTTCCCTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgg                                                   >  2:915774/1‑139 (MQ=255)
                                                                              aaTAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                        <  2:720695/86‑1 (MQ=255)
                                                                              aaTAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                        >  1:720695/1‑86 (MQ=255)
                                                                                 aaaCACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAggggg                                                <  2:696504/139‑1 (MQ=255)
                                                                                 aaaCACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAggggg                                                <  2:655515/139‑1 (MQ=255)
                                                                                 aaaCACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAggggg                                                <  1:754880/139‑1 (MQ=255)
                                                                                   acacTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGa                                                                                                                         <  1:484641/64‑1 (MQ=255)
                                                                                   acacTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGa                                                                                                                         >  2:484641/1‑64 (MQ=255)
                                                                                   acacTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgca                                              <  1:570907/139‑1 (MQ=255)
                                                                                     acTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATa                                                                                                    >  2:1042816/1‑83 (MQ=255)
                                                                                     acTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATa                                                                                                    <  1:1042816/83‑1 (MQ=255)
                                                                                        tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGcc                                                              <  2:996598/118‑1 (MQ=255)
                                                                                        tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGcc                                                              >  1:996598/1‑118 (MQ=255)
                                                                                        tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                         >  2:362881/1‑139 (MQ=255)
                                                                                         cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgca                                              <  2:1022097/133‑1 (MQ=255)
                                                                                         cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgca                                              >  1:1022097/1‑133 (MQ=255)
                                                                                                      ggTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGt                           >  1:882267/1‑139 (MQ=255)
                                                                                                      ggTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGt                           >  2:444338/1‑139 (MQ=255)
                                                                                                        tCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGTGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgc                                               >  1:147508/1‑117 (MQ=255)
                                                                                                        tCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgc                                               <  2:147508/117‑1 (MQ=255)
                                                                                                          cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaa                                                                 <  1:961429/97‑1 (MQ=255)
                                                                                                          cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaa                                                                 >  2:961429/1‑97 (MQ=255)
                                                                                                          cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                       <  1:362881/139‑1 (MQ=255)
                                                                                                          cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                       <  1:444338/139‑1 (MQ=255)
                                                                                                          cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATCCCGCAGCAGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                       <  1:1127938/139‑1 (MQ=255)
                                                                                                             aCGCTTATTACGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                  >  1:5831/1‑60 (MQ=255)
                                                                                                             aCGCTTATTACGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                  <  2:5831/60‑1 (MQ=255)
                                                                                                                  tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGt                                                                                              >  2:174297/1‑60 (MQ=255)
                                                                                                                  tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGt                                                                                              <  1:174297/60‑1 (MQ=255)
                                                                                                                  tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGcccc                                                            <  2:1294431/94‑1 (MQ=255)
                                                                                                                  tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGcccc                                                            >  1:1294431/1‑94 (MQ=255)
                                                                                                                   attaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtata              >  1:881867/1‑139 (MQ=255)
                                                                                                                         ggCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                       >  2:57239/1‑124 (MQ=255)
                                                                                                                         ggCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                       <  1:57239/124‑1 (MQ=255)
                                                                                                                             gAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                  >  2:581748/1‑45 (MQ=255)
                                                                                                                             gAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                  <  1:581748/45‑1 (MQ=255)
                                                                                                                             aaaGTAAAACGCGCGCAGGGGAAATTTCCTTCACGGGAAGATAGCTGGTCTAATTCCGCCGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAGCAATGGCAATGTTTCGTTGGTTATACCTGGtatta    <  2:128249/138‑5 (MQ=255)
                                                                                                                              aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                        <  2:74293/38‑1 (MQ=255)
                                                                                                                              aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                        >  1:74293/1‑38 (MQ=255)
                                                                                                                              aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAg   <  2:1220163/139‑1 (MQ=255)
                                                                                                                               aGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGc  >  1:1121881/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
ATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGC  >  NC_000913/4424255‑4424520

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