Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,064,594 T→A 11.1% Q135L (CAG→CTG)  PP_1841 ← cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,064,5940TA11.1% 92.2 / 3.8 36Q135L (CAG→CTG) PP_1841cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base T (18/14);  new base A (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCA  >  NC_002947/2064454‑2064740
                                                                                                                                            |                                                                                                                                                  
gggCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCAGGCCAgcggca                                                                                                                                              >  7:185481/1‑146 (MQ=255)
gggCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCAGGCCAgcggca                                                                                                                                              <  8:185481/147‑2 (MQ=255)
 ggCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAg                                                                                                                                           >  6:269196/1‑149 (MQ=255)
     aaGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCca                                                                                                                                       <  2:102190/149‑1 (MQ=255)
               gcACCGGGCTGTCGGCCGGCAGATTCGCGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGCGCCAGCTTGCGAAGGGTGGCAATCTGCAATGGCCACTGTGCAGGGCGGATCCCGCTGCCCTGCGGCGCAGTCCCCACCAgg                                                                                                                                >  2:211039/1‑146 (MQ=255)
                cACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAg                                                                                                                             >  8:233436/1‑148 (MQ=255)
                 aCCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGaa                                                                                                                           <  7:233436/149‑1 (MQ=255)
                      gCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAg                                                                                                                             >  6:109402/1‑142 (MQ=255)
                      gCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAg                                                                                                                             <  5:109402/142‑1 (MQ=255)
                          tCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCg                                                                                                                    <  5:71165/147‑1 (MQ=255)
                           cggccggcAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCa                                                                                                                 >  6:271752/1‑149 (MQ=255)
                            ggccggcAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCAGGCCAgcggca                                                                                                                                              <  3:193306/119‑2 (MQ=255)
                            ggccggcAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCAGGCCAgcggca                                                                                                                                              >  4:193306/1‑118 (MQ=255)
                               cggcAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGCACGGCCCGGCCCtgg                                                                                                             >  4:159683/1‑147 (MQ=255)
                                          ggtcgggCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGa                                                                                                  >  1:246117/1‑149 (MQ=255)
                                               ggCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAg                                                                                             >  4:48301/1‑149 (MQ=255)
                                                  agcagcGCGGGAAACCGCTGGTCGACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTgg                                                                                          >  7:90318/1‑149 (MQ=255)
                                                      gcgcAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTggcggc                                                                                      <  1:50999/149‑1 (MQ=255)
                                                            gAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGc                                                                                >  2:21500/1‑149 (MQ=255)
                                                              aaCCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCtt                                                                              >  4:212458/1‑149 (MQ=255)
                                                                            cGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCa                                                                >  7:159063/1‑149 (MQ=255)
                                                                                     cTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAAcc                                                       <  5:271752/149‑1 (MQ=255)
                                                                                        gCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAg                                                     <  3:48301/148‑1 (MQ=255)
                                                                                          gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGcc                                                  >  5:81039/1‑149 (MQ=255)
                                                                                               gTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGGTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCAGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTg                                              <  1:211039/148‑1 (MQ=255)
                                                                                                   cAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGc                                         <  6:81039/149‑1 (MQ=255)
                                                                                                         ggAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGc                                   >  6:184743/1‑149 (MQ=255)
                                                                                                          gAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCt                                      <  8:159063/145‑1 (MQ=255)
                                                                                                          gAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTgg                                    <  2:246117/147‑1 (MQ=255)
                                                                                                             tGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGCCCAACCAGGCCCGTGCTGGCCCTGGCTgca                               <  1:21500/149‑1 (MQ=255)
                                                                                                                  aCTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAAt                          <  8:90318/149‑1 (MQ=255)
                                                                                                                  aCTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCACGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAAt                          >  3:99059/1‑149 (MQ=255)
                                                                                                                         aGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCACTGGCTGCAGCAATGCCGCAg                   >  3:56636/1‑149 (MQ=255)
                                                                                                                                  tggtgtGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCAttt          >  6:60330/6‑149 (MQ=255)
                                                                                                                                      cTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggc      <  5:269196/149‑1 (MQ=255)
                                                                                                                                           cTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCa  >  8:129110/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCA  >  NC_002947/2064454‑2064740

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