Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,912,591 T→C H354R (CAC→CGC)  relA ← GDP/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,912,5910TC97.1% 111.1 / ‑4.0 34H354R (CAC→CGC) relAGDP/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/15);  total (19/15)
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.38e-01

ACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATC  >  NC_000913/2912455‑2912727
                                                                                                                                        |                                                                                                                                        
aCTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGcgcg                                                                                                                                       >  2:137448/1‑140 (MQ=255)
aCTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGcgcg                                                                                                                                       >  2:91840/1‑140 (MQ=255)
aCTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGcgcg                                                                                                                                       <  2:829642/140‑1 (MQ=255)
aCTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGcgcg                                                                                                                                       >  1:829642/1‑140 (MQ=255)
aCTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGcgcg                                                                                                                                       >  1:719590/1‑140 (MQ=255)
   aCGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCgcag                                                                                                                                    <  1:91840/140‑1 (MQ=255)
       aCTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAac                                                                                                                                >  2:574535/1‑140 (MQ=255)
       aCTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAac                                                                                                                                >  2:158685/1‑140 (MQ=255)
          tcgtcgAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACAcc                                                                                                                             >  2:365016/1‑140 (MQ=255)
          tcgtcgAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACAcc                                                                                                                             >  1:707043/1‑140 (MQ=255)
          tcgtcgAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACAcc                                                                                                                             >  2:127483/1‑140 (MQ=255)
                     ttCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCAt                                                                                                                  <  2:707043/140‑1 (MQ=255)
                     ttCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCAt                                                                                                                  >  1:836555/1‑140 (MQ=255)
                            gAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATg                                                                                                           >  1:297608/1‑140 (MQ=255)
                                    cATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGtt                                                                                                   <  1:539582/140‑1 (MQ=255)
                                       ctctTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTgg                                                                                                >  2:195279/1‑140 (MQ=255)
                                        tctTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGt                                                                                               <  1:127483/140‑1 (MQ=255)
                                         ctTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTg                                                                                              <  2:297608/140‑1 (MQ=255)
                                         ctTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTg                                                                                              <  2:204843/140‑1 (MQ=255)
                                                  aCGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTgg                                                                                     >  2:223799/1‑140 (MQ=255)
                                                        tCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCa                                                                               <  1:137448/140‑1 (MQ=255)
                                                                     ccAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCAc                                                                  >  1:661615/1‑140 (MQ=255)
                                                                     ccAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCAc                                                                  >  2:552829/1‑140 (MQ=255)
                                                                      cAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCAcc                                                                 <  1:158685/140‑1 (MQ=255)
                                                                          cAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGc                                                             <  1:145852/140‑1 (MQ=255)
                                                                            aTCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGccc                                                           <  1:365016/140‑1 (MQ=255)
                                                                               cGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAg                                                        >  1:504305/1‑140 (MQ=255)
                                                                               cGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAg                                                        >  2:178826/1‑140 (MQ=255)
                                                                                  tCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAAc                                                     <  1:552829/140‑1 (MQ=255)
                                                                                   cTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAAcc                                                    <  2:661615/140‑1 (MQ=255)
                                                                                   cTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAAcc                                                    <  2:91059/140‑1 (MQ=255)
                                                                                      cATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACg                                                 >  2:718603/1‑140 (MQ=255)
                                                                                             cGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGaa                                          <  2:719590/140‑1 (MQ=255)
                                                                                              gAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAAt                                         >  2:739760/1‑140 (MQ=255)
                                                                                               aaCGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATa                                        <  1:739760/140‑1 (MQ=255)
                                                                                                                   gcgcCCTCTTTATATTTCCAGCGCGCAGCAACACCCAACTCTGCACCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGttt                    >  2:78863/1‑140 (MQ=255)
                                                                                                                                 tttCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGt      >  1:441008/1‑140 (MQ=255)
                                                                                                                                   tCCAGCGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTaa    >  2:742456/1‑140 (MQ=255)
                                                                                                                                     cAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATc  >  1:440505/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
ACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATC  >  NC_000913/2912455‑2912727

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