Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,709,006 A→G 21.5% A696A (GCA→GCG rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,709,0060AG21.5% 46.3 / 14.5 42A696A (GCA→GCGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (4/4);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAG  >  NC_000913/1708860‑1709152
                                                                                                                                                  |                                                                                                                                                  
gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAACGCGGCCGTGGCAGCGGCGAGTGCACGCGCTAAGGCCCga                                                                                                                                                 >  1:214084/1‑149 (MQ=14)
    aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGc                                                                                                                                             >  1:402801/1‑150 (MQ=255)
         cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGaa                                                                                                                                        >  1:947775/1‑150 (MQ=255)
          gcAAAGCCGCCGTCTAAGCCGCTATTGCCCGATCCAAAGCGCGCAACCGGGAACAGCAACAGGCTAATGCGGATCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAAc                                                                                                                                       <  1:977172/150‑1 (MQ=14)
                     gTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt                                                                                                                            >  1:1107302/1‑150 (MQ=255)
                        gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTaaa                                                                                                                         >  2:879677/1‑149 (MQ=255)
                               cTATTTCCCGCGCCAAAGCGCTCACTCGGTAACAGCAACCGCCTATTGCGGAGCCACTAGCTCAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAc                                                                                                                  <  1:1071516/150‑1 (MQ=17)
                                      ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCCACACGCTAATGCGGTACCagcaga                                                                                                           >  1:139692/1‑150 (MQ=255)
                                             acaGCGCGCAACGGGGAACAGCAACCGGCTAATGAGGACGCAGAACACCTGGTTGATCCGCGAAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGtt                                                                                                    <  1:1188501/148‑1 (MQ=17)
                                                   cgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcg                                                                                              >  1:1055256/1‑150 (MQ=255)
                                                     cAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGACGCGGCTATTGCACGCGCAAAAGCACGCAAGCAGGACCAGCAGCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCagcc                                                                                            >  2:878263/1‑149 (MQ=37)
                                                                  cAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGtt                                                                               >  1:258391/1‑150 (MQ=255)
                                                                       ggcaactGCCGAGCCAGCACCCCAGGCCGAGCCGCGCAAACCCGCCGTCGACGCCGCTATTGCACGCGCCACAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGGACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgc                                                                          <  1:322144/144‑1 (MQ=14)
                                                                       ggCTAATTCGGACGGAGCATAGCAGGTTGATCCTCGTAAATCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTATTGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgc                                                                          <  1:696677/150‑1 (MQ=17)
                                                                       ggCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCTCAAATGCAAGAATGATGTTAAAGCATCAGGTTAATGAGGTAGCAGAAGAGCAGGTTGATCCGCGCAGAGCGGCAGTTGCcgc                                                                          >  2:661388/1‑150 (MQ=25)
                                                                               cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTg                                                                  <  1:1197025/150‑1 (MQ=255)
                                                                                ggcgCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCCAAGCACGCAAGCTGGAACAGCAACAGGCTAATGAGGTACCAGAAGACCAGGTAGCTCCGCGCAACGCGGCAGTTGCCGCGGCTATTGc                                                                 >  2:674791/4‑150 (MQ=18)
                                                                                     cAGAAGCCCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgcg                                                            <  1:391526/150‑1 (MQ=255)
                                                                                     cAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAATAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCAATTGCCcgcg                                                            >  2:101123/1‑150 (MQ=255)
                                                                                          cccgttgTTGATCCTCGTAAAGCTGCCGTAGAAGCCGCTATTGCACGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAg                                                       <  2:60705/144‑1 (MQ=18)
                                                                                             cAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGACAAAGCACGCAAGCTGGAACAGAAACAGGCTAATTCGGTACCAGATGAACAGGTTCAACCGCGCAAAGCGGCAGTTGACGCGGAAATTGCACGCGCTCAGGcc                                                    >  1:946206/1‑150 (MQ=21)
                                                                                              tggTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCa                                                   <  2:1090721/149‑1 (MQ=11)
                                                                                                  tGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaa                                               >  2:281782/1‑150 (MQ=255)
                                                                                                            aaaGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCAGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcag                                     >  1:116265/1‑150 (MQ=255)
                                                                                                              tgCTGCCGGCAAAGACGATATTGCACCTGCCTAAGAACGCAAGCTGGAACAGCAACAGGCTAATGCGGCCCCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagca                                   <  1:66908/149‑1 (MQ=18)
                                                                                                              aGCTGCCTTCTAAGCCGCTTTTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagca                                   <  1:757306/150‑1 (MQ=255)
                                                                                                              aGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGGAAGCAGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGCTCCGCGCAAAGCGGCAGTTCCCGCGGCTATTGCACGCGCTCAGGCCAAAAAAGCCGCCcagca                                   >  1:336859/1‑150 (MQ=25)
                                                                                                                cTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACCGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGa                                 <  1:869914/150‑1 (MQ=255)
                                                                                                                             ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACgt                    >  2:237556/1‑149 (MQ=255)
                                                                                                                             cgCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGCCGCGGCTATTGCCCGAGCTAAGGCCAAAAAACCCGCACAGCAGCAGGCTGTAAACGa                    >  2:1027145/2‑150 (MQ=17)
                                                                                                                              gCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAACAGAAGAACAGGTTGATCCTCTCAAAGCGGCAGTTGCCGCTGCTATTGCCCTCGCTCAGGCAAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAg                   <  1:1209455/150‑1 (MQ=25)
                                                                                                                               cTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg                  <  1:1027145/150‑1 (MQ=255)
                                                                                                                                  ttccgctcgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACt               <  1:68345/143‑1 (MQ=255)
                                                                                                                                   tGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTa              >  1:1180697/1‑150 (MQ=255)
                                                                                                                                       cgcgCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAAGGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATg          <  1:56229/150‑1 (MQ=255)
                                                                                                                                        gcgcCAAAGCACGCAAGCTGGAACCGCAACTGGCTAATTCGGTTCCAGAAGTACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATgg         <  2:214084/150‑1 (MQ=38)
                                                                                                                                           ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcgcg                                                                                            >  2:177806/1‑64 (MQ=25)
                                                                                                                                           ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcgcg                                                                                            <  1:177806/64‑1 (MQ=25)
                                                                                                                                            cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTTGCCGCGGCTATTGCCCGCGCACAGGCCAAAAAAGCCGAACAGCAGAAGGATGAAAACGAGGACTAAAAGGTAtt     >  1:892279/1‑150 (MQ=18)
                                                                                                                                            cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTAGCCGCGGCTATTGCCCGCGCGCAGGCCAACAAAGCCGCCCAGCAGCAGGGTGTAACCGAGGACGAAAAGGTagt     >  1:469875/1‑148 (MQ=17)
                                                                                                                                              aaGCCCGCACGCGCGAACAGCTACCTTCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCa   <  1:892725/150‑1 (MQ=37)
                                                                                                                                               aGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAg  <  2:947775/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAG  >  NC_000913/1708860‑1709152

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