Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 633,010 T→C 9.5% E144G (GAA→GGA)  aldB‑I ← aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947633,0100TC 9.5% 112.8 / 3.1 42E144G (GAA→GGA) aldB‑Ialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base C (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TTCGGCAGGCTTGAGCACCACGCAGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCGGCTTGCCGTTGTCCCAGGTTTCGGTAATGGCCAG  >  NC_002947/632868‑633149
                                                                                                                                              |                                                                                                                                           
ttCGGCAGGCTTGAGCACCACGCAGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGa                                                                                                                                       >  7:92674/1‑149 (MQ=255)
        gCTTGAGCACCACGCAGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGa                                                                                                                                       <  4:174677/141‑1 (MQ=255)
        gCTTGAGCACCACGCAGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGa                                                                                                                                       >  3:174677/1‑141 (MQ=255)
                  cacGCAGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCtt                                                                                                                     >  5:42383/1‑149 (MQ=255)
                       aGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTgg                                                                                                                <  7:93973/149‑1 (MQ=255)
                            ccAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGgcgc                                                                                                            <  8:92674/148‑1 (MQ=255)
                             cAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGcgcg                                                                                                           >  6:241474/1‑148 (MQ=255)
                                  gCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATg                                                                                                       >  4:220090/1‑147 (MQ=255)
                                   ccAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGcag                                                                                                    >  7:46306/1‑149 (MQ=255)
                                         ccGGGGCAAGCTTCCAGGGAGCCATCAGGTTCGGGAAGTTCGAGGGGATGATCTGCGCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCg                                                                                              <  7:128939/149‑1 (MQ=255)
                                               cAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGACACCGAAATAg                                                                                        >  3:21885/1‑149 (MQ=255)
                                                      ccAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGt                                                                                 >  7:21287/1‑149 (MQ=255)
                                                       cAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTg                                                                                <  2:57880/149‑1 (MQ=255)
                                                           cAGCCATCAGTATCTTGAAGTTCCAGGGGATGATCTGCCCGACCACGCGCAGTGGCTCGTGGATGTGGTAGGCCACGTTGCCTTCATTGATTTCGGCAGCGCGGCCTTCCTGTGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCg                                                                            <  6:149663/149‑1 (MQ=255)
                                                              ccATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCgccgcc                                                                                                                       <  3:230816/103‑1 (MQ=255)
                                                              ccATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCgccgcc                                                                                                                       >  4:230816/1‑103 (MQ=255)
                                                                   aGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCTCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAg                                                                    <  3:220090/149‑1 (MQ=255)
                                                                       tCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTCCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGaa                                                                                   >  6:53209/1‑130 (MQ=255)
                                                                       tCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTCCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGaa                                                                                   >  6:49304/1‑130 (MQ=255)
                                                                       tCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTCCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGaa                                                                                   <  5:53209/130‑1 (MQ=255)
                                                                       tCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTCCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGaa                                                                                   <  5:49304/130‑1 (MQ=255)
                                                                          ggAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATg                                                             >  8:215477/1‑149 (MQ=255)
                                                                               ttCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTg                                                                                >  8:46228/1‑125 (MQ=255)
                                                                               ttCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTg                                                                                <  7:46228/125‑1 (MQ=255)
                                                                                    ggggATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGa                                                   <  6:42383/149‑1 (MQ=255)
                                                                                     gggATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAg                                                  <  8:200980/149‑1 (MQ=255)
                                                                                        atgatCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGt                                               <  5:241474/149‑1 (MQ=255)
                                                                                          gatCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGa                                                   <  1:47545/143‑1 (MQ=255)
                                                                                          gatCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGa                                                   >  2:47545/1‑143 (MQ=255)
                                                                                            tCTGCCCTACCACTCCCAGTTGCTCGTGGATGTTGTATGCCAAGGTGCCTTCATTGATTTCTGCAGCGCCTCCTTCCTTGGCTCGGATGCAGCCATCGAAATAGCTGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTc                                            <  5:271122/148‑1 (MQ=255)
                                                                                                cccGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGc                                                                                       >  5:103198/1‑101 (MQ=255)
                                                                                                cccGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGc                                                                                       <  6:103198/101‑1 (MQ=255)
                                                                                                cccGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTcgc                                          <  7:215477/146‑1 (MQ=255)
                                                                                                cccGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTcgc                                          <  8:21287/146‑1 (MQ=255)
                                                                                                    aCCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCgg                                   <  8:46306/149‑1 (MQ=255)
                                                                                                      cACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCGGCt                                 >  8:159637/1‑149 (MQ=255)
                                                                                                      cACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCGGCt                                 >  1:143302/1‑149 (MQ=255)
                                                                                                        cGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCGGCTTg                               >  5:50900/1‑149 (MQ=255)
                                                                                                            cAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCg                                                       <  4:158210/121‑1 (MQ=255)
                                                                                                            cAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCg                                                       >  3:158210/1‑121 (MQ=255)
                                                                                                                     gTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCGGCTTGCCGTTGTCCCAgg                  >  5:169322/1‑149 (MQ=255)
                                                                                                                                     aCGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCGAGCGGAATGTCGGCGTTTAGGGTTTCGCGGATCGGCTTGCCGTTGTCCCAGGTTTCGGTATTGGCcag  <  4:21885/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
TTCGGCAGGCTTGAGCACCACGCAGTTACCAGCTGCCAGTGCCGGGGCAAGCTTCCAGGCAGCCATCAGGATCGGGAAGTTCCAGGGGATGATCTGCCCGACCACGCCCAGTGGCTCGTGGATGTGGTAGGCCACGGTGCCTTCATTGATTTCGGCAGCGCCGCCTTCCTGGGCGCGGATGCAGCCAGCGAAATAGCGGAAGTGGTCGACCGCCAGCGGAATGTCGGCGTTGAGGGTTTCGCGGATCGGCTTGCCGTTGTCCCAGGTTTCGGTAATGGCCAG  >  NC_002947/632868‑633149

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