Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,153:1 +C 9.3% coding (424/519 nt) PP_2340 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,673,1531.C 9.3% 153.9 / 11.3 43coding (424/519 nt)PP_2340hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (20/19);  new base C (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

TCCGCTTTGTCGAGCAGCAGTCGGCGTATTGCTGGGCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGT  >  NC_002947/2673004‑2673299
                                                                                                                                                      |                                                                                                                                                  
tCCGCTTTGTCGAGCAGCAGTCGGCGTATTGCTGGGCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGcgatcg                                                                                                                                                     >  1:40745/1‑150 (MQ=255)
       tGTCGAGCAGCAGTCGGCGTATTGCTGGGCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGa                                                                                                                                             <  1:903385/150‑1 (MQ=255)
           gAGCAGCAGTCGGCGTATTGCTGGGCTTACCATCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGttc                                                                                                                                         <  1:786205/150‑1 (MQ=255)
              cagcagTCGGCGTATTGCTGGGCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTGTGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGttcttc                                                                                                                                      <  2:688772/150‑1 (MQ=255)
                               cTGGGCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCa                                                                                                                     <  2:692304/150‑1 (MQ=255)
                                   gCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTAc                                                                                                                 <  1:362313/150‑1 (MQ=255)
                                         cAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGttct                                                                                                                                        >  1:956835/1‑121 (MQ=255)
                                         cAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGttct                                                                                                                                        <  2:956835/121‑1 (MQ=255)
                                          aGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGAGCt                                                                                                          >  1:368963/1‑150 (MQ=255)
                                           gCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCta                                                                                                         >  1:381147/1‑150 (MQ=255)
                                                ttGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGGTTTCAGAGGCGCACTACCAGCGCTATATCg                                                                                                    >  2:536515/1‑150 (MQ=255)
                                                  gCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGcc                                                                                                  <  2:40745/150‑1 (MQ=255)
                                                      gcgGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCa                                                                                              <  2:936978/150‑1 (MQ=255)
                                                             gcgGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCt                                                                                       <  2:1045129/150‑1 (MQ=255)
                                                                     tCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATg                                                                               <  1:209783/150‑1 (MQ=255)
                                                                             cTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATAGCCCCCAGCTTCCTGCGCGATGCCATGGGc                                                                       >  1:311187/1‑150 (MQ=255)
                                                                              tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                      <  2:1101029/150‑1 (MQ=255)
                                                                                     cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCa                                                               >  1:621925/1‑150 (MQ=255)
                                                                                     cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCa                                                               >  1:929588/1‑150 (MQ=255)
                                                                                             tGGTGTTGCTGCCGGAGTTCCGTGAGTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATANCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                       <  1:270266/150‑1 (MQ=255)
                                                                                                tgtTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAg                                                    >  1:1105493/1‑150 (MQ=255)
                                                                                                    gctgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCgg                                                >  1:973041/1‑150 (MQ=255)
                                                                                                      agcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGGGAGCCGGGc                                              >  1:318327/2‑150 (MQ=255)
                                                                                                        ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGa                                            >  1:944775/1‑150 (MQ=255)
                                                                                                          ggAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                          <  2:973041/150‑1 (MQ=255)
                                                                                                          ggAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGCGCTTCCTGCTTTCACAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCGGCGCGCAGGGCCAGTGAGCCGGGCGaac                                          >  2:715843/1‑149 (MQ=255)
                                                                                                           gAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATc                                         >  2:196153/1‑150 (MQ=255)
                                                                                                              ttCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTatat                                                                                                      <  1:979576/86‑1 (MQ=255)
                                                                                                              ttCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTatat                                                                                                      >  2:979576/1‑86 (MQ=255)
                                                                                                                  gTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGcc                                  <  2:1105493/150‑1 (MQ=255)
                                                                                                                   tGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATTCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAATCGCCg                                 >  2:155065/1‑150 (MQ=255)
                                                                                                                            aCTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGATGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCg                        >  1:577606/1‑150 (MQ=255)
                                                                                                                             cTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGCGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGc                       <  1:1029457/150‑1 (MQ=255)
                                                                                                                              tGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCa                      >  2:1109680/1‑150 (MQ=255)
                                                                                                                                 cccAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCt                                                                                           >  2:586241/1‑78 (MQ=255)
                                                                                                                                 cccAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCt                                                                                           <  1:586241/78‑1 (MQ=255)
                                                                                                                                 cccAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAAcc                   >  2:416430/1‑150 (MQ=255)
                                                                                                                                        tcgcgATCGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATgg                                                                         <  1:705476/88‑1 (MQ=255)
                                                                                                                                        tcgcgATCGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATgg                                                                         >  2:705476/2‑89 (MQ=255)
                                                                                                                                        tcgcgATCGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGgt                               <  2:939648/130‑1 (MQ=255)
                                                                                                                                        tcgcgATCGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGgt                               >  1:939648/2‑131 (MQ=255)
                                                                                                                                          gcgatgcgatCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCa          >  1:616138/1‑150 (MQ=255)
                                                                                                                                           cgatgcgatCGATCGTGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg         <  2:1049718/150‑1 (MQ=255)
                                                                                                                                           cgatgcgatCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg         <  1:681189/150‑1 (MQ=255)
                                                                                                                                                  atcnatcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt  >  2:699873/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCCGCTTTGTCGAGCAGCAGTCGGCGTATTGCTGGGCTTACCAGCTGATTGCCCGCGGGGCGCGGATGATCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGT  >  NC_002947/2673004‑2673299

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: