Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,039 G→T 15.4% G60* (GGA→TGA)  PP_5376 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,129,0390GT15.4% 61.6 / 3.5 26G60* (GGA→TGA) PP_5376hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/14);  new base T (2/2);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

TGATTCTGATGGCTTCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGC  >  NC_002947/6128893‑6129185
                                                                                                                                                  |                                                                                                                                                  
tGATTCTGATGGCTTCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGaa                                                                                                                                                 <  2:509868/150‑1 (MQ=255)
        aTGGCTTCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTa                                                                                                                                         <  1:370513/150‑1 (MQ=255)
             ttCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGt                                                                                                                                    <  2:191000/150‑1 (MQ=255)
              tCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGtt                                                                                                                                   <  1:295331/150‑1 (MQ=255)
              tCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGtt                                                                                                                                   <  2:522510/150‑1 (MQ=255)
                        cgATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCAGTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGa                                                                                                                         <  2:89703/150‑1 (MQ=255)
                               tGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCtgt                                                                                                                  <  2:506530/150‑1 (MQ=255)
                                  aGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGAc                                                                                                               <  1:299890/150‑1 (MQ=255)
                                       gATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGa                                                                                                          >  1:505640/1‑150 (MQ=255)
                                                    gCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCtg                                                                                             <  1:275818/150‑1 (MQ=255)
                                                       cGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGa                                                                                          >  2:452358/1‑150 (MQ=255)
                                                                              aaCAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCt                                                                   >  2:94217/1‑150 (MQ=255)
                                                                                   gCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGaa                                                              >  1:27081/1‑150 (MQ=255)
                                                                                      aTGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTg                                                                                                                      <  1:390420/91‑1 (MQ=255)
                                                                                      aTGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTg                                                                                                                      >  2:390420/1‑91 (MQ=255)
                                                                                        gcgcgcATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGCCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTc                                                         <  2:160200/150‑1 (MQ=255)
                                                                                                        tCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACTGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATTAAGACTCGACAg                                                    >  2:130852/1‑139 (MQ=255)
                                                                                                        tCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACTGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATTAAGACTCGACAg                                                    <  1:130844/139‑1 (MQ=255)
                                                                                                        tCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACTGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATTAAGACTCGACAg                                                    >  2:130844/1‑139 (MQ=255)
                                                                                                        tCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACTGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATTAAGACTCGACAg                                                    <  1:130852/139‑1 (MQ=255)
                                                                                                                             gAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCt                    <  1:452358/150‑1 (MQ=255)
                                                                                                                                 ttGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGGCGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGAGCAGCTGGCa                >  1:68633/1‑150 (MQ=255)
                                                                                                                                      cACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAAt           <  2:505640/150‑1 (MQ=255)
                                                                                                                                        cTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGa                                                               <  2:414353/96‑1 (MQ=255)
                                                                                                                                        cTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGa                                                               >  1:414353/1‑96 (MQ=255)
                                                                                                                                               tACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGc  >  2:313233/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGATTCTGATGGCTTCCTCGGCCGCGATGGCTGAAGGTGGATCTGACCGCCTGCACGGGAAGATGATTCAGGCGAATGAACAGGCCATGCGCGCATACGCAGCCTCCATCGGTAAAAATCCGCCTGAGGTTGTCCACTACCGCTACGGAATGGATCTAGATGTTGCAAAAGTGATTGCTGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGC  >  NC_002947/6128893‑6129185

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: