Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,603,520 G→T 15.4% P393T (CCG→ACG)  PP_3175 ← ferredoxin reductase dioxygenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,603,5200GT15.4% 66.8 / 2.2 26P393T (CCG→ACG) PP_3175ferredoxin reductase dioxygenase subunit
Reads supporting (aligned to +/- strand):  ref base G (8/13);  new base T (2/2);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

GGCCATGCTGGTAGGTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAG  >  NC_002947/3603385‑3603652
                                                                                                                                       |                                                                                                                                    
ggCCATGCTGGTAGGTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCgg                                                                                                                                      >  1:482153/1‑136 (MQ=255)
  ccATGCTGGTAGGTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGTGa                                                                                                                                    >  2:317582/1‑136 (MQ=255)
  ccATGCTGGTAGGTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGTGa                                                                                                                                    >  2:181032/1‑136 (MQ=255)
  ccATGCTGGTAGGTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGCGa                                                                                                                                    >  2:319554/1‑136 (MQ=255)
             ggTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGata                                                                                                                         >  2:300625/1‑136 (MQ=255)
                    aaaaaCGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCcgg                                                                                                                  >  1:491585/1‑136 (MQ=255)
                    aaaaaCGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCcgg                                                                                                                  >  1:522360/1‑136 (MQ=255)
                           gTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCca                                                                                                           <  2:90565/136‑1 (MQ=255)
                            tACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCcac                                                                                                          <  1:445780/136‑1 (MQ=255)
                                    cATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTc                                                                                                  >  2:10762/1‑136 (MQ=255)
                                            cGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTgcgc                                                                                          >  2:94602/1‑136 (MQ=255)
                                                                                     ttGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGc                                                 <  1:334267/136‑1 (MQ=255)
                                                                                          gggCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCAc                                            <  2:389334/136‑1 (MQ=255)
                                                                                              cAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTCCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCATCGCCAAGCGGCACGcca                                        >  2:61713/1‑135 (MQ=255)
                                                                                                   gcCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAg                                   >  2:50040/1‑136 (MQ=255)
                                                                                                           tCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCgg                           <  1:161465/136‑1 (MQ=255)
                                                                                                            cGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGt                          <  1:94602/136‑1 (MQ=255)
                                                                                                            cGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGt                          <  1:48214/136‑1 (MQ=255)
                                                                                                            cGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGt                          <  2:466788/136‑1 (MQ=255)
                                                                                                                gcgcACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCt                                        >  2:274434/1‑118 (MQ=255)
                                                                                                                gcgcACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCt                                        <  1:274434/118‑1 (MQ=255)
                                                                                                                    aCCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGcc                  <  2:324571/136‑1 (MQ=255)
                                                                                                                    aCCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGcc                  <  2:467831/136‑1 (MQ=255)
                                                                                                                         gcgcCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTAcc             <  2:329526/136‑1 (MQ=255)
                                                                                                                                  gAGCGTGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCg    <  1:317582/136‑1 (MQ=255)
                                                                                                                                  gAGCGTGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCg    <  1:181032/136‑1 (MQ=255)
                                                                                                                                    gCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAg  <  2:139203/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GGCCATGCTGGTAGGTCCAGAAAAACGGTACATCGGCATAACGGCGTTGTTCGCCAAGCATGTTGGCGGCGGCGATCACCCCGTGTTGCTGGGCCAGGCGCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAG  >  NC_002947/3603385‑3603652

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