Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 572,528 C→T 14.3% V189M (GTG→ATG)  PP_0486 ← GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947572,5280CT14.3% 79.3 / 6.5 35V189M (GTG→ATG) PP_0486GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (13/17);  new base T (3/2);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

GCCTCCGGGATGAGCCCGTGTTCGTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGAGTACCTGCCCCAAGTCCGA  >  NC_002947/572386‑572675
                                                                                                                                              |                                                                                                                                                   
gCCTCCGGGATGAGCCCGTGTTCGTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGGTGGCCATCagcagc                                                                                                                                              <  1:508110/150‑1 (MQ=255)
   tCCGGGATGAGCCCGTGTTCGTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGGTGGCCATCagcagc                                                                                                                                              <  2:129187/147‑1 (MQ=255)
   tCCGGGATGAGCCCGTGTTCGTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGGTGGCCATCagcagc                                                                                                                                              >  1:129187/1‑147 (MQ=255)
       ggATGAGCCCGTGTTCGTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGtt                                                                                                                                       <  1:227603/150‑1 (MQ=255)
                       gTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAg                                                                                                                       <  1:812287/150‑1 (MQ=255)
                         ccATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGta                                                                                                                     >  1:446600/1‑150 (MQ=255)
                         ccATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGACCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGta                                                                                                                     <  1:119127/150‑1 (MQ=255)
                                  cgccgcGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACctgct                                                                                                            <  1:808489/150‑1 (MQ=255)
                                          aGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGa                                                                                                    <  2:77182/150‑1 (MQ=255)
                                                  gCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAgg                                                                                            <  1:186351/150‑1 (MQ=255)
                                                    ccAGGCTGTGGGGCATGGATATGCTGCCTGGTTAGCAGAGCGACTCGCACAGCACCACATCGTCACTATTCGACAGGCTGCTGCTGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGca                                                                                          >  2:54111/1‑150 (MQ=255)
                                                    ccAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGca                                                                                          <  1:110624/150‑1 (MQ=255)
                                                       cagcctgGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGcagcc                                                                                       >  2:814695/6‑150 (MQ=255)
                                                            tgGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAgcgc                                                                                  <  2:206021/150‑1 (MQ=255)
                                                                ccAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTa                                                                              >  2:232276/1‑150 (MQ=255)
                                                                            tgCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCAc                                                                  >  1:411998/1‑150 (MQ=255)
                                                                                 gggTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCGGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCt                                                             >  1:152431/1‑150 (MQ=255)
                                                                                      gCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGGGGTCAGtactac                                                                                                                 >  1:270912/1‑93 (MQ=255)
                                                                                      gCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGtactac                                                                                                                 <  2:270912/93‑1 (MQ=255)
                                                                                         gagCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGcagc                                                                                        <  2:350763/115‑1 (MQ=255)
                                                                                         gagCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGcagc                                                                                        >  1:350763/1‑115 (MQ=255)
                                                                                                       gCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCa                                                                                              <  2:226309/95‑1 (MQ=255)
                                                                                                       gCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCa                                                                                              >  1:226309/1‑95 (MQ=255)
                                                                                                              aTCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTgc                                <  2:295381/150‑2 (MQ=255)
                                                                                                                 gTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTcac                             <  2:446600/150‑1 (MQ=255)
                                                                                                                  tCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTcaca                            >  2:774260/1‑150 (MQ=255)
                                                                                                                    aCTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAg                                                                                             >  1:198536/1‑83 (MQ=255)
                                                                                                                    aCTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAg                                                                                             <  2:198536/83‑1 (MQ=255)
                                                                                                                        ttCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCg                      >  2:655808/1‑150 (MQ=255)
                                                                                                                         tCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCGGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGa                     <  2:152431/150‑1 (MQ=255)
                                                                                                                             cAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCttgtttgt                                                      >  2:176479/1‑113 (MQ=255)
                                                                                                                             cAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCttgtttgt                                                      <  1:176479/113‑1 (MQ=255)
                                                                                                                                gctgctGGTGGCCATCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGAGTACCTg              >  1:723840/1‑150 (MQ=255)
                                                                                                                                  tgctgGTGGCCATCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGAGTACCTGcc            >  2:309019/1‑150 (MQ=255)
                                                                                                                                      gATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGAGTACCTGCCCCaa        <  1:232276/150‑1 (MQ=255)
                                                                                                                                            caccaGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGAGTACCTGCCCCAAGTCCGa  <  1:774260/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCCTCCGGGATGAGCCCGTGTTCGTCCATGGCGACGCCGCGCAGGCGCCGGCCCAGGCTGTGGGCCAGGGATATGATGCCCGGGTAGCAGAGCGACTCGCACAGCACCACATCGTCACTCTTCGACAGGCTGCTGATGGCCACCAGCAGCCCATGTTGGGCGCCAGCGGTCAGTACTACCTGCTGCCAGCGAGCATCAGGCAGCCAGCGCTGTATCCAGGTCGCACCGGCTTGTTTGTGGGACGTGTGACCACCCTCTGTCACATAGTCGAGTACCTGCCCCAAGTCCGA  >  NC_002947/572386‑572675

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: