Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,793,863 A→G C83R (TGT→CGT)  btuE ← thioredoxin/glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,8630AG97.4% 134.4 / NA 39C83R (TGT→CGT) btuEthioredoxin/glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/19);  minor base T (0/1);  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 = 1.00e+00

AATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAG  >  NC_000913/1793733‑1793996
                                                                                                                                  |                                                                                                                                     
aaTCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTGACTGAACACCGGGAACTTCACCCCCGATGTCGTGGTACGGTAAGtttt                                                                                                                              >  2:367426/1‑140 (MQ=255)
   ccGCTCTCTTCCGGCGCGACTGGGGCCGGCGCTGCGGTAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAAt                                                                                                                           <  1:363954/140‑1 (MQ=255)
   ccGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAAt                                                                                                                           <  1:1012002/140‑1 (MQ=255)
    cGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAAtc                                                                                                                          >  2:998055/1‑140 (MQ=255)
    cGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTACTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAAtc                                                                                                                          >  2:515004/1‑140 (MQ=255)
                   cgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGccc                                                                                                           <  2:1019429/140‑1 (MQ=255)
                     aCTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCgg                                                                                                         >  2:18902/1‑140 (MQ=255)
                     aCTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCgg                                                                                                         >  1:684164/1‑140 (MQ=255)
                       tGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGtt                                                                                                       >  2:1003971/1‑140 (MQ=255)
                                                 tttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACggg                                                                             <  1:18902/140‑1 (MQ=255)
                                                      tACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATcc                                                                        <  1:1003971/140‑1 (MQ=255)
                                                         aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAg                                                                     <  1:998055/140‑1 (MQ=255)
                                                         aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAg                                                                     >  1:938023/1‑140 (MQ=255)
                                                           cGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCa                                                                   >  2:926084/1‑140 (MQ=255)
                                                               tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCAt                                                               <  2:668519/140‑1 (MQ=255)
                                                                    tCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCGAGAAACTAGTTACATGGGAATACCAGCACCATAAAAc                                                          >  1:473101/1‑140 (MQ=255)
                                                                    tCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAc                                                          >  2:105122/1‑140 (MQ=255)
                                                                     ccTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcc                                                         >  1:134845/1‑140 (MQ=255)
                                                                     ccTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcc                                                         >  2:1012234/1‑140 (MQ=255)
                                                                         cGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                                     >  1:569182/1‑140 (MQ=255)
                                                                          gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                    <  2:569182/140‑1 (MQ=255)
                                                                               tAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAcc                                               <  1:515004/140‑1 (MQ=255)
                                                                                    tCAATCTTACTGAACTTCGGGACCGTCACCCCCCATGCGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGc                                          <  1:388871/140‑1 (MQ=255)
                                                                                    tCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGc                                          <  1:556146/140‑1 (MQ=255)
                                                                                      aaTCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCtt                                        >  2:827757/1‑140 (MQ=255)
                                                                                            aCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                  <  1:926084/140‑1 (MQ=255)
                                                                                            aCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                  <  2:684164/140‑1 (MQ=255)
                                                                                                 aCATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                             <  1:55238/140‑1 (MQ=255)
                                                                                                  cATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGCATCCCAGCCCCATAAACCCTCGATCAACCCAGGCTTTCTGAATATtct                            >  2:341951/1‑140 (MQ=255)
                                                                                                  cATCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACAATAAAACATCGTTCGACCCAGGCTTTCTGAATATtct                            >  1:317785/1‑140 (MQ=255)
                                                                                                    tCGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGGCCCAGGCTTTCTGAATATTCTcc                          <  2:938023/140‑1 (MQ=255)
                                                                                                     cGGGAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCAATCTACCCAGGCTTTCTGAATATTCTCCa                         >  1:424323/1‑140 (MQ=255)
                                                                                                     cGGGAACGTCACCCCCCATGTGGTGCTACTGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCa                         <  1:367426/140‑1 (MQ=255)
                                                                                                      gggAACGTCCCCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCaa                        <  1:583980/140‑1 (MQ=255)
                                                                                                      gggAACGTCACCCCCCATGTGGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCaa                        <  1:105122/140‑1 (MQ=255)
                                                                                                          aCGTCACCCCCCATGTTCTGGTACGGTAAGGTTCAATCTCTTCATCGCCGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                    <  2:473101/140‑1 (MQ=255)
                                                                                                                       tgtgGTGGTACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggcg       <  1:827757/140‑1 (MQ=255)
                                                                                                                           gtggtACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTGTCTGAATTTTCGCCAACTGCTGATATTGCGGCGTTaa   >  1:276504/1‑140 (MQ=255)
                                                                                                                            tggtACGGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg  >  2:252819/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                     
AATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAG  >  NC_000913/1793733‑1793996

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