Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,469,455 A→G 8.3% T611A (ACC→GCC) ‡ mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,469,4550AG 8.3% 148.8 / 3.2 48T611A (ACC→GCC) ‡mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base G (2/2);  total (23/25)
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

GAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGT  >  NC_000913/4469309‑4469597
                                                                                                                                                  |                                                                                                                                              
gAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTAca                                                                                                                                              >  4:9946/1‑148 (MQ=255)
   ttACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACg                                                                                                                                           >  2:321216/1‑149 (MQ=255)
       cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGt                                                                                                                                       <  4:168174/149‑1 (MQ=255)
       cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGt                                                                                                                                       >  2:43126/1‑149 (MQ=255)
         tAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGttt                                                                                                                                     <  2:34048/149‑1 (MQ=255)
            aaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGcc                                                                                                                                  <  4:34138/149‑1 (MQ=255)
                    aCTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGAc                                                                                                                          <  3:433882/149‑1 (MQ=255)
                        gCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaa                                                                                                                                   <  1:135219/136‑5 (MQ=255)
                        gCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaa                                                                                                                                   >  2:135219/1‑132 (MQ=255)
                          gACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGAt                                                                                                                    <  1:132532/149‑1 (MQ=255)
                          gACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGAt                                                                                                                    <  1:106105/149‑1 (MQ=255)
                                     agttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACAACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAAc                                                                                                         <  3:9946/149‑1 (MQ=255)
                                                          tGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAg                                                                                                            <  4:283846/125‑1 (MQ=255)
                                                          tGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAg                                                                                                            >  3:283846/1‑125 (MQ=255)
                                                                ggggCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAg                                                                              <  2:252332/149‑1 (MQ=255)
                                                                 gggCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAgg                                                                             >  2:35846/1‑149 (MQ=255)
                                                                    cTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCa                                                                          <  3:405614/149‑1 (MQ=255)
                                                                           cgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGtt                                                                   <  2:421318/149‑1 (MQ=255)
                                                                                    tggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATg                                                          >  1:83494/1‑149 (MQ=255)
                                                                                       tCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                       >  3:462212/1‑149 (MQ=255)
                                                                                       tCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                       >  1:92899/1‑149 (MQ=255)
                                                                                         aTTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGGAGGGCATGTGGTTGGCTTTATGGGCGa                                                     <  1:224087/149‑1 (MQ=255)
                                                                                         aTTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGa                                                     >  1:111331/1‑149 (MQ=255)
                                                                                          ttGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGAt                                                    >  4:196651/1‑149 (MQ=255)
                                                                                          ttGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGAt                                                    <  1:43126/149‑1 (MQ=255)
                                                                                          ttGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGAt                                                    <  2:196693/149‑1 (MQ=255)
                                                                                           tGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATg                                                   <  2:181140/149‑1 (MQ=255)
                                                                                             gTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGtt                                                                   <  1:64299/131‑1 (MQ=255)
                                                                                             gTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGtt                                                                   >  2:64299/1‑131 (MQ=255)
                                                                                                  gATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTaa                                            <  2:108897/149‑1 (MQ=255)
                                                                                                     atTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatga                                         >  4:396832/1‑149 (MQ=255)
                                                                                                                       gacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTAc                                                                                          <  3:354740/82‑1 (MQ=255)
                                                                                                                       gacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTAc                                                                                          >  4:354740/1‑82 (MQ=255)
                                                                                                                       gacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACgcgc                       >  1:36361/1‑149 (MQ=255)
                                                                                                                        acgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCt                      <  2:86041/149‑1 (MQ=255)
                                                                                                                          gaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGc                    >  3:355800/1‑149 (MQ=255)
                                                                                                                          gaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGc                    <  3:396832/149‑1 (MQ=255)
                                                                                                                            aCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGAt                                                    <  3:337802/115‑1 (MQ=255)
                                                                                                                            aCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGAt                                                    >  4:337802/1‑115 (MQ=255)
                                                                                                                                 cAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATc             >  3:44418/1‑149 (MQ=255)
                                                                                                                                  aaaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCg            >  4:346378/1‑149 (MQ=255)
                                                                                                                                   aaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                        >  2:183065/1‑104 (MQ=255)
                                                                                                                                   aaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                        <  1:183065/104‑1 (MQ=255)
                                                                                                                                         aaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGGATGGGGTTGGCTTTGTGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc     >  2:294113/3‑149 (MQ=255)
                                                                                                                                          cacaGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTAc                           <  1:387339/126‑1 (MQ=255)
                                                                                                                                          cacaGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTAc                           >  2:387339/1‑126 (MQ=255)
                                                                                                                                           acaGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtg   >  4:262149/1‑149 (MQ=255)
                                                                                                                                            caGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtgt  <  1:90019/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGT  >  NC_000913/4469309‑4469597

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑