Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,091 T→A 16.0% intergenic (+21/+28) pepT → / ← roxA peptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,0910TA16.0% 139.0 / 19.4 50intergenic (+21/+28)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base A (3/5);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGT  >  NC_000913/1186948‑1187236
                                                                                                                                               |                                                                                                                                                 
ccGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCTGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                              >  3:1785408/1‑149 (MQ=255)
ccGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                              <  4:882776/149‑1 (MQ=255)
     ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                                                                         <  3:1549998/149‑1 (MQ=255)
      cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCa                                                                                                                                          <  1:674504/147‑1 (MQ=255)
      cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCa                                                                                                                                          >  2:674504/1‑147 (MQ=255)
            aCTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                              >  4:19921/1‑137 (MQ=255)
            aCTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                              <  3:19921/137‑1 (MQ=255)
                      aCAACTATCATGGTAAGCATGCGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAAc                                                                                                                        <  4:1785408/149‑1 (MQ=255)
                        aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCaaa                                                                                                                         <  4:538871/146‑1 (MQ=255)
                        aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCaaa                                                                                                                         >  3:538871/1‑146 (MQ=255)
                        aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                                                                         <  1:279817/130‑1 (MQ=255)
                        aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                                                                         >  2:279817/1‑130 (MQ=255)
                             tCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                                                              <  1:1762035/136‑1 (MQ=255)
                             tCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                                                              >  2:1762035/1‑136 (MQ=255)
                                tGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCg                                                                                                              >  2:259114/1‑149 (MQ=255)
                                   tAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCa                                                                                                                                          >  4:1783197/1‑118 (MQ=255)
                                   tAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCa                                                                                                                                          <  3:1783197/118‑1 (MQ=255)
                                       cATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                                                             <  1:814473/111‑1 (MQ=255)
                                       cATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                                                             >  2:814473/1‑111 (MQ=255)
                                        aTGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                                                                        >  3:1866760/1‑115 (MQ=255)
                                        aTGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                                                                        <  4:1866760/115‑1 (MQ=255)
                                         tGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCaa                                                                                                     >  1:1823559/1‑149 (MQ=255)
                                            gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAc                                                                                                  >  2:659352/1‑149 (MQ=255)
                                              ttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAccc                                                                                                <  2:1823559/149‑1 (MQ=255)
                                                     tctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                              >  3:897123/1‑96 (MQ=255)
                                                     tctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                              <  4:897123/96‑1 (MQ=255)
                                                         gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGAt                                                                                                                                                <  4:36134/90‑1 (MQ=255)
                                                         gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGAt                                                                                                                                                >  3:36134/1‑90 (MQ=255)
                                                              tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCaa                                                                                                     >  3:1440586/1‑128 (MQ=255)
                                                              tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCaa                                                                                                     <  4:1440586/128‑1 (MQ=255)
                                                                          ggTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgcg                                                                    >  2:679889/1‑149 (MQ=255)
                                                                           gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGa                                                                   <  1:259114/149‑1 (MQ=255)
                                                                                       gTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGAt                                                       >  4:1636660/1‑149 (MQ=255)
                                                                                           gTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                                                                   >  4:1518082/1‑53 (MQ=255)
                                                                                           gTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                                                                   <  3:1518082/53‑1 (MQ=255)
                                                                                              ttGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAg                                                >  4:132078/1‑149 (MQ=255)
                                                                                                  cGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCa                                            >  2:1541794/1‑149 (MQ=255)
                                                                                                  cGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCa                                            <  1:679889/149‑1 (MQ=255)
                                                                                                  cGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCa                                            <  1:659352/149‑1 (MQ=255)
                                                                                                      ttAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGc                                        >  2:44425/1‑149 (MQ=255)
                                                                                                           ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaaaa                                   <  1:44425/149‑1 (MQ=255)
                                                                                                               cAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtct                               <  1:1541794/149‑1 (MQ=255)
                                                                                                                  cGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAtt                                                                                          <  3:693977/87‑1 (MQ=255)
                                                                                                                  cGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAtt                                                                                          >  4:693977/1‑87 (MQ=255)
                                                                                                                    gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCa                          <  3:132078/149‑1 (MQ=255)
                                                                                                                             gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCtt                                                    <  3:797767/114‑1 (MQ=255)
                                                                                                                             gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCtt                                                    >  4:797767/1‑114 (MQ=255)
                                                                                                                                     aTGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt         >  4:1070742/1‑149 (MQ=255)
                                                                                                                                     aTGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgc                                                                     <  1:1110306/89‑1 (MQ=39)
                                                                                                                                         ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTgg       >  3:188422/1‑147 (MQ=255)
                                                                                                                                         ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTgg       <  4:188422/147‑1 (MQ=255)
                                                                                                                                           cATGTGATGCCGCATCCGGCTTAAATCCAAACTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGt  >  2:688450/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGT  >  NC_000913/1186948‑1187236

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