Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,557,009 A→T 8.1% intergenic (+61/‑32) purU‑II → / → PP_5468 formyltetrahydrofolate hydrolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,557,0090AT 8.1% 140.4 / 2.1 49intergenic (+61/‑32)purU‑II/PP_5468formyltetrahydrofolate hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base T (2/2);  total (25/24)
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.68e-01

TGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGCACCGACTTC  >  NC_002947/1556875‑1557151
                                                                                                                                      |                                                                                                                                              
tGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCaa                                                                                                                                  >  2:93017/1‑149 (MQ=255)
    ccGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCTCTGCGcccctgtc                                                                                                                                   <  8:178922/144‑5 (MQ=255)
    ccGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCTCTGCGcccctgtc                                                                                                                                   >  7:178922/1‑140 (MQ=255)
           cTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGa                                                                                                                        >  5:345075/1‑148 (MQ=255)
           cTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCACTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGa                                                                                                                        >  6:295020/1‑148 (MQ=255)
                   tCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAAt                                                                                                               >  3:269349/1‑149 (MQ=255)
                          gaggaCCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAg                                                                                                        >  1:195853/1‑149 (MQ=255)
                             gaCCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAgcaa                                                                                                     >  6:357554/1‑149 (MQ=255)
                               ccGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAgc                                                                                                       <  8:145698/145‑1 (MQ=255)
                               ccGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAgc                                                                                                       >  7:145698/1‑145 (MQ=255)
                                      cTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGcaacaa                                                                                             <  7:178296/148‑1 (MQ=255)
                                         gTGCATGGTAACAAGCCTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAAcc                                                                                           <  8:303011/147‑1 (MQ=255)
                                             aTGGTAACAAGACTGTGGTGGTCGACTGACGCCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGGGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCGCGGCAAGCCGGGGGGGCGCGGCAGGGGGCAAGCACCGCGTCAGCCCCGTCg                                                                                     >  2:273963/1‑149 (MQ=255)
                                              tGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGc                                                                                    <  5:357554/149‑1 (MQ=255)
                                                    cAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTggg                                                                                                                           <  6:30353/104‑1 (MQ=255)
                                                    cAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTggg                                                                                                                           >  5:30353/1‑104 (MQ=255)
                                                     aaGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCt                                                                             >  4:249279/1‑149 (MQ=255)
                                                      aGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTc                                                                            <  6:345075/149‑1 (MQ=255)
                                                      aGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTc                                                                            <  5:295020/149‑1 (MQ=255)
                                                              gtgtTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATc                                                                    >  5:68944/1‑149 (MQ=255)
                                                                 ttCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCtt                                                                 <  1:93017/149‑1 (MQ=255)
                                                                       tGACACCGTCTGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCg                                                             >  8:216788/1‑147 (MQ=255)
                                                                       tGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCg                                                             >  7:49479/1‑147 (MQ=255)
                                                                           acCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTg                                                                >  2:155094/1‑140 (MQ=255)
                                                                           acCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTg                                                                <  1:155094/140‑1 (MQ=255)
                                                                           acCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCt                                                       <  3:249279/149‑1 (MQ=255)
                                                                                   gCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGcc                                               <  2:195853/149‑1 (MQ=255)
                                                                                         tgtgAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGc                                                <  5:63185/142‑1 (MQ=255)
                                                                                         tgtgAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGc                                                >  6:63185/1‑142 (MQ=255)
                                                                                          gtgAGATTCGGTGCTGCTGCGCAGCCCTGTCACGACTGGGTCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCGGTCTCCTCACCTCCACCGATCCTTGGCGATGGCTGCCTGGCCCGc                                            <  5:192730/145‑1 (MQ=255)
                                                                                          gtgAGATTCCGCGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCAACACCTCCACCGATCATTGGCGAAGGCTGCATGTCCCGCTag                                         >  6:192730/1‑147 (MQ=255)
                                                                                           tgAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTAc                                         <  6:68944/147‑1 (MQ=255)
                                                                                                    gggCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCt                                                       <  2:215308/124‑1 (MQ=255)
                                                                                                    gggCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCt                                                       >  1:215308/1‑124 (MQ=255)
                                                                                                        tgctgcGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCt                                   <  8:79787/140‑1 (MQ=255)
                                                                                                        tgctgcGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCt                                   >  7:79787/1‑140 (MQ=255)
                                                                                                        tgctgcGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTg                          >  8:326529/1‑149 (MQ=255)
                                                                                                           tgcGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATc                                                                    >  6:41034/1‑104 (MQ=255)
                                                                                                           tgcGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATc                                                                    <  5:41034/104‑1 (MQ=255)
                                                                                                               cAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGTCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGaca                   <  7:216788/149‑1 (MQ=255)
                                                                                                                     tGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTAc                                         <  2:148172/121‑1 (MQ=255)
                                                                                                                     tGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTAc                                         >  1:148172/1‑121 (MQ=255)
                                                                                                                       tCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGc           <  8:49479/149‑1 (MQ=255)
                                                                                                                           gACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGa                                       >  6:11801/1‑117 (MQ=255)
                                                                                                                           gACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGa                                       <  5:11801/117‑1 (MQ=255)
                                                                                                                             cAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGCACCGAc     <  7:326529/149‑1 (MQ=255)
                                                                                                                              cggggCGCTCTGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGCACCGAc     <  7:330232/147‑1 (MQ=255)
                                                                                                                              cggggCGCTCTGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGCACCGAc     >  8:330232/2‑148 (MQ=255)
                                                                                                                                gggCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGCACCGACTTc  >  5:343902/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                              
TGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGACACGGATGGCACCGACTTC  >  NC_002947/1556875‑1557151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑