Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 633,281 C→T 9.7% A55T (GCC→ACC)  aldB‑I ← aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947633,2810CT 9.7% 107.5 / 3.5 41A55T (GCC→ACC) aldB‑Ialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (14/23);  new base T (2/2);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

CCAGGTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCGCGCCTTCGGTACCGGGA  >  NC_002947/633132‑633429
                                                                                                                                                     |                                                                                                                                                    
ccAGGTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGc                                                                                                                                                      >  2:1139632/1‑150 (MQ=255)
ccAGGTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGc                                                                                                                                                      >  1:357741/1‑150 (MQ=255)
ccAGGTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGCCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGc                                                                                                                                                      >  2:168501/1‑150 (MQ=255)
  aGGTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAg                                                                                                                                                    <  1:1186216/150‑1 (MQ=255)
    gTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTa                                                                                                                                                  <  1:529559/150‑1 (MQ=255)
    gTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTa                                                                                                                                                  <  1:162944/150‑1 (MQ=255)
         ggTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCCTCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCg                                                                                                                                             >  1:1036167/1‑150 (MQ=255)
               ggCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGaa                                                                                                                                       <  2:302414/150‑1 (MQ=255)
               ggCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGaa                                                                                                                                       <  1:395178/150‑1 (MQ=255)
                     cagTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGc                                                                                                                                 <  2:474080/150‑1 (MQ=255)
                          cGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGTAGTAGAGCGGGGGAACTCAGCGATCa                                                                                                                            >  1:96032/1‑150 (MQ=255)
                               ttCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGTAGTAGAGCGGGGGAACTCAGCGATCAGCTTg                                                                                                                       <  2:96032/150‑1 (MQ=255)
                                          gCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCgg                                                                                                            <  1:234154/150‑1 (MQ=255)
                                             gTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCCGCGATCAGCTTGCCATTCACCGGGGa                                                                                                         <  1:817914/150‑1 (MQ=255)
                                                   gATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTAtt                                                                                                   >  1:679995/1‑150 (MQ=255)
                                                      cTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTc                                                                                                <  1:1139632/150‑1 (MQ=255)
                                                           gcggcaCGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGt                                                                                           >  1:1203864/4‑150 (MQ=255)
                                                                      gAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGaa                                                                                             <  1:363315/137‑1 (MQ=255)
                                                                      gAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGaa                                                                                             >  2:363315/1‑137 (MQ=255)
                                                                        aCGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTAc                                                                              <  1:1128922/150‑1 (MQ=255)
                                                                                  aCTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGa                                                                                                         <  2:236225/113‑1 (MQ=255)
                                                                                  aCTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGa                                                                                                         >  1:236225/1‑113 (MQ=255)
                                                                                    tGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGc                                                                                      <  1:138438/130‑1 (MQ=255)
                                                                                    tGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGc                                                                                      >  2:138438/1‑130 (MQ=255)
                                                                                             gCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAg                                                                                                                           <  1:110867/84‑1 (MQ=255)
                                                                                             gCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAg                                                                                                                           >  2:110867/1‑84 (MQ=255)
                                                                                                  aGGCGTCGGCCGCAGCATGGGCGGCATGCAGGGCTTTGTCGATGTCTTCGGCAGTAAAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCTGAGTACCGAACTCACCACCGATGt                                                    >  2:397668/1‑150 (MQ=255)
                                                                                                             gcagcaTGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAg                                         <  2:1185039/150‑1 (MQ=255)
                                                                                                             gcagcaTGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAg                                         >  2:447750/1‑150 (MQ=255)
                                                                                                              cagcaTGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGc                                        >  1:230303/1‑150 (MQ=255)
                                                                                                                    gggCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCCGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGctct                                  <  2:357741/150‑1 (MQ=255)
                                                                                                                        ggcATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGaa                              <  1:981737/150‑1 (MQ=255)
                                                                                                                           aTCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGa                           <  2:1203864/150‑1 (MQ=255)
                                                                                                                            tCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGaa                          <  1:168501/150‑1 (MQ=255)
                                                                                                                               aGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAATcc                       <  2:230303/150‑1 (MQ=255)
                                                                                                                               aGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAAcc                       <  2:679995/150‑1 (MQ=255)
                                                                                                                                   cTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTcg                   >  1:358106/1‑150 (MQ=255)
                                                                                                                                      tttctaTGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCgcgc                <  2:1036167/145‑1 (MQ=255)
                                                                                                                                          gATGTCTTCGGTAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGaa                                                                                             >  2:204702/1‑69 (MQ=255)
                                                                                                                                          gATGTCTTCGGTAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGaa                                                                                             <  1:204702/69‑1 (MQ=255)
                                                                                                                                           aTGTCTTTGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCGCGCCTTCg           <  2:305746/150‑1 (MQ=255)
                                                                                                                                               cTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCGCGCCTTCGGTAc       >  2:1204156/1‑150 (MQ=255)
                                                                                                                                                ttCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCGCGCCTTCGGTAcc      >  2:308821/1‑150 (MQ=255)
                                                                                                                                                    gCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCTATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCGCGCCTTCGGTACCGGGa  >  1:200078/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGGTTTCGGTAATGGCCAGCAGTTCGAGGTTCTGCTCGATGCGGTCGGCGATCTTCAGCAGCACGTTGGAACGATCCTGCACTGAAGTGCGGCCCCAGGCGTCGGCCGCAGCATGGGCGGCATCCAGGGCTTTGTCGATGTCTTCGGCAGTAGAGCGGGGGAACTCAGCGATCAGCTTGCCATTCACCGGGGAGGTATTTTCGAAGTACTGCCCCTTTACCGGAGTAACGAACTCACCACCGATGTAGTTGCCGTAGCGGCTCTTGAAGGAAACCTTCGCGCCTTCGGTACCGGGA  >  NC_002947/633132‑633429

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