Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,308,146 G→C 10.4% D104H (GAC→CAC)  csiD → carbon starvation induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,308,1460GC10.4% 137.0 / 4.1 46D104H (GAC→CAC) csiDcarbon starvation induced protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (2/3);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACAGCGACAACTCCGACAGCTACCTGCG  >  NC_002947/3308004‑3308288
                                                                                                                                              |                                                                                                                                              
aTTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAg                                                                                                                                         >  2:1347998/1‑150 (MQ=255)
aTTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAg                                                                                                                                         >  1:475042/1‑150 (MQ=255)
    gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAg                                                                                                                                     <  2:238756/150‑1 (MQ=255)
      gtaaaaGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGCTCACGTCATCCAGGc                                                                                                                                   <  2:1348822/146‑1 (MQ=255)
      gTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGCTCACGTCATCCAGGc                                                                                                                                   <  2:642193/150‑1 (MQ=255)
      gTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGc                                                                                                                                   <  2:815174/150‑1 (MQ=255)
       tACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGCTCACGtc                                                                                                                                           <  2:47089/141‑2 (MQ=255)
       tACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGCTCACGtc                                                                                                                                           >  1:47089/1‑140 (MQ=255)
               gTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACAt                                                                                                                          >  1:138911/1‑150 (MQ=255)
                    tGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCa                                                                                                                     >  1:950453/1‑150 (MQ=255)
                             gtgtAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGACCGaa                                                                                                                               <  2:611584/131‑1 (MQ=255)
                             gtgtAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGACCGaa                                                                                                                               >  1:611584/1‑131 (MQ=255)
                                     gAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTgc                                                                                                    >  2:1355077/1‑150 (MQ=255)
                                          ccTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCa                                                                                               >  2:188728/1‑150 (MQ=255)
                                                       tgCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGt                                                                                  <  1:1067970/150‑1 (MQ=255)
                                                       tgCCAGGGCACGCTGAGTAAGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGt                                                                                  >  2:837713/1‑150 (MQ=255)
                                                            gggCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCaa                                                                             <  1:78655/150‑1 (MQ=255)
                                                                  gCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGa                                                                       <  1:1355077/150‑1 (MQ=255)
                                                                  gCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGa                                                                       >  1:1321285/1‑150 (MQ=255)
                                                                        tCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCAt                                                                 <  1:1347998/150‑1 (MQ=255)
                                                                             tgctgcTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGACGGAGCGATCATCAGGCAGTGGACTTCCGGCTCTCACCTGATCGGCCGTTCCAGCTATGACGCCATGAGTg                                                            >  1:1236569/1‑150 (MQ=255)
                                                                                  ctCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAg                                                       >  1:402023/1‑150 (MQ=255)
                                                                                   tCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAACCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGt                                                      >  1:137960/1‑150 (MQ=255)
                                                                                    cAACACCATCCTCGACCGCGCCACAGGCGGCATGCTACAAAGTGCAAATGACCTATGTCACGGGAGCCTTGCCTTAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGtt                                                     >  1:387253/1‑150 (MQ=255)
                                                                                      acacCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATg                                                                                                                         <  2:121712/80‑1 (MQ=255)
                                                                                      acacCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATg                                                                                                                         >  1:121712/1‑80 (MQ=255)
                                                                                        acCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTAt                                                 <  1:188728/150‑1 (MQ=255)
                                                                                        acCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTAt                                                 <  2:138911/150‑1 (MQ=255)
                                                                                          cATCCTCGACCGCGCGACAGGCGGTATGCTGATTACCCCGATCCGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTACACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATCACGCCATGAGTGGCCAGTTCTATgc                                               <  2:1356972/150‑1 (MQ=255)
                                                                                            tCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATgcgc                                             >  2:809781/1‑150 (MQ=255)
                                                                                            tCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATgcgc                                             <  1:996597/150‑1 (MQ=255)
                                                                                            tCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATgcgc                                             <  2:1321285/150‑1 (MQ=255)
                                                                                               tCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCt                                          <  2:950453/150‑1 (MQ=255)
                                                                                                  aCCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCg                                       >  1:855590/1‑150 (MQ=255)
                                                                                                   ccGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTACGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCgt                                      <  1:1046657/150‑1 (MQ=255)
                                                                                                      cgcCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCgtggt                                   <  2:475042/150‑1 (MQ=255)
                                                                                                        ccACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGtggtgg                                 <  2:971181/150‑1 (MQ=255)
                                                                                                            aGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCaa                             <  2:387253/150‑1 (MQ=255)
                                                                                                                        gATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACa                           >  2:1366814/1‑140 (MQ=255)
                                                                                                                        gATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACa                           <  1:1366814/140‑1 (MQ=255)
                                                                                                                              cccGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCaa                                                                             >  1:987521/1‑84 (MQ=255)
                                                                                                                              cccGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCaa                                                                             <  2:987521/84‑1 (MQ=255)
                                                                                                                              cccGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACAGCGACAACTCCGACAg           >  1:443990/1‑150 (MQ=255)
                                                                                                                                cGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACAGCGACAACTCCGACAGCt         <  2:402023/150‑1 (MQ=255)
                                                                                                                                      gCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACAGCGACAACTCCGACAGCTACCTgc   >  1:644798/1‑150 (MQ=255)
                                                                                                                                       ccTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACAGCGACAACTCCGACAGCTACCTgcg  >  1:1259559/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCATGAGTGGCCAGTTCTATGCGCGCTTCGTGGTGGTCAACAGCGACAACTCCGACAGCTACCTGCG  >  NC_002947/3308004‑3308288

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: ‑