Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,279 G→T 5.9% D396Y (GAC→TAC)  ygeH → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,2790GT 5.9% 242.8 / 5.4 68D396Y (GAC→TAC) ygeHputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base T (2/2);  total (34/34)
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.14e-01

ACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAG  >  NC_000913/2993149‑2993414
                                                                                                                                  |                                                                                                                                       
aCAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTc                                                                                                                        <  4:418882/148‑1 (MQ=255)
     cAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAg                                                                                                                                         >  2:359733/1‑126 (MQ=255)
     cAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAg                                                                                                                                         <  1:359733/126‑1 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 <  4:865931/116‑1 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 <  2:314258/116‑1 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 >  1:314258/1‑116 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 >  3:865931/1‑116 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCAc                                                                                                                                >  2:808645/1‑117 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCAc                                                                                                                                <  1:808645/117‑1 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACt                                                                                                                               >  2:906806/1‑118 (MQ=255)
                       tATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACt                                                                                                                               <  1:906806/118‑1 (MQ=255)
                          ttGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTCTATACAAATc                                                                                                                                  >  3:695867/1‑112 (MQ=255)
                          ttGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTCTATACAAATc                                                                                                                                  <  4:695867/112‑1 (MQ=255)
                             aaCAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACt                                                                                                                         >  4:816929/1‑118 (MQ=255)
                             aaCAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACt                                                                                                                         <  3:816929/118‑1 (MQ=255)
                             aaCAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGAt                                                                                           >  1:519014/1‑148 (MQ=255)
                              aCAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGa                                                                                                                                        <  4:163430/102‑1 (MQ=255)
                              aCAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGa                                                                                                                                        >  3:163430/1‑102 (MQ=255)
                                     aaGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 >  2:801830/1‑102 (MQ=255)
                                     aaGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 <  2:499590/102‑1 (MQ=255)
                                     aaGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 >  1:499590/1‑102 (MQ=255)
                                     aaGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCa                                                                                                                                 <  1:801830/102‑1 (MQ=255)
                                             cTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACa                                                                                                                            <  4:155322/99‑1 (MQ=255)
                                             cTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACa                                                                                                                            >  3:155322/1‑99 (MQ=255)
                                              tCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGata                                                                          <  2:519014/148‑1 (MQ=255)
                                                aaCTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCaga                                                                                                               >  2:1083793/1‑109 (MQ=255)
                                                aaCTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCaga                                                                                                               <  1:1083793/109‑1 (MQ=255)
                                                             tctctctACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATaac                                                                                        <  2:1079174/119‑2 (MQ=255)
                                                             tctctctACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATaaa                                                                                        >  1:1079174/1‑119 (MQ=255)
                                                                 tctACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGa                                                                                     <  4:1102552/118‑1 (MQ=255)
                                                                 tctACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGa                                                                                     >  3:1102552/1‑118 (MQ=255)
                                                                    actactaTAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGaa                                                                                                                     >  2:1003483/1‑83 (MQ=255)
                                                                    actactaTAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGaa                                                                                                                     <  1:1003483/83‑1 (MQ=255)
                                                                    actactaTAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCaa                                                    >  4:411816/1‑148 (MQ=255)
                                                                     ctactaTAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCaaa                                                   <  3:411816/148‑1 (MQ=255)
                                                                               gCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAAtgt                                                                                 <  3:57647/108‑1 (MQ=255)
                                                                               gCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAAtgt                                                                                 >  4:57647/1‑108 (MQ=255)
                                                                                 aCTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAg                                                                             <  2:897309/110‑1 (MQ=255)
                                                                                 aCTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAg                                                                             >  1:897309/1‑110 (MQ=255)
                                                                                  cTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGa                                                                                            >  2:950344/1‑94 (MQ=255)
                                                                                  cTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGa                                                                                            <  1:950344/94‑1 (MQ=255)
                                                                                       ccACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGc                                                                                                    >  4:433390/1‑81 (MQ=255)
                                                                                       ccACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGc                                                                                                    <  3:433390/81‑1 (MQ=255)
                                                                                       ccACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAAt                                                                                   <  1:648260/98‑1 (MQ=255)
                                                                                       ccACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAAt                                                                                   >  2:648260/1‑98 (MQ=255)
                                                                                                 aaTGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGata                                                                          >  2:1121758/1‑97 (MQ=255)
                                                                                                 aaTGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGata                                                                          <  1:1121758/97‑1 (MQ=255)
                                                                                                          aTTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTaa                                        >  2:569317/1‑122 (MQ=255)
                                                                                                          aTTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCGAACCCGCTCAAAAACAACATTaa                                        <  1:569317/122‑1 (MQ=255)
                                                                                                           ttGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTtgtg                                                                                             <  4:499004/68‑1 (MQ=255)
                                                                                                           ttGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTtgtg                                                                                             >  3:499004/1‑68 (MQ=255)
                                                                                                                 agGCGAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGaa                        <  3:613443/131‑1 (MQ=255)
                                                                                                                 agGCAAGGATTTGTATAGACAAATCACTACAACTCGa                                                                                                                      >  1:807479/1‑37 (MQ=255)
                                                                                                                 agGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGaa                        >  4:613443/1‑131 (MQ=255)
                                                                                                                 agGCAAGGAATTGTATAGACAAATCACTACAACTCGa                                                                                                                      <  2:807479/37‑1 (MQ=255)
                                                                                                                    cAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTcta                                  >  2:531155/1‑118 (MQ=255)
                                                                                                                    cAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTcta                                  <  1:531155/118‑1 (MQ=255)
                                                                                                                      aGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTAtgtg                                                                  >  1:616594/1‑84 (MQ=255)
                                                                                                                      aGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAg  >  4:717666/1‑148 (MQ=255)
                                                                                                                      aGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATCTAtgtg                                                                  <  2:616594/84‑1 (MQ=255)
                                                                                                                        gATTTGTCTATACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGa   <  3:630444/145‑1 (MQ=255)
                                                                                                                        gATTTGTCTATACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGa   >  4:630444/1‑145 (MQ=255)
                                                                                                                           ttGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAAcc                                                           <  2:1004734/86‑1 (MQ=255)
                                                                                                                           ttGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAAcc                                                           >  1:1004734/1‑86 (MQ=255)
                                                                                                                               ataGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGAtatgt                                                                       >  3:829747/1‑70 (MQ=255)
                                                                                                                               ataGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGAtatgt                                                                       <  4:829747/70‑1 (MQ=255)
                                                                                                                                 aGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAaaca                                               >  4:104733/1‑92 (MQ=255)
                                                                                                                                 aGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAaaca                                               <  3:104733/92‑1 (MQ=255)
                                                                                                                                 aGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTga                   >  2:705513/1‑120 (MQ=255)
                                                                                                                                 aGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTga                   <  1:705513/120‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
ACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATAGCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAG  >  NC_000913/2993149‑2993414

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: