Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,404,402 G→A S71L (TCG→TTG)  nuoB ← NADH:quinone oxidoreductase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,404,4020GA93.3% 87.6 / ‑1.3 30S71L (TCG→TTG) nuoBNADH:quinone oxidoreductase subunit B
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (17/11);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CATCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACGAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGATCTATGCGGGTGAGCGTATAATCCATCTTAATGCCTCGCGGTTAGCGTTGACGATTAGCGATACT  >  NC_000913/2404157‑2404652
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
caTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTc                                                                                                                                                                                                                                                       <  1:259640/251‑1 (MQ=255)
          gATGAATTTATCGAAGCAATGAACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATAAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCAACGTCATGCACCGCGGAAAACAAAGTAACCAtctaa                                                                                                                                                                                                                                              >  2:136632/1‑248 (MQ=255)
                 ttATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAAc                                                                                                                                                                                                                                       <  1:277079/250‑1 (MQ=255)
                      gACGCCCTGCACAACGGAATAAATATCGTACATAACACCAGAGTTGGCACAGGCAACCATTGAGATAACTCATTTTGGTTACAGAATCTGGTCATACAGACGCTGAATAACAGGTGCCATTTTGGTAAAGCAGGTTACTGCAACCACCATCAGGGCAGCCTGACGCGGCGAAGCACGAAATACTTCTGCGCCAAAACGAGCCACGTCATGCACCGCGGCAAACAATGTCACCATCTCAACGTAACAGCAgg                                                                                                                                                                                                                                 >  1:156628/1‑251 (MQ=255)
                           ccTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCGCCATCTCAACGTAACAGCAGGAAAGa                                                                                                                                                                                                                            <  1:262693/251‑1 (MQ=255)
                                     ggAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTCAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCACTACTTCTGCGCCAAAACGCGCCACGTACTGAACCGCGGTAAACAAAGTCACCa                                                                                                                                                                                                                                                   >  2:147270/1‑218 (MQ=255)
                                     ggAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCa                                                                                                                                                                                                                                                   <  1:147270/218‑1 (MQ=255)
                                                  tACATACCAACAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGATCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCACCGTAACAGCAGGAAAACCCGAAGTTATACGGCCAAATTgg                                                                                                                                                                                                     >  2:270955/1‑250 (MQ=255)
                                                                     cacaGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAACCGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGCAACAGCAGGAAAGACCGAAGATATACGGCCAAATTGAGTTTTTACGAcccacgtt                                                                                                                                                                                   >  2:259987/1‑245 (MQ=255)
                                                                     cacaGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGtt                                                                                                                                                                                   >  1:160620/1‑250 (MQ=255)
                                                                              ccATTGAGATATCCCATTTTGGTTCCAGCAGTTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCTCCATCAGGTCAGTCTGACGCGGCTAAGCACGCAATACTTCTTCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAATTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCa                                                                                                                                                                         <  1:131346/251‑1 (MQ=255)
                                                                                   gagaTAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAACGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAAc                                                                                                                                                                                >  2:113542/1‑239 (MQ=255)
                                                                                   gagaTAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAAc                                                                                                                                                                                <  1:113542/239‑1 (MQ=255)
                                                                                               tttGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATaa                                                                                                                                                         >  1:205830/1‑250 (MQ=255)
                                                                                                      ccAGCATCTGGTCATACAGACGCTGAATAACCGGAGCCATTTTGGAAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTATGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAAGACGttt                                                                                                                                                 >  1:236571/1‑251 (MQ=255)
                                                                                                         gCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGACTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAAATGAGTTTTTACGACCCCAGTTACCAATGTCCTTAAGCTTGCCCATAAACACGTTTTTg                                                                                                                                              >  2:80008/1‑251 (MQ=255)
                                                                                                                gTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActt                                                                                                                                       <  1:98942/251‑1 (MQ=255)
                                                                                                                       aGACGCTGAATAACCGGTGCCATTTTGGTACAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAAtt                                                                                                                                                                                                       >  2:137378/1‑180 (MQ=255)
                                                                                                                       aGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAAtt                                                                                                                                                                                                       <  1:137378/180‑1 (MQ=255)
                                                                                                                                     cGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACCCGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGtt                                                                                                                   >  2:43020/1‑250 (MQ=255)
                                                                                                                                                                   aaCCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGt                                                                                    >  1:79844/1‑251 (MQ=255)
                                                                                                                                                                          aTCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTg                                                                                                                                                                                                      >  2:3179/1‑130 (MQ=255)
                                                                                                                                                                          aTCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTg                                                                                                                                                                                                      <  1:3179/130‑1 (MQ=255)
                                                                                                                                                                                     cTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGAt                                                                  >  1:116592/1‑251 (MQ=255)
                                                                                                                                                                                       gACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGATCt                                                                >  1:79881/1‑251 (MQ=255)
                                                                                                                                                                                            ggcgAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACGCAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGAACCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGAAGGGGGTAACGGTCGTTCTCACCGTTGGGATCTATGCg                                                           >  2:43094/1‑251 (MQ=255)
                                                                                                                                                                                                  gCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACGAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCCGGGGGTAACGGTCGTTCTCACCGTTGGGATCTATGCGGGTg                                                       >  1:1450/1‑249 (MQ=255)
                                                                                                                                                                                                        aaTACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGATCTATGCGGGTGAGCGtata                                               <  1:43020/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccAAAACGCGCCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGATCTATGCGGGTGAGCGTATAATCCATCTTAAt                                   <  1:136632/251‑1 (MQ=255)
                                                                                                                                                                                                                           gcgcCACGTCATGCACCGCGGTAAACAAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGCTCTATGCGGGTGAGCGTATAATCCATCTTAATGCCTcgc                            >  2:226839/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     gAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGATCGATGCGGGTGAGCGTATAATCCATCTTAATGCCTCGCGGTTAGCGTTGACGATTAGCGATACt  <  1:43094/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
CATCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGACGCTGAATAACCGGTGCCATTTTGGTAAAGCAGGTTCCTGCAACCACCATCAGGTCAGCCTGACGCGGCGAAGCACGCAATACTTCTGCGCCAAAACGCGCCACGTCATGCACCGCGGTAAACGAAGTCACCATCTCAACGTAACAGCAGGAAAGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGGTCGTTCTCACCGTTGGGATCTATGCGGGTGAGCGTATAATCCATCTTAATGCCTCGCGGTTAGCGTTGACGATTAGCGATACT  >  NC_000913/2404157‑2404652

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