Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,950,940 T→C T508A (ACG→GCG)  recD ← exodeoxyribonuclease V subunit RecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,950,9400TC97.5% 138.7 / ‑4.7 40T508A (ACG→GCG) recDexodeoxyribonuclease V subunit RecD
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/19);  total (20/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

AAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTG  >  NC_000913/2950804‑2951072
                                                                                                                                        |                                                                                                                                    
aaaaTCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGcccc                                                                                                                                   <  2:1007695/140‑1 (MQ=255)
aaaaTCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGcccc                                                                                                                                   <  2:103224/140‑1 (MQ=255)
aaaaTCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGcccc                                                                                                                                   <  1:55758/140‑1 (MQ=255)
  aaTCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCt                                                                                                                                 <  1:1110706/140‑1 (MQ=255)
            gcATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATc                                                                                                                       <  1:485723/140‑1 (MQ=255)
            gcATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATc                                                                                                                       <  2:685918/140‑1 (MQ=255)
                  tCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAgcgc                                                                                                                 <  2:868953/140‑1 (MQ=255)
                      aCTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATa                                                                                                             <  1:1092779/140‑1 (MQ=255)
                            aTCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGACCCAGCGCAATACCGata                                                                                                       >  2:165786/1‑140 (MQ=255)
                             tCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGatat                                                                                                      <  1:1030985/140‑1 (MQ=255)
                             tCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGatat                                                                                                      <  1:516018/140‑1 (MQ=255)
                               ccTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGGACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCg                                                                                                    >  2:1014620/1‑140 (MQ=255)
                                cTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGc                                                                                                   <  2:658477/140‑1 (MQ=255)
                                  gCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCa                                                                                                 <  2:479107/140‑1 (MQ=255)
                                                gTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAg                                                                                   >  1:100153/1‑140 (MQ=255)
                                                gTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAg                                                                                   >  2:765643/1‑140 (MQ=255)
                                                gTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAg                                                                                   >  1:286135/1‑140 (MQ=255)
                                                  cATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAgcg                                                                                 <  1:552612/140‑1 (MQ=255)
                                                    tCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAgcgcg                                                                               <  1:765643/140‑1 (MQ=255)
                                                           cGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCAt                                                                        >  1:55682/1‑140 (MQ=255)
                                                           cGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCAt                                                                        >  2:1048590/1‑140 (MQ=255)
                                                           cGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCAt                                                                        >  1:679881/1‑140 (MQ=255)
                                                                      gCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAatc                                                             >  1:102286/1‑140 (MQ=255)
                                                                                   aCTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGAc                                                <  2:100153/140‑1 (MQ=255)
                                                                                    cTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGAcc                                               >  2:684626/1‑140 (MQ=255)
                                                                                     tCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCt                                              <  1:684626/140‑1 (MQ=255)
                                                                                     tCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCt                                              <  2:55682/140‑1 (MQ=255)
                                                                                     tCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCt                                              <  2:679881/140‑1 (MQ=255)
                                                                                         ttGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGt                                          <  1:77341/140‑1 (MQ=255)
                                                                                               aGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAAc                                    <  2:1103939/140‑1 (MQ=255)
                                                                                                 aCTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACga                                  >  1:621251/1‑140 (MQ=255)
                                                                                                           gCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                        >  1:909794/1‑140 (MQ=255)
                                                                                                           gCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                        >  1:772795/1‑140 (MQ=255)
                                                                                                           gCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                        >  1:580691/1‑140 (MQ=255)
                                                                                                           gCCGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                        >  1:379936/1‑140 (MQ=255)
                                                                                                            ccGTCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGTGGATg                       >  1:735362/1‑140 (MQ=255)
                                                                                                            ccGTCCGGCATCACAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGCTCTCGCCATTAAATAACCCCAGTGCGCTGTCATTACGGGCAATCATCACCGGACGACCTTCGTACTAACGAGAGTGCggcta                       >  2:790667/1‑137 (MQ=255)
                                                                                                               tCCGGCATCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACg                    <  2:580691/140‑1 (MQ=255)
                                                                                                                      tCGCAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAAtt             >  2:91401/1‑140 (MQ=255)
                                                                                                                         cAAACCAGACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTg          <  1:1014620/140‑1 (MQ=255)
                                                                                                                                gACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt   <  2:286135/140‑1 (MQ=255)
                                                                                                                                gACGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt   >  1:720386/1‑140 (MQ=255)
                                                                                                                                 aCGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTg  >  2:344865/1‑140 (MQ=255)
                                                                                                                                 aCGCGCGCCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTg  >  1:545068/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGTAGTTTCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTG  >  NC_000913/2950804‑2951072

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