Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,019,126 C→A 9.1% L120M (CTG→ATG) ‡ PP_1797 → HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,019,1260CA 9.1% 134.2 / 3.0 44L120M (CTG→ATG) ‡PP_1797HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base A (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

TCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACC  >  NC_002947/2018986‑2019273
                                                                                                                                            |                                                                                                                                                   
tCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcct                                                                                                                                            <  1:429984/150‑2 (MQ=255)
 cGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGa                                                                                                                                           <  1:314261/150‑1 (MQ=255)
             gCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                               <  1:569890/150‑1 (MQ=255)
               aTGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTgg                                                                                                                             >  1:612183/1‑150 (MQ=255)
                 gCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGAc                                                                                                                           >  1:783941/1‑150 (MQ=255)
                   ccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTg                                                                                                                              >  2:245859/1‑145 (MQ=255)
                   ccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTg                                                                                                                              <  1:245859/145‑1 (MQ=255)
                   ccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTc                                                                                                                         >  2:756940/1‑150 (MQ=255)
                      cgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATCAGTCGCAGCTCAGCGCCGAGTATGAGTTCCTCTTCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGAt                                                                                                                      <  1:377265/150‑1 (MQ=255)
                                               gATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTAcc                                                                                             >  2:357522/1‑150 (MQ=255)
                                                     cTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCt                                                                                       >  2:545918/1‑150 (MQ=255)
                                                         ttGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACAAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGAcc                                                                                   <  1:45428/150‑1 (MQ=255)
                                                               tGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCgcgcgc                                                                             >  1:56749/1‑150 (MQ=255)
                                                                     aGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaagaag                                                                       >  1:389281/1‑150 (MQ=255)
                                                                       tCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAgcg                                                                     <  2:56749/150‑1 (MQ=255)
                                                                            gCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                   >  2:485823/1‑115 (MQ=255)
                                                                            gCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                   <  1:485823/115‑1 (MQ=255)
                                                                                aGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                            <  2:385324/150‑1 (MQ=255)
                                                                                    ccGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagc                                                        <  2:612183/150‑1 (MQ=255)
                                                                                         tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                   <  1:756940/150‑1 (MQ=255)
                                                                                         tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGc                                                                                                                                              >  1:52041/1‑59 (MQ=255)
                                                                                         tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGc                                                                                                                                              <  2:52041/59‑1 (MQ=255)
                                                                                           tGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGgt                                                 >  1:617338/1‑150 (MQ=255)
                                                                                                ccACCCGCAACCAACTGATTGTTGCCGGGTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCa                                            >  1:107855/1‑150 (MQ=255)
                                                                                                    ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaagaa                                                                        <  1:158441/118‑1 (MQ=255)
                                                                                                    ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaagaa                                                                        >  2:158441/1‑118 (MQ=255)
                                                                                                    ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGc                                        >  2:482932/1‑150 (MQ=255)
                                                                                                     cGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGcc                                       >  2:341433/1‑150 (MQ=255)
                                                                                                                 aTTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCg                           >  2:342563/1‑150 (MQ=255)
                                                                                                                     ttgCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGc                       >  1:505138/1‑150 (MQ=255)
                                                                                                                       gCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCTGCCACGATTCGTTGCTg                     <  1:341433/150‑1 (MQ=255)
                                                                                                                         ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAgg                                         <  1:141890/128‑1 (MQ=255)
                                                                                                                         ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAgg                                         >  2:141890/1‑128 (MQ=255)
                                                                                                                           gTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCGCGCCACGATTCGTTGCTGGGCg                 >  1:360169/1‑150 (MQ=255)
                                                                                                                              ataAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGaaa              >  2:782631/1‑150 (MQ=255)
                                                                                                                                aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGCTCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCGCGCCACGATTCGTTGCTGGGCGAAAtc            >  1:368840/1‑150 (MQ=255)
                                                                                                                                aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                   <  2:4285/111‑1 (MQ=255)
                                                                                                                                aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                   >  1:4285/1‑111 (MQ=255)
                                                                                                                                   tcggCACGCAGGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         >  1:695530/3‑150 (MQ=255)
                                                                                                                                   gagGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         <  1:482932/150‑1 (MQ=255)
                                                                                                                                   gagGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         <  1:357522/150‑1 (MQ=255)
                                                                                                                                    agGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCgg        <  1:782631/150‑1 (MQ=255)
                                                                                                                                         cGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAAc   <  1:343507/150‑1 (MQ=255)
                                                                                                                                          gCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAAcc  >  2:463127/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACC  >  NC_002947/2018986‑2019273

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: