Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,174 T→A 15.6% intergenic (‑371/+26) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,1740TA15.6% 111.9 / 13.3 45intergenic (‑371/+26)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (4/3);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAA  >  NC_002947/5079032‑5079318
                                                                                                                                              |                                                                                                                                                
cAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGc                                                                                                                                           <  1:919279/150‑1 (MQ=255)
 aGTGCGACGATTGTTGTGAATGGGAAAGCACATTCGATGTTTGAGGCCGCTCAGCGGGTCAAAACCAAACCGCAATCATTACACAAGCGCATAGATCGACAGGCTCAAGAATTCATTAGCGAACAAACAAAAAGCCCGCTCTATTGAGCg                                                                                                                                          <  2:1212257/150‑1 (MQ=255)
   tGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTACAAACTACACCGCAATAAATACATAATAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCggg                                                                                                                                        >  1:173570/1‑150 (MQ=255)
              ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCtt                                                                                                                                     >  2:1369497/1‑142 (MQ=255)
              ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCtt                                                                                                                                     <  1:1369497/142‑1 (MQ=255)
               tgtgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                            <  1:790217/150‑1 (MQ=255)
                   gATGGGAAAGCACATTTGATTATTGAGGCCGCCCCGCGCCAACCATCCGCCCCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATCGGCGCACAACAAAACAGCCCGTTATATTGGTCTGGGCTTTTGATGGGGGAt                                                                                                                        <  2:1020731/150‑1 (MQ=255)
                        gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGc                                                                                                                     >  2:228282/1‑148 (MQ=255)
                        gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGc                                                                                                                     <  1:228282/148‑1 (MQ=255)
                        gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCgg                                                                                                                   >  2:1250487/1‑150 (MQ=255)
                                               ccGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCg                                                                                            >  1:966560/1‑150 (MQ=255)
                                                cGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAt                                                                                                                                                <  2:936949/97‑1 (MQ=255)
                                                cGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAt                                                                                                                                                >  1:936949/1‑97 (MQ=255)
                                                   tCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCg                                                                                        >  1:1198226/1‑150 (MQ=255)
                                                         gCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGt                                                                                  <  1:1250487/150‑1 (MQ=255)
                                                                   acacCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                        <  1:1257320/150‑1 (MQ=255)
                                                                     acCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                        <  2:697784/148‑1 (MQ=255)
                                                                     acCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                        >  1:697784/1‑148 (MQ=255)
                                                                              aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgg                                                                                                                             >  2:67846/1‑86 (MQ=255)
                                                                              aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgg                                                                                                                             <  1:67846/86‑1 (MQ=255)
                                                                              aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCt                                                                                                                                      >  1:632881/1‑77 (MQ=38)
                                                                              aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCt                                                                                                                                      <  2:632881/77‑1 (MQ=38)
                                                                                         gCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                  <  2:966560/150‑1 (MQ=255)
                                                                                           gTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                >  2:579532/1‑150 (MQ=255)
                                                                                                tCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACATGCAACTGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCg                                           >  1:797788/1‑150 (MQ=255)
                                                                                                   aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGAGCGCTTGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGGCCGTGa                                        >  1:86789/1‑150 (MQ=255)
                                                                                                   aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                        <  2:1198226/150‑1 (MQ=255)
                                                                                                   aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                        <  2:1060727/150‑1 (MQ=255)
                                                                                                     aGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATg                                      >  2:297959/1‑150 (MQ=255)
                                                                                                        cTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagc                                                                                                            <  1:850799/77‑1 (MQ=255)
                                                                                                        cTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagc                                                                                                            >  2:850799/1‑77 (MQ=255)
                                                                                                                       gcgcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                        <  2:1079264/98‑1 (MQ=255)
                                                                                                                       gcgcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                        >  1:1079264/1‑98 (MQ=255)
                                                                                                                         gcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGGCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTAc                  <  2:124935/150‑1 (MQ=255)
                                                                                                                         gcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTAc                  >  1:1192570/1‑150 (MQ=255)
                                                                                                                           acaaacaaAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCa                <  2:1192570/150‑1 (MQ=255)
                                                                                                                             aaacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                      >  2:456633/1‑110 (MQ=255)
                                                                                                                             aaacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                      <  1:456633/110‑1 (MQ=255)
                                                                                                                                caaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGc           <  1:579532/150‑1 (MQ=255)
                                                                                                                                  aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCTGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCa                        >  2:1249527/1‑135 (MQ=255)
                                                                                                                                  aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCTGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCa                        <  1:1249527/135‑1 (MQ=255)
                                                                                                                                  aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCtt         >  2:967431/1‑150 (MQ=255)
                                                                                                                                        ccGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                    <  1:950756/69‑1 (MQ=255)
                                                                                                                                        ccGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                    >  2:950756/1‑69 (MQ=255)
                                                                                                                                         cGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTGGCCGAGTACTCTATCCAGCTGAGCAATGAGTCCGGGACGGTAGGTTTGGGCCAGGTTCTCAATGGCTTACTGaca  >  2:589363/1‑148 (MQ=21)
                                                                                                                                              |                                                                                                                                                
CAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAA  >  NC_002947/5079032‑5079318

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