Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,234,507 A→T 15.6% intergenic (‑52/+232) higB ← / ← rlmG mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,234,5070AT15.6% 181.5 / 24.6 64intergenic (‑52/+232)higB/rlmGmRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base T (5/5);  total (32/32)
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.94e-01

ATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGA  >  NC_000913/3234367‑3234648
                                                                                                                                            |                                                                                                                                              
aTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                        >  1:1276449/1‑148 (MQ=255)
aTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                        <  2:1276449/148‑1 (MQ=255)
 tCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACg                                                                                                                                      <  2:170455/149‑1 (MQ=255)
                  accaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                        >  3:1194384/1‑130 (MQ=255)
                  accaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                        <  4:1194384/130‑1 (MQ=255)
                   ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACgg                                                                                                                                     >  3:387203/1‑132 (MQ=255)
                   ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACgg                                                                                                                                     <  4:387203/132‑1 (MQ=255)
                          cGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaag                                                                                                             >  2:53942/1‑149 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑Tgg                                                                                                                                           >  3:1218148/1‑117 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑Tgg                                                                                                                                           <  4:1218148/117‑1 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgt                                                                                                                       >  2:1420540/1‑137 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgt                                                                                                                       <  1:1420540/137‑1 (MQ=255)
                                   ttGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                               <  4:402054/122‑1 (MQ=255)
                                   ttGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                               >  3:402054/1‑122 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAttt                                                                                                                             >  3:548937/1‑106 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAttt                                                                                                                             <  4:548937/106‑1 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                              >  2:133164/1‑121 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                              <  1:133164/121‑1 (MQ=255)
                                                      gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                 <  1:53942/149‑1 (MQ=255)
                                                                 tGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGa                                                                                                                                    <  1:546446/87‑1 (MQ=39)
                                                                 tGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCA‑ATCTGGCAACGGa                                                                                                                                    >  2:546446/1‑87 (MQ=38)
                                                                            ttATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                            <  1:528304/116‑1 (MQ=255)
                                                                            ttATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                            >  2:528304/1‑116 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAAc                                                                                                                                       >  1:268799/1‑65 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAAc                                                                                                                                       <  2:268799/65‑1 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATa                                                                                                         <  4:1182878/95‑1 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATa                                                                                                         >  3:1182878/1‑95 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                      >  4:482220/1‑98 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                      <  3:482220/98‑1 (MQ=255)
                                                                                             catcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTa                                                                          <  2:184460/117‑1 (MQ=255)
                                                                                             catcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTa                                                                          >  1:184460/1‑117 (MQ=255)
                                                                                               tcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGaa                                        >  3:624289/1‑149 (MQ=255)
                                                                                                   gTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTAt                                                                                                                           >  2:129050/1‑62 (MQ=255)
                                                                                                   gTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTAt                                                                                                                           <  1:129050/62‑1 (MQ=255)
                                                                                                       gTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                <  4:624289/149‑1 (MQ=255)
                                                                                                               aTTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACt                        >  4:369222/1‑149 (MQ=255)
                                                                                                                              aCATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAAt                                                                             <  4:1407014/81‑1 (MQ=255)
                                                                                                                              aCATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAAt                                                                             >  3:1407014/1‑81 (MQ=255)
                                                                                                                                aTCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                               >  2:632848/1‑109 (MQ=255)
                                                                                                                                aTCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                               <  1:632848/109‑1 (MQ=255)
                                                                                                                                 tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              <  1:1075672/109‑1 (MQ=255)
                                                                                                                                 tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              <  2:448789/109‑1 (MQ=255)
                                                                                                                                 tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              >  2:1075672/1‑109 (MQ=255)
                                                                                                                                 tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              >  1:448789/1‑109 (MQ=255)
                                                                                                                                 tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAt              >  2:621787/1‑141 (MQ=255)
                                                                                                                                 tCCGTTGCCA‑ATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAt              <  1:621787/141‑1 (MQ=255)
                                                                                                                                  ccGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                >  3:491726/1‑122 (MQ=255)
                                                                                                                                  ccGTTGCCA‑ATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                <  4:491726/122‑1 (MQ=255)
                                                                                                                                   cGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGCTTTATGGAAGCGCGa                                  >  4:501443/1‑119 (MQ=255)
                                                                                                                                   cGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGCTTTATGGAAGCGCGa                                  <  3:501443/119‑1 (MQ=255)
                                                                                                                                      tGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaaa                                                                                             <  2:461973/57‑1 (MQ=255)
                                                                                                                                      tGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaaa                                                                                             >  1:461973/1‑57 (MQ=255)
                                                                                                                                        ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATg                                                                   >  3:125218/1‑81 (MQ=255)
                                                                                                                                        ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATg                                                                   <  4:125218/81‑1 (MQ=255)
                                                                                                                                        ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            <  3:706670/104‑1 (MQ=255)
                                                                                                                                        ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            >  2:485863/1‑104 (MQ=255)
                                                                                                                                        ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            <  1:485863/104‑1 (MQ=255)
                                                                                                                                        ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            >  4:706670/1‑104 (MQ=255)
                                                                                                                                          aGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                 >  3:601051/1‑97 (MQ=255)
                                                                                                                                          aGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                 <  4:601051/97‑1 (MQ=255)
                                                                                                                                          aGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            >  4:337317/1‑102 (MQ=255)
                                                                                                                                          aGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            <  3:337317/102‑1 (MQ=255)
                                                                                                                                          aGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGa  <  4:1279411/144‑1 (MQ=255)
                                                                                                                                          aGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGa  >  3:1279411/1‑144 (MQ=255)
                                                                                                                                            |                                                                                                                                              
ATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGA  >  NC_000913/3234367‑3234648

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