Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,212,513 A→G 7.0% I21I (ATT→ATC yidC ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,212,5130AG 7.0% 171.4 / 3.6 57I21I (ATT→ATCyidChypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GCGCACGTTTCACTTTTTCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTCACAACCAGCCAAAAAG  >  CP000730/2212374‑2212648
                                                                                                                                           |                                                                                                                                       
gcgcACGTTTCACTTTTTCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGataataa                                                                                                                                        <  2:554250/141‑1 (MQ=255)
 cgcACGCTTCACTTTTTCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATAAGAATTACCAATACAAGGagaataaa                                                                                                                                       >  2:526169/1‑140 (MQ=255)
   cACGTTTCACTTTTTCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAg                                                                                                                                     >  1:403118/1‑141 (MQ=255)
       tttCACTTTTACTTGAATTTTTTCTACATATGTCTTTGCAACTTTCAATTTTTGACGCATCATATGACTATTTTTATAGTTTGAGAACATGAATTGTAATAAAATAAAACGAATTACCAATACAAGTATAATAATAGCTaa                                                                                                                                 <  2:12074/141‑1 (MQ=255)
       gttCACTTTTTCTTGAATTTTTTCAACTTCTTGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTaa                                                                                                                                 <  1:84977/140‑1 (MQ=255)
             tttttCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATa                                                                                                                           <  1:651005/141‑1 (MQ=255)
                 tCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTg                                                                                                                       >  1:497170/1‑141 (MQ=255)
                 tCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTg                                                                                                                       <  1:434603/141‑1 (MQ=255)
                  cTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGt                                                                                                                      >  1:182846/1‑141 (MQ=255)
                  cTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGt                                                                                                                      <  1:174751/141‑1 (MQ=255)
                  cTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGt                                                                                                                      <  1:253032/141‑1 (MQ=255)
                      aaTTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATATCTAAACCATAATTGTCGtt                                                                                                                  >  2:469394/1‑141 (MQ=255)
                      aaTTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGtt                                                                                                                  >  1:343249/1‑141 (MQ=255)
                           tttCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTtaata                                                                                                             >  1:35414/1‑141 (MQ=255)
                              cAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGt                                                                                                          <  2:434549/141‑1 (MQ=255)
                               tatttttttcTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATTATACGTATTACTAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGtt                                                                                                         >  2:427017/10‑141 (MQ=255)
                               aaCTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGtt                                                                                                         >  2:276129/1‑141 (MQ=255)
                                    cTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTc                                                                                                    >  1:301331/1‑141 (MQ=255)
                                     tGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTcg                                                                                                   <  1:335280/141‑2 (MQ=255)
                                     tGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTcc                                                                                                   <  2:204370/141‑1 (MQ=255)
                                     tGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTcc                                                                                                   <  1:550405/141‑1 (MQ=255)
                                        cTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCccaa                                                                                                <  2:544642/141‑1 (MQ=255)
                                                cTTTCATTTTTTTACTCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAAccaatcca                                                                                        >  2:455855/1‑141 (MQ=255)
                                                cTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAAccaatcca                                                                                        >  2:455959/1‑141 (MQ=255)
                                                  ttCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAAccaatccaat                                                                                      <  1:188558/141‑1 (MQ=255)
                                                    cATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAACCCAATAc                                                                                    <  1:561282/141‑1 (MQ=255)
                                                      ttttttttttttttatatGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTTTCGTTTAATAAGTTATTTCCCAACCAATCCAATACAt                                                                                  >  2:132300/15‑141 (MQ=255)
                                                       tttttGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACAtt                                                                                 >  1:620662/1‑141 (MQ=255)
                                                              ttttttatatGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCATTACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCAtt                                                                          >  2:401219/7‑141 (MQ=255)
                                                              cGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTGAGAAAATAATACGAATTACCAATACAAGGATAATGATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCAtt                                                                          >  1:695854/1‑141 (MQ=255)
                                                              cGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAAccaatccaatgaaaaatgtatt                                                                          <  1:167807/141‑13 (MQ=255)
                                                              cGCATCATATGACTATTTTTATAGTTTGACAACATGAATGATTATTGAATAATACGAATTACCAATACAAGGATAATGATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCAtt                                                                          >  1:695829/1‑141 (MQ=255)
                                                                      atGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTAc                                                                  >  2:33100/1‑141 (MQ=255)
                                                                         aCTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGaa                                                               <  2:219972/141‑1 (MQ=255)
                                                                                 tataGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAAAGAATTga                                                       <  2:377522/141‑2 (MQ=255)
                                                                                    aGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGa                                                    >  1:55801/1‑141 (MQ=255)
                                                                                    aGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGa                                                    >  1:517706/1‑141 (MQ=255)
                                                                                            aacaTGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCa                                            >  2:621058/1‑141 (MQ=255)
                                                                                              caTGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACt                                          >  2:578085/1‑141 (MQ=255)
                                                                                                 gAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACg                                       <  1:517020/141‑1 (MQ=255)
                                                                                                  aaTTGTTATATAAGAATGCGAATTATCAATACAAGGATACTAATAGCTAAACCATAATTGTCTTTTAGTAAGTTTTTTCCCTACCAATCCAATTCTTTTTTCATTGGTTCTACTAATGTATTGTAGAAAAACCCACTACGt                                      <  2:316746/141‑1 (MQ=255)
                                                                                                   aTGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTTGATCTACGAATGTATTGTAGAAAAACCCACTACGtt                                     >  2:180285/1‑141 (MQ=255)
                                                                                                         ataaaataaTACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAg                               <  1:111785/141‑1 (MQ=255)
                                                                                                         ataaaataaTACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAg                               >  1:325351/1‑141 (MQ=255)
                                                                                                            aaataaTACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGtt                            <  1:585001/141‑1 (MQ=255)
                                                                                                                aataCGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGa                        <  1:436264/141‑1 (MQ=255)
                                                                                                                aataCGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGa                        >  1:299684/1‑141 (MQ=255)
                                                                                                                aataCGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGa                        >  2:519857/1‑141 (MQ=255)
                                                                                                                 ataCGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGaa                       >  1:568836/1‑141 (MQ=255)
                                                                                                                  taCGAATTACCAATACAAGGATAATGATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAAt                      <  2:695854/141‑1 (MQ=255)
                                                                                                                  taCGAATTACCAATACAAGGATAATGATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAAt                      <  2:695829/141‑1 (MQ=255)
                                                                                                                    cGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAg                    <  2:35414/141‑1 (MQ=255)
                                                                                                                        ttACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTcac                >  2:70476/1‑141 (MQ=255)
                                                                                                                           ccAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCACAAACAAACAAATACATTTTGCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGtaaaaaa             <  2:403118/141‑7 (MQ=255)
                                                                                                                            cAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCACAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAAAATAgaaaaaaaa            <  2:238093/141‑9 (MQ=255)
                                                                                                                                aCAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTCACAACCAGcc        >  1:303206/1‑141 (MQ=255)
                                                                                                                                 cAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTCACAACCAGCCa       <  1:526169/141‑1 (MQ=255)
                                                                                                                                  aaGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTATAAAAACCCACTACGTTTTTCAGGTTTAGAATAGtaaaaaaaaaaaaa      <  2:591830/141‑14 (MQ=255)
                                                                                                                                    ggATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTCACAACCAGCCaaaa    >  2:629104/1‑141 (MQ=255)
                                                                                                                                      ataataataGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTCACAACCAGCCAAAAAg  >  2:208887/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GCGCACGTTTCACTTTTTCTTGAATTTTTTCAACTTCTGGCTTTGCAACTTTCATTTTTTGACGCATCATATGACTATTTTTATAGTTTGACAACATGAATGGTAATAAAATAATACGAATTACCAATACAAGGATAATAATAGCTAAACCATAATTGTCGTTTAATAAGTTATTTCCCAACCAATCCAATACATTTTTCATTGGATCTACGAATGTATTGTAGAAAAACCCACTACGTTTTTCAGGTTTAGAATAGTCACAACCAGCCAAAAAG  >  CP000730/2212374‑2212648

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