Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,645,712 C→A 15.4% E110* (GAA→TAA)  PP_3212 ← Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,645,7120CA15.4% 98.2 / 7.2 39E110* (GAA→TAA) PP_3212Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base A (3/3);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

AAGGAGTTCTCCATGACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAG  >  NC_002947/3645564‑3645840
                                                                                                                                                    |                                                                                                                                
aaGGAGTTCTCCATGACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTcc                                                                                                                                 >  1:715541/1‑150 (MQ=255)
   gAGTTCTCCATGACCCGCAGCTGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTgg                                                                                                                              <  1:1030330/150‑1 (MQ=255)
     gTTCTCCATGACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCa                                                                                                                            <  2:465376/150‑1 (MQ=255)
          ccATGACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGAt                                                                                                                       >  2:130774/1‑150 (MQ=255)
            aTGACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATaa                                                                                                                     <  1:130774/150‑1 (MQ=255)
              gACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGc                                                                                                                   >  2:526215/1‑150 (MQ=255)
                                 ttCCACTGTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCgga                                                                                                <  2:449149/150‑1 (MQ=255)
                                  tCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCggag                                                                                               <  1:1216376/150‑1 (MQ=255)
                                  tCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCggag                                                                                               <  1:676323/150‑1 (MQ=255)
                                                 aCTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTgg                                                                                <  1:822865/150‑1 (MQ=255)
                                                   tGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGa                                                                              >  2:1105101/1‑150 (MQ=255)
                                                     gaggATGCGTGGGGAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGAct                                                                            >  2:922062/4‑150 (MQ=255)
                                                     gAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCTGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGAct                                                                            >  2:827719/1‑150 (MQ=255)
                                                        gATGCGTCGGTAACTGGGGTCGCTGTCTGCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTg                                                                         >  2:151008/1‑150 (MQ=255)
                                                         aTGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGc                                                                        <  1:526215/150‑1 (MQ=255)
                                                                 gTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCtt                                                                >  2:905018/1‑150 (MQ=255)
                                                                  tAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTc                                                               <  1:827719/150‑1 (MQ=255)
                                                                   aaCCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTca                                                              <  2:715541/150‑1 (MQ=255)
                                                                   aaCCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTAATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTca                                                              >  2:1138116/1‑150 (MQ=255)
                                                                           gcgcCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTg                                                      >  1:83146/1‑150 (MQ=255)
                                                                            cgcCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGc                                                     <  1:1151890/150‑1 (MQ=255)
                                                                              ccGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCg                                                   >  2:188393/1‑150 (MQ=255)
                                                                                             tCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTcc                                    >  1:227405/1‑150 (MQ=255)
                                                                                               aGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTACTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGAt                                                                   <  2:1178657/117‑1 (MQ=255)
                                                                                               aGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTACTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGAt                                                                   <  2:1179001/117‑1 (MQ=255)
                                                                                               aGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTACTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGAt                                                                   >  1:1179001/1‑117 (MQ=255)
                                                                                               aGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTACTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGAt                                                                   >  1:1178657/1‑117 (MQ=255)
                                                                                                   aTCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGAGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCg                              <  1:1105101/150‑1 (MQ=255)
                                                                                                        ttGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTACTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATa                         >  1:396594/1‑150 (MQ=255)
                                                                                                                     tcCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGt                             <  1:1161446/133‑1 (MQ=255)
                                                                                                                     tcCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGt                             >  2:1161446/1‑133 (MQ=255)
                                                                                                                       cAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGt                             >  1:6645/1‑131 (MQ=255)
                                                                                                                       cAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGt                             <  2:6645/131‑1 (MQ=255)
                                                                                                                              aCCCAGGCATAGCCGTACTTTTACTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCa   <  2:396594/150‑1 (MQ=255)
                                                                                                                                ccAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAg  <  1:272059/149‑1 (MQ=255)
                                                                                                                                ccAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGGCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAg  >  2:272059/1‑149 (MQ=255)
                                                                                                                                        aGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATAccc                      <  2:417975/121‑1 (MQ=255)
                                                                                                                                        aGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATAccc                      >  1:417975/1‑121 (MQ=255)
                                                                                                                                                ttttCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCCCGCACGTGCCGGTGCAGTCGTATTcc                                    >  2:38293/1‑99 (MQ=255)
                                                                                                                                                ttttCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTcc                                    <  1:38293/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                
AAGGAGTTCTCCATGACCCGCAGCGGGCTGGTCTTCCACTCTTCGTAAAACTGGAAGATGCGTCGGTAACCGGGGGCGCCGTCTTCCGGGAACTCAGGGATCGGTTGTAGCGGGTCGTCCAGGCTGACCCAGGCATAGCCGTACTTTTCCTGGCAGTGATAAGCCTTCACCGCGGCCCCCGGAGGAATAGCGCCGGTGGGACTCTGCGGGATCTTCACGCACGTGCCGGTGCAGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAG  >  NC_002947/3645564‑3645840

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: ‑