Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,723,174 C→A 9.2% L130M (CTG→ATG)  PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,723,1740CA 9.2% 132.8 / 2.1 43L130M (CTG→ATG) PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base A (2/2);  total (21/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 = 9.95e-01

GCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATC  >  NC_002947/3723026‑3723317
                                                                                                                                                    |                                                                                                                                               
gcgcAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCt                                                                                                                                                >  2:430194/1‑150 (MQ=255)
                gTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTtgctg                                                                                                                                <  1:1076619/150‑1 (MQ=255)
                        cGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAAc                                                                                                                        >  1:172914/1‑150 (MQ=255)
                         gAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGat                                                                                                                                            >  1:175853/1‑129 (MQ=255)
                         gAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGat                                                                                                                                            <  2:175853/129‑1 (MQ=255)
                           ggtggtGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAAcgcg                                                                                                                     <  2:172914/150‑1 (MQ=255)
                                               ccTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTac                                                                                                 <  2:714988/150‑1 (MQ=255)
                                               ccTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTac                                                                                                 >  2:628981/1‑150 (MQ=255)
                                                    gaaaagaTCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATAACCAATTGCTGCGCGAAACGCAAATACCCTTGCACCTGGTACACAAg                                                                                            >  1:41900/5‑150 (MQ=255)
                                                       gaaaTCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCCCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACGTCGCTGGATATNCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGgcc                                                                                         <  1:668090/147‑1 (MQ=255)
                                                             cTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTgg                                                                                                    >  1:1190443/1‑133 (MQ=255)
                                                             cTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTgg                                                                                                    <  2:1190443/133‑1 (MQ=255)
                                                              tGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCt                                                                                                             >  2:34107/1‑123 (MQ=255)
                                                              tGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCt                                                                                                             <  1:34107/123‑1 (MQ=255)
                                                              tGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATg                                                                                  >  1:547496/1‑150 (MQ=255)
                                                                      ttGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCa                                                                          >  2:88448/1‑150 (MQ=255)
                                                                             aGCAACCGTTTGCCACGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAggg                                                                   <  1:423332/150‑1 (MQ=255)
                                                                                   cGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTgg                                                             >  2:1193523/1‑150 (MQ=255)
                                                                                          cATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTgg                                                                                                    <  2:110693/104‑1 (MQ=255)
                                                                                          cATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTgg                                                                                                    >  1:110693/1‑104 (MQ=255)
                                                                                          cATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGt                                                      <  2:547496/150‑1 (MQ=255)
                                                                                              cTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGAc                                                  <  1:88448/150‑1 (MQ=255)
                                                                                               tGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATAACCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCATGTGGGGTTCTGGCGGGGGTGGAcg                                                 >  2:116702/1‑149 (MQ=255)
                                                                                               tGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGAGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGAcc                                                 >  1:800741/1‑150 (MQ=255)
                                                                                                          cTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCg                                      <  1:1193523/150‑1 (MQ=255)
                                                                                                           tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTacac                                                                                               >  1:355521/1‑92 (MQ=255)
                                                                                                           tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTacac                                                                                               <  2:355521/92‑1 (MQ=255)
                                                                                                           tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCg                                                                                <  1:253551/107‑1 (MQ=255)
                                                                                                           tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCg                                                                                >  2:253551/1‑107 (MQ=255)
                                                                                                           tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTc                                             <  1:303718/142‑1 (MQ=255)
                                                                                                           tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTc                                             >  2:303718/1‑142 (MQ=255)
                                                                                                             gAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCAt                                                                                   >  2:132608/1‑102 (MQ=255)
                                                                                                             gAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCAt                                                                                   <  1:132608/102‑1 (MQ=255)
                                                                                                                    aGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGt                            <  1:116702/150‑1 (MQ=255)
                                                                                                                         tggtggATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGaa                       <  1:628981/150‑1 (MQ=255)
                                                                                                                              gATGCGCGCAATGTTCACTTCGATGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAgg                                                                    <  2:220513/100‑1 (MQ=255)
                                                                                                                              gATGCGCGCAATGTTCACTTCGATGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAgg                                                                    <  2:220481/100‑1 (MQ=255)
                                                                                                                              gATGCGCGCAATGTTCACTTCGATGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAgg                                                                    >  1:220481/1‑100 (MQ=255)
                                                                                                                              gATGCGCGCAATGTTCACTTCGATGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAgg                                                                    >  1:220513/1‑100 (MQ=255)
                                                                                                                                tGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACaa                                                                                             <  1:429064/73‑1 (MQ=255)
                                                                                                                                tGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACaa                                                                                             >  2:429064/1‑73 (MQ=255)
                                                                                                                                          gTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATGAGTTTGAAGGCTTCAACGACCAGat      <  1:1068658/150‑1 (MQ=255)
                                                                                                                                              aCTCCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGatcatc  >  2:760758/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATC  >  NC_002947/3723026‑3723317

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