Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,232 A→G 10.8% intergenic (+38/+29) PP_5376 → / ← PP_5377 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,129,2320AG10.8% 120.9 / 5.2 46intergenic (+38/+29)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base G (3/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGG  >  NC_002947/6129083‑6129380
                                                                                                                                                     |                                                                                                                                                    
ccGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTaaaa                                                                                                                                                      <  1:427384/150‑1 (MQ=255)
         cTGGGACGTGATGCCTGCGTGAATGGCCTATGAACAGTCGACGGGCAAGCTCAACGTTCTCTAATACCGCGCAGCTGGCATCAATTTGTCCGGCCAAAACTGATCTCTGACTGCGGTCCAATCGAAAACGACTCCTTACAATCGCCGttt                                                                                                                                             <  2:486033/150‑1 (MQ=255)
               cGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTg                                                                                                                                       <  2:998189/150‑1 (MQ=255)
                gTGATGCCTGCGCAAATGACCTATGAAAACTAACCACGAAAATCGAACATTCTCGAATACCGCGCAGCTCGCCTCAAGTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTgg                                                                                                                                      <  1:519372/150‑1 (MQ=255)
                      ccTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCAGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAAc                                                                                                                                >  2:367542/1‑150 (MQ=255)
                        tGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACgg                                                                                                                              >  1:1030396/1‑150 (MQ=255)
                         gcgcAAATGACCTATGAAGCCTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGt                                                                                                                             <  2:896418/150‑1 (MQ=255)
                          cgcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCg                                                                                                                                          >  2:528314/1‑136 (MQ=39)
                          cgcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCg                                                                                                                                          <  1:528314/136‑1 (MQ=39)
                          cgcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTg                                                                                                                            <  2:148397/150‑1 (MQ=255)
                          cgcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAAGTGCCGCGGCCAAAATTGTGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGGTTGCGTTGGGAAATCGGTg                                                                                                                            <  2:378524/150‑1 (MQ=255)
                                   ccTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGAGTTGGGAAAACGGTGGGTCAGTgg                                                                                                                   >  1:1078022/1‑150 (MQ=255)
                                    cTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGc                                                                                                                  >  1:997653/1‑150 (MQ=255)
                                    cTATGAAGACTCGACAGGCAAGCGCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGGCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGc                                                                                                                  <  2:538290/150‑1 (MQ=255)
                                      aTGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGTCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCt                                                                                                                <  1:202317/150‑1 (MQ=255)
                                          aGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGAt                                                                                                            >  2:449830/1‑150 (MQ=255)
                                             cTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCTTTTAAGGCGCCGttttc                                                                                                                                           >  1:1013346/1‑114 (MQ=255)
                                             cTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCTTTTAAGGCGCCGttttc                                                                                                                                           <  2:1013346/116‑3 (MQ=255)
                                              tCGACAGGCAAGCCCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTg                                                                                                                                       <  1:445106/119‑1 (MQ=255)
                                              tCGACAGGCAAGCCCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTg                                                                                                                                       >  2:445106/1‑119 (MQ=255)
                                                gACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTg                                                                                                      >  1:681036/1‑150 (MQ=255)
                                                   aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGtt                                                                                                                                              >  1:84389/1‑107 (MQ=255)
                                                   aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGtt                                                                                                                                              <  2:84389/107‑1 (MQ=255)
                                                     gCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGt                                                                                                                             >  1:368265/1‑122 (MQ=255)
                                                     gCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGt                                                                                                                             <  2:368265/122‑1 (MQ=255)
                                                       aaGCTCACCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCTGTTTGCGTTGGGAAAACGGTGGGTTAGTGGCCTGGATACTCTGCTACCAc                                                                                               <  2:324829/150‑1 (MQ=255)
                                                          cTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTc                                                                                            >  1:943509/1‑150 (MQ=255)
                                                           tCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTTGCCTGGATACTCTGCAACCACTTCa                                                                                           >  1:550515/1‑150 (MQ=255)
                                                                aTTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCa                                                                                                                       <  1:624520/117‑1 (MQ=255)
                                                                aTTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCa                                                                                                                       >  2:624520/1‑117 (MQ=255)
                                                                                      gCATCAATTGCCGCGGCCAAACTTGAGCTCTGAATGCTCTCCAAACGAAAACGGCGCCTAAAAAGCGCCGTTTGCGTCGGGAAAACGGTGGGTCAGTAGCCTGGATACTCTGCAACCACTTCATGCCGTCCGCTCTTGTTGAGACCGATg                                                                >  2:356886/1‑150 (MQ=255)
                                                                                            aTTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTgg                                                          <  1:442858/150‑1 (MQ=255)
                                                                                             ttGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCACTCTTTTTGAGACCGATGACCTGGt                                                         <  2:98861/150‑1 (MQ=255)
                                                                                               gCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTAt                                                       <  1:449830/150‑1 (MQ=255)
                                                                                                 cgcgGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGc                                                     <  2:250993/150‑1 (MQ=255)
                                                                                                   cgGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGAGACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCAt                                                   <  2:997653/150‑1 (MQ=255)
                                                                                                                gCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGCCCTGGTATGCAGGAGAGACACCATc                                      <  2:943509/150‑1 (MQ=255)
                                                                                                                  tctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTCGTATGCATGAGAGACACCATCCa                                    <  2:550515/150‑1 (MQ=255)
                                                                                                                             tCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGt                                                                                                                         <  1:235733/54‑1 (MQ=255)
                                                                                                                             tCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGt                                                                                                                         >  2:235733/1‑54 (MQ=255)
                                                                                                                                    gcaaaCGGCGCTTTTAAGGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGc                  >  2:857702/3‑150 (MQ=255)
                                                                                                                                     aaaaCGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATAc                           <  1:809049/140‑1 (MQ=255)
                                                                                                                                     aaaaCGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATAc                           >  2:809049/1‑140 (MQ=255)
                                                                                                                                      aaaCGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCt                >  1:545047/1‑150 (MQ=255)
                                                                                                                                             gcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCa         >  2:838060/1‑150 (MQ=255)
                                                                                                                                                 tAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                                      <  1:7462/117‑1 (MQ=255)
                                                                                                                                                 tAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                                      >  2:7462/1‑117 (MQ=255)
                                                                                                                                                    aaGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCAGACCTgg  <  2:545047/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGG  >  NC_002947/6129083‑6129380

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