Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,324,838 T→A 17.1% Y31* (TAT→TAA) ‡ yciQ → enhancer of membrane protein expression; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,324,8380TA17.1% 163.4 / 22.8 58Y31* (TAT→TAA) ‡yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base A (6/4);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAG  >  NC_000913/1324703‑1324982
                                                                                                                                        |                                                                                                                                                
gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGa                                                                                                                                        >  2:117897/1‑146 (MQ=255)
gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGa                                                                                                                                        <  1:117897/146‑1 (MQ=255)
 aaGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaaga                                                                                                                                             >  1:930147/1‑140 (MQ=255)
 aaGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaaga                                                                                                                                             <  2:930147/140‑1 (MQ=255)
            ataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                          <  1:49384/148‑1 (MQ=255)
              aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                                  >  2:376292/1‑122 (MQ=255)
              aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                                  <  1:376292/122‑1 (MQ=255)
              aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaaga                                                                                                                                             <  2:267712/127‑1 (MQ=255)
              aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaaga                                                                                                                                             >  1:267712/1‑127 (MQ=255)
              aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                        >  1:907657/1‑148 (MQ=255)
                cAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAgg                                                                                                                                          >  2:863867/1‑128 (MQ=255)
                          tCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAgg                                                                                                                                          <  3:876060/118‑1 (MQ=255)
                          tCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAgg                                                                                                                                          >  4:876060/1‑118 (MQ=255)
                             aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTg                                                                                                                            <  1:910227/129‑1 (MQ=255)
                             aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTg                                                                                                                            >  2:910227/1‑129 (MQ=255)
                              cGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaaga                                                                                                                                             <  3:972495/111‑1 (MQ=255)
                              cGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaaga                                                                                                                                             >  4:972495/1‑111 (MQ=255)
                              cGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGaagaag                                                                                                                                           <  1:165962/113‑1 (MQ=255)
                              cGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGaagaag                                                                                                                                           >  2:165962/1‑113 (MQ=255)
                                    gAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                                  <  2:374070/100‑1 (MQ=255)
                                    gAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                                  >  1:374070/1‑100 (MQ=255)
                                       gACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaagaa                                                                                                                                            <  1:178990/103‑1 (MQ=255)
                                       gACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaagaa                                                                                                                                            >  2:178990/1‑103 (MQ=255)
                                                 gggAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAgg                                                                                                                                          <  1:935020/95‑1 (MQ=255)
                                                 gggAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCGCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAgg                                                                                                                                          >  2:935020/1‑95 (MQ=255)
                                                       tttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCAT‑AAGAAGGGa                                                                                                                                        <  2:631831/91‑1 (MQ=39)
                                                       tttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGa                                                                                                                                        >  1:631831/1‑91 (MQ=255)
                                                                            aTAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGaa                                                                                                                    >  4:492996/1‑90 (MQ=255)
                                                                            aTAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGaa                                                                                                                    <  3:492996/90‑1 (MQ=255)
                                                                                  gggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAgg                                                             >  4:420449/1‑139 (MQ=255)
                                                                                  gggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAgg                                                             <  3:420449/139‑1 (MQ=255)
                                                                                          tGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa                                                                              >  2:675447/1‑114 (MQ=255)
                                                                                          tGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa                                                                              <  1:675447/114‑1 (MQ=255)
                                                                                                     tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                >  2:335930/1‑85 (MQ=255)
                                                                                                     tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                <  1:335930/85‑1 (MQ=255)
                                                                                                                  aaaCACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTa                    >  3:121752/1‑148 (MQ=255)
                                                                                                                       cGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCt                                                                                          >  1:1101765/1‑73 (MQ=255)
                                                                                                                       cGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCt                                                                                          <  2:1101765/73‑1 (MQ=255)
                                                                                                                         gagaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATcag                                                  <  2:264783/111‑1 (MQ=255)
                                                                                                                         gagaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATcag                                                  >  1:264783/1‑111 (MQ=255)
                                                                                                                           gaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc                                                                                      >  1:750058/1‑73 (MQ=255)
                                                                                                                           gaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc                                                                                      <  2:750058/73‑1 (MQ=255)
                                                                                                                           gaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTa           <  2:907657/148‑1 (MQ=255)
                                                                                                                             tCCCTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg         <  3:182333/148‑1 (MQ=255)
                                                                                                                             tCCCTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg         >  4:182333/1‑148 (MQ=255)
                                                                                                                              cccTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAAt                                                     >  1:1084540/1‑103 (MQ=255)
                                                                                                                              cccTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAAt                                                     <  2:1084540/103‑1 (MQ=255)
                                                                                                                              cccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCaa                                                                                 <  3:550614/75‑1 (MQ=255)
                                                                                                                              cccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCaa                                                                                 >  4:550614/1‑75 (MQ=255)
                                                                                                                               ccTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCtttt                               >  3:399829/1‑124 (MQ=255)
                                                                                                                               ccTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCtttt                               <  4:399829/124‑1 (MQ=255)
                                                                                                                               ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATca                                                   >  3:612690/1‑104 (MQ=255)
                                                                                                                               ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATca                                                   <  4:612690/104‑1 (MQ=255)
                                                                                                                               ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTg                         >  1:86277/1‑130 (MQ=255)
                                                                                                                               ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTg                         <  2:86277/130‑1 (MQ=255)
                                                                                                                               ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg        <  2:547834/147‑1 (MQ=255)
                                                                                                                               ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg        >  1:547834/1‑147 (MQ=255)
                                                                                                                                 ttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                                >  4:1115738/1‑121 (MQ=255)
                                                                                                                                 ttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                                <  3:1115738/121‑1 (MQ=255)
                                                                                                                                   ctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCCGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTc                                                                >  4:375387/1‑87 (MQ=255)
                                                                                                                                   ctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCCGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTc                                                                <  3:375387/87‑1 (MQ=255)
                                                                                                                                    tt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAg  <  4:121752/148‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                
GAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAG  >  NC_000913/1324703‑1324982

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