Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,990,126 C→T 10.8% Q123* (CAG→TAG)  PP_3516 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,990,1260CT10.8% 92.9 / 3.1 37Q123* (CAG→TAG) PP_3516AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/2);  total (19/18)
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.74e-01

CGGCCAGGACCATGAAGCGTTCTACCTGGTGACCGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCC  >  NC_002947/3989978‑3990270
                                                                                                                                                    |                                                                                                                                                
cGGCCAGGACCATGAAGCGTTCTACCTGGTGACCGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCa                                                                                                                                                 <  1:453574/150‑1 (MQ=255)
             gAAGCGTTCTACCTGGTGACCGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTgg                                                                                                                                    <  2:116851/150‑1 (MQ=255)
                  gTTCTACCTGGTGACCGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCt                                                                                                                               >  2:316182/1‑150 (MQ=255)
                        ccTGGTGACCGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCt                                                                                                                         <  1:251330/150‑1 (MQ=255)
                                 cGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCg                                                                                                                >  1:462486/1‑150 (MQ=255)
                                         aTGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCa                                                                                                        >  2:21789/1‑150 (MQ=255)
                                            gCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAg                                                                                                     >  1:91231/1‑150 (MQ=255)
                                                         cTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTggg                                                                                        <  2:462486/150‑1 (MQ=255)
                                                         cTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGCGTTGTTTCATCGTAGAGCGCGGCGATGCCGGTTCCTTCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTggg                                                                                        <  1:391529/150‑1 (MQ=255)
                                                                  ggACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACtt                                                                                                                                        <  2:364177/93‑1 (MQ=255)
                                                                  ggACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACtt                                                                                                                                        >  1:364177/1‑93 (MQ=255)
                                                                        ccaccaGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATAc                                                                         >  2:384851/1‑150 (MQ=255)
                                                                         caccaGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCACCCTGCTGGGGCACAAGCGCTATAcc                                                                        >  1:87781/1‑150 (MQ=255)
                                                                             aGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGc                                                                    <  2:392547/150‑1 (MQ=255)
                                                                                   gTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTgg                                                                                                                                    <  1:206812/80‑1 (MQ=255)
                                                                                   gTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTgg                                                                                                                                    >  2:206812/1‑80 (MQ=255)
                                                                                                 ggTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAg                                                <  2:91231/150‑1 (MQ=255)
                                                                                                       aTCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTg                                          >  1:219805/1‑150 (MQ=255)
                                                                                                              aGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCCCAAGCGCTATACCCGGCCTTGCTTCGACGCCCAACAAGCCCCGGGGCGCa                                   >  2:411096/1‑150 (MQ=255)
                                                                                                                cgcgGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATc                                 <  1:316182/150‑1 (MQ=255)
                                                                                                                 gcggcgATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGgc                                      <  1:19009/144‑1 (MQ=255)
                                                                                                                 gcggcgATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGgc                                      >  2:19009/1‑144 (MQ=255)
                                                                                                                        tGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGa                         <  1:186379/150‑1 (MQ=255)
                                                                                                                         gCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAg                        <  2:132046/150‑1 (MQ=255)
                                                                                                                           aCCCTACCGCTTCCACTACGGCACGTAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCa                      >  2:214355/1‑150 (MQ=255)
                                                                                                                             ccTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAAc                    >  2:101061/1‑150 (MQ=255)
                                                                                                                              cTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACt                   >  2:10402/1‑150 (MQ=255)
                                                                                                                               tACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTc                  <  1:10402/150‑1 (MQ=255)
                                                                                                                               tACCGCTTCCACTACGGCACGCAGAACGACTGGTGGGTGCTCACGCTCACCGAGCGGGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCAGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTc                  >  1:395174/1‑150 (MQ=255)
                                                                                                                                   gCTTCCACTACGGCACGTAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGAc               >  2:43737/1‑149 (MQ=255)
                                                                                                                                   gCTTCCACTACGGCACGTAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGAc               <  1:43737/149‑1 (MQ=255)
                                                                                                                                    cTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCa                                   >  2:178410/1‑128 (MQ=255)
                                                                                                                                    cTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCa                                   <  1:178410/128‑1 (MQ=255)
                                                                                                                                      tCCACTACGGCACGTAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCtgt           <  1:214355/150‑1 (MQ=255)
                                                                                                                                               gCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCGACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGcc  >  1:375752/1‑150 (MQ=255)
                                                                                                                                                   gCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAgcgc                                                                              <  1:368426/70‑1 (MQ=255)
                                                                                                                                                   gCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAgcgc                                                                              >  2:368426/1‑70 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGGCCAGGACCATGAAGCGTTCTACCTGGTGACCGTGCCGCATGGCAGCGAGGTGCACTTTGAACAGGACGGCCACCAGATCAGTTGCAGCCCGGGCGGTTTCATCGTCGAGCGCGGCGATGCACCCTACCGCTTCCACTACGGCACGCAGAACGACTTGTGGGTGCTCAAGCTGCCCGAGCGTGCGCTCAAAGCCAACCTGCTGGGGCACAAGCGCTATACCCGGCATTGCTTCGACGCCCGTCAGGGCCTGGGGCGCATCTTCGTCGAGCAACTCGACCTGTGTGCCCGCC  >  NC_002947/3989978‑3990270

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