Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,326,898 G→T 12.9% V536V (GTG→GTT mqo‑III → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,326,8980GT12.9% 79.9 / 3.8 32V536V (GTG→GTTmqo‑IIImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base T (2/2);  total (15/17)
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.68e-01

CAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGAACACCACCCGGTCGGC  >  NC_002947/3326749‑3327027
                                                                                                                                                     |                                                                                                                                 
cAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTg                                                                                                                                   >  2:390163/1‑150 (MQ=255)
      atcCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGc                                                                                                                             >  2:27581/1‑150 (MQ=255)
               tATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGCAGCTGCTGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCa                                                                                                                    >  2:425770/1‑150 (MQ=255)
                         gCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCaa                                                                                                          >  1:715829/1‑150 (MQ=255)
                                cGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGgcg                                                                                                   <  1:390163/150‑1 (MQ=255)
                                     accCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTTGCCAGCCGGTAATCATTGcc                                                                                                               >  1:270201/1‑133 (MQ=255)
                                     accCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTTGCCAGCCGGTAATCATTGcc                                                                                                               <  2:270201/133‑1 (MQ=255)
                                       ccGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCtt                                                                                            >  2:558989/1‑150 (MQ=255)
                                         ggATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTgg                                                                                          <  1:27581/150‑2 (MQ=255)
                                                    agaAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATg                                                                               >  1:705077/1‑150 (MQ=255)
                                                        ggAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCAGGATGTGGATGGCGGTGTCGTGTTCGGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCgt                                                                           >  1:181408/1‑150 (MQ=255)
                                                           gTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAAt                                                                                               >  2:726145/1‑127 (MQ=255)
                                                           gTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAAt                                                                                               <  1:726145/127‑1 (MQ=255)
                                                                 atatACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAagcca                                                                  <  2:738645/150‑1 (MQ=255)
                                                                       cAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTg                                                            >  2:196582/1‑150 (MQ=255)
                                                                        aaCGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAagccagcgcaag                                                           <  1:425770/150‑6 (MQ=255)
                                                                            aGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAg                                                       <  2:745848/150‑1 (MQ=255)
                                                                                gctgcAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGc                                                                                                                                 >  1:288032/1‑72 (MQ=255)
                                                                                gctgcAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGc                                                                                                                                 <  2:288032/72‑1 (MQ=255)
                                                                                  tgcAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCAtg                                                 >  1:628020/1‑150 (MQ=255)
                                                                                         tGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCgg                                          <  2:715829/150‑1 (MQ=255)
                                                                                         tGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCgg                                          <  2:451434/150‑1 (MQ=255)
                                                                                          gAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCggg                                         <  1:57353/150‑1 (MQ=255)
                                                                                                    tCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTTGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAgg                               >  1:195166/1‑150 (MQ=255)
                                                                                                        tCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGa                                                        <  2:281660/121‑1 (MQ=255)
                                                                                                        tCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGa                                                        >  1:281660/1‑121 (MQ=255)
                                                                                                           gCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGc                        <  2:628020/150‑1 (MQ=255)
                                                                                                                      aGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGAAcacc             <  2:181408/150‑1 (MQ=255)
                                                                                                                        tCGATAAGGCCGGCTGACAGCCAGGGAGTTGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGAAcaccac           <  2:195166/150‑1 (MQ=255)
                                                                                                                                gCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTAAAATCTGCGGCCATCGAACACCACCCGGTCgg   >  2:582240/1‑150 (MQ=255)
                                                                                                                                 ccGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGAACACCACCCGGTCGGc  <  1:558989/150‑1 (MQ=255)
                                                                                                                                            ccAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGaa                 >  2:591582/1‑124 (MQ=255)
                                                                                                                                            ccAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGaa                 <  1:591582/124‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                 
CAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGCCTGTCAGTCCGCAGCGCCCTGAGTCGATAAGGCCGGCTGACAGCCAGGGAGTGGCCAGCCGGTAATCATTGCCGTCAAGCTGGCGCAATCTTGCCACGTCCGATGCCGTAGTAAGCCAGCCCTGCCGAGGCCATGTGCTCGGGTTCGTACAGGTTGCGGCCATCGAACACCACCCGGTCGGC  >  NC_002947/3326749‑3327027

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