Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,994,311 G→T 7.1% F132L (TTC→TTA uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,994,3110GT 7.1% 187.5 / 3.0 56F132L (TTC→TTAuvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base T (2/2);  total (31/25)
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.34e-01

GCACGGTCCCAGACAGCGCCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTT  >  NC_000913/1994173‑1994450
                                                                                                                                          |                                                                                                                                           
gCACGGTCCCAGACAGCGCCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGa                                                                                                                                            <  2:1158567/140‑1 (MQ=255)
               gcgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATAtt                                                                                                                             <  1:293651/140‑1 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGTAACGGGCCGAAATATTc                                                                                                                            >  2:834468/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGTAACGGGCCGAAATATTc                                                                                                                            >  2:834455/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGACCGGGCCGAAATATTc                                                                                                                            >  1:322076/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:1201464/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  2:957304/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  2:915644/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  2:915516/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  2:1070597/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  2:1047159/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:937474/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:786830/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:624983/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:514069/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:463252/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:455046/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:137951/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:322056/1‑140 (MQ=255)
                cgcCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTATTGCAGTAGCGCCAGTGTTTCACGTACTGCATAGCCATTCGGGAACGGGCCGAAATATTc                                                                                                                            >  1:455621/1‑140 (MQ=255)
                         agctATTGCAGACACGGTCGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGc                                                                                                                   <  2:1233876/137‑1 (MQ=255)
                          ggTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCa                                                                                                                  <  2:588504/140‑1 (MQ=255)
                          ggTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCa                                                                                                                  <  2:10953/140‑1 (MQ=255)
                          ggTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCa                                                                                                                  <  2:10819/140‑1 (MQ=255)
                            tATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATg                                                                                                                >  2:200661/1‑140 (MQ=255)
                                           cgcgAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATg                                                                                                 >  1:308772/1‑140 (MQ=255)
                                                          tAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCgg                                                                                  <  2:944086/140‑1 (MQ=255)
                                                             acacTTTGTTCGCACTGGCGAATGGGGAAAATCATTTGCAGTAGCGCCTGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCgggtg                                                                               <  2:1121909/140‑1 (MQ=255)
                                                             acacTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGGACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCCCGATGCATCGCCAGACGCgggtg                                                                               <  2:576513/140‑1 (MQ=255)
                                                                   ttttCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATc                                                                         >  2:708311/1‑140 (MQ=255)
                                                                     aacGCACTGGCGAATGGGGAAAATCTAATGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATcac                                                                       <  2:76890/138‑1 (MQ=255)
                                                                        gCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATctccac                                                                    >  2:1187122/1‑140 (MQ=255)
                                                                             ggCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAgg                                                               >  1:440957/1‑140 (MQ=255)
                                                                             ggCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAgg                                                               <  1:1231145/140‑1 (MQ=255)
                                                                             ggCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAgg                                                               <  1:461796/140‑1 (MQ=255)
                                                                                     gggAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATaaa                                                       <  2:822621/140‑1 (MQ=255)
                                                                                       gAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAg                                                     <  2:721260/140‑1 (MQ=255)
                                                                                                     gTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTatca                                       <  1:672379/140‑1 (MQ=255)
                                                                                                        gcgcCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATcg                                    <  1:20294/140‑1 (MQ=255)
                                                                                                        gcgcCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATcg                                    <  1:887917/140‑1 (MQ=255)
                                                                                                          gcCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTCATCATcgcg                                  <  2:624983/140‑1 (MQ=255)
                                                                                                                  ttCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCaaaa                          <  1:1108394/140‑1 (MQ=255)
                                                                                                                         aCGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTa                   >  2:33051/1‑140 (MQ=255)
                                                                                                                               tAGCCATTCGGTAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCgg             <  1:834468/140‑1 (MQ=255)
                                                                                                                               tAGCCATTCGGTAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCgg             <  1:834455/140‑1 (MQ=255)
                                                                                                                                aGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGc            <  2:170528/140‑1 (MQ=255)
                                                                                                                                aGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGc            <  2:170536/140‑1 (MQ=255)
                                                                                                                                    aTTCGGGAACTGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGTGTGGTTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGAt        <  2:985440/140‑1 (MQ=255)
                                                                                                                                    aTTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGTGTGGGTATCACCACTCATGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGAt        <  1:200661/140‑1 (MQ=255)
                                                                                                                                     ttCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCCGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCAGCGCGTAGCAAAACGGTGTCACGCGGCGGAtc       >  2:629003/1‑139 (MQ=255)
                                                                                                                                     ttCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATa       >  2:1170187/1‑140 (MQ=255)
                                                                                                                                       cGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga     >  1:213992/1‑140 (MQ=255)
                                                                                                                                       cGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga     >  1:247833/1‑140 (MQ=255)
                                                                                                                                       cGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga     >  1:628552/1‑140 (MQ=255)
                                                                                                                                       cGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga     >  1:510203/1‑140 (MQ=255)
                                                                                                                                       cGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga     >  1:501244/1‑140 (MQ=255)
                                                                                                                                       agggAACGGGCCGAAATATTCACCTTTGGCATGATTCGCACCACGATGCATCTCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGTTaga     >  1:90549/2‑140 (MQ=255)
                                                                                                                                          gAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGtt  <  1:63083/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GCACGGTCCCAGACAGCGCCCTATCTGGTATTGCAGACACGGACGCGAGCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTT  >  NC_000913/1994173‑1994450

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