Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,226 C→T 11.0% intergenic (+32/+35) 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,2260CT11.0% 114.4 / 4.4 45intergenic (+32/+35)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (3/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCA  >  NC_002947/6129083‑6129373
                                                                                                                                               |                                                                                                                                                   
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                                                                                                                                   <  1:528314/136‑1 (MQ=39)
                          cgcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCg                                                                                                                                   >  2:528314/1‑136 (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                                                                                                                  >  2:235733/1‑54 (MQ=255)
                                                                                                                             tCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGt                                                                                                                  <  1:235733/54‑1 (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)
                                                                                                                                               |                                                                                                                                                   
CCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCA  >  NC_002947/6129083‑6129373

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