Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,521,919 C→T 10.3% A7V (GCA→GTA)  PP_4855 → osmotically‑inducible lipoprotein OsmE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,521,9190CT10.3% 96.7 / 3.7 39A7V (GCA→GTA) PP_4855osmotically‑inducible lipoprotein OsmE
Reads supporting (aligned to +/- strand):  ref base C (21/14);  new base T (2/2);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

CTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTGGCGGCAGCCCGT  >  NC_002947/5521771‑5522068
                                                                                                                                                    |                                                                                                                                                     
ctgcATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCa                                                                                                                                                      >  1:512317/1‑150 (MQ=255)
        ccGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCt                                                                                                                                              <  1:36472/150‑1 (MQ=255)
                          tACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGgc                                                                                                                            >  1:350853/1‑150 (MQ=255)
                                  gCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCg                                                                                                                    >  1:458734/1‑150 (MQ=255)
                                      acacTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAg                                                                                                                >  1:423927/1‑150 (MQ=255)
                                          tACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAgccgc                                                                                                            >  2:516578/1‑150 (MQ=255)
                                                     tcatGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAAccc                                                                                                 >  1:538816/1‑150 (MQ=255)
                                                          aCCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTgg                                                                                            >  1:50311/1‑150 (MQ=255)
                                                            cATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTg                                                                                                                              >  2:96105/1‑114 (MQ=255)
                                                            cATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTg                                                                                                                              <  1:96105/114‑1 (MQ=255)
                                                                  aaTCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTc                                                                                    >  1:434089/1‑150 (MQ=255)
                                                                                   gcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGcc                                                                   <  2:350853/150‑1 (MQ=255)
                                                                                            ttCAAAACCACACAGACCCTTTGAAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCCGAACACCCGGTGGATTACGTCACCGACCGCGATAAGCCGCTGGTCaa                                                          <  2:697445/150‑1 (MQ=255)
                                                                                             tCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt                                                                        <  2:157279/135‑1 (MQ=255)
                                                                                             tCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt                                                                        >  1:157279/1‑135 (MQ=255)
                                                                                              aaaaaCCACAAAGACCCTTTGCAAGGAGCGGGACCATGTATCAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCCCCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGTTTACGTCACCTACCGAGATGTGCCGCTGGTCAAGc                                                        <  2:847583/149‑1 (MQ=255)
                                                                                               aaaaCCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGATGGCCAAGCa                                                       >  2:700768/1‑150 (MQ=255)
                                                                                                aaaCCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCcagcag                                                      >  2:282404/1‑150 (MQ=255)
                                                                                                 aaCCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAgc                                                     <  2:585986/150‑2 (MQ=255)
                                                                                                             ccTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTa                                         <  2:538816/150‑1 (MQ=255)
                                                                                                             ccTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTa                                         >  2:305423/1‑150 (MQ=255)
                                                                                                              cTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTAt                                        <  1:513742/150‑1 (MQ=255)
                                                                                                               tttGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGTAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGc                                                                 >  1:803849/1‑124 (MQ=255)
                                                                                                               tttGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGTAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGc                                                                 <  2:803847/124‑1 (MQ=255)
                                                                                                               tttGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGTAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGc                                                                 <  2:803849/124‑1 (MQ=255)
                                                                                                               tttGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGTAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGc                                                                 >  1:803847/1‑124 (MQ=255)
                                                                                                                ttGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGa                                      <  1:516578/150‑1 (MQ=255)
                                                                                                                          cGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGaa                            <  2:423927/150‑1 (MQ=255)
                                                                                                                          cGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGaa                            >  2:569202/1‑150 (MQ=255)
                                                                                                                                 aTGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATc                     >  2:631913/1‑150 (MQ=255)
                                                                                                                                   gTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGc                                                                    >  2:509494/1‑101 (MQ=255)
                                                                                                                                   gTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGc                                                                    <  1:509494/101‑1 (MQ=255)
                                                                                                                                         gCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTgg             <  2:512317/150‑1 (MQ=255)
                                                                                                                                         gCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTgg             >  2:43204/1‑150 (MQ=255)
                                                                                                                                             aCCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAg                                                      >  2:533844/1‑105 (MQ=255)
                                                                                                                                             aCCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAg                                                      <  1:533844/105‑1 (MQ=255)
                                                                                                                                               ccTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTGGCGGCAg       >  2:131363/1‑150 (MQ=255)
                                                                                                                                                 tGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTGGCGGCAGcc     >  2:802553/1‑150 (MQ=255)
                                                                                                                                                   gCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGAAAAACCCGGTGAATTACGTCACCTCACGCAGTAATACGCTAGGCAAGCAGGTTGAGTACGTTATGACCATGCAGAAGGTCATCGCCATTGGCGGCAGCCCg   >  2:736451/1‑150 (MQ=255)
                                                                                                                                                    cAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGTCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTGGCGGCAGCCCGt  >  2:118613/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAAGCAGGTGGAGAACGGTATGACCATGCAGAAGGTCATCGCCATTGGCGGCAGCCCGT  >  NC_002947/5521771‑5522068

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