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

ATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGA  >  NC_000913/3234367‑3234648
                                                                                                                                            |                                                                                                                                              
aTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCaa                                                                                                                                        >  1:1276449/1‑148 (MQ=255)
aTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCaa                                                                                                                                        <  2:1276449/148‑1 (MQ=255)
 tCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACg                                                                                                                                      <  2:170455/149‑1 (MQ=255)
                  accaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCaa                                                                                                                                        >  3:1194384/1‑130 (MQ=255)
                  accaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCaa                                                                                                                                        <  4:1194384/130‑1 (MQ=255)
                   ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACgg                                                                                                                                     >  3:387203/1‑132 (MQ=255)
                   ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACgg                                                                                                                                     <  4:387203/132‑1 (MQ=255)
                          cGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaag                                                                                                             >  2:53942/1‑149 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTgg                                                                                                                                           >  3:1218148/1‑117 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTgg                                                                                                                                           <  4:1218148/117‑1 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCtgt                                                                                                                       >  2:1420540/1‑137 (MQ=255)
                            ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCtgt                                                                                                                       <  1:1420540/137‑1 (MQ=255)
                                   ttGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAt                                                                                                                               <  4:402054/122‑1 (MQ=255)
                                   ttGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAt                                                                                                                               >  3:402054/1‑122 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAttt                                                                                                                             >  3:548937/1‑106 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAAttt                                                                                                                             <  4:548937/106‑1 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                              >  2:133164/1‑121 (MQ=255)
                                                     agcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                              <  1:133164/121‑1 (MQ=255)
                                                      gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                 <  1:53942/149‑1 (MQ=255)
                                                                 tGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACGGa                                                                                                                                    <  1:546446/87‑1 (MQ=39)
                                                                 tGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACGGa                                                                                                                                    >  2:546446/1‑87 (MQ=38)
                                                                            ttATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                            <  1:528304/116‑1 (MQ=255)
                                                                            ttATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                            >  2:528304/1‑116 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAAc                                                                                                                                       >  1:268799/1‑65 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAAc                                                                                                                                       <  2:268799/65‑1 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATa                                                                                                         <  4:1182878/95‑1 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATa                                                                                                         >  3:1182878/1‑95 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                      >  4:482220/1‑98 (MQ=255)
                                                                                    tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                      <  3:482220/98‑1 (MQ=255)
                                                                                             catcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTa                                                                          <  2:184460/117‑1 (MQ=255)
                                                                                             catcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTa                                                                          >  1:184460/1‑117 (MQ=255)
                                                                                               tcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGaa                                        >  3:624289/1‑149 (MQ=255)
                                                                                                   gTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTAt                                                                                                                           >  2:129050/1‑62 (MQ=255)
                                                                                                   gTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTAt                                                                                                                           <  1:129050/62‑1 (MQ=255)
                                                                                                       gTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                <  4:624289/149‑1 (MQ=255)
                                                                                                               aTTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACt                        >  4:369222/1‑149 (MQ=255)
                                                                                                                              aCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAAt                                                                             <  4:1407014/81‑1 (MQ=255)
                                                                                                                              aCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAAt                                                                             >  3:1407014/1‑81 (MQ=255)
                                                                                                                                aTCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                               >  2:632848/1‑109 (MQ=255)
                                                                                                                                aTCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                               <  1:632848/109‑1 (MQ=255)
                                                                                                                                 tCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              <  1:1075672/109‑1 (MQ=255)
                                                                                                                                 tCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              <  2:448789/109‑1 (MQ=255)
                                                                                                                                 tCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              >  2:1075672/1‑109 (MQ=255)
                                                                                                                                 tCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtttt                                              >  1:448789/1‑109 (MQ=255)
                                                                                                                                 tCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAt              >  2:621787/1‑141 (MQ=255)
                                                                                                                                 tCCGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAt              <  1:621787/141‑1 (MQ=255)
                                                                                                                                  ccGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                >  3:491726/1‑122 (MQ=255)
                                                                                                                                  ccGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGAtt                                <  4:491726/122‑1 (MQ=255)
                                                                                                                                   cGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGCTTTATGGAAGCGCGa                                  >  4:501443/1‑119 (MQ=255)
                                                                                                                                   cGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGCTTTATGGAAGCGCGa                                  <  3:501443/119‑1 (MQ=255)
                                                                                                                                      tGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaaa                                                                                             <  2:461973/57‑1 (MQ=255)
                                                                                                                                      tGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaaa                                                                                             >  1:461973/1‑57 (MQ=255)
                                                                                                                                        ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATg                                                                   >  3:125218/1‑81 (MQ=255)
                                                                                                                                        ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATg                                                                   <  4:125218/81‑1 (MQ=255)
                                                                                                                                        ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            <  3:706670/104‑1 (MQ=255)
                                                                                                                                        ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            >  2:485863/1‑104 (MQ=255)
                                                                                                                                        ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            <  1:485863/104‑1 (MQ=255)
                                                                                                                                        ccAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            >  4:706670/1‑104 (MQ=255)
                                                                                                                                          aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                 >  3:601051/1‑97 (MQ=255)
                                                                                                                                          aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                 <  4:601051/97‑1 (MQ=255)
                                                                                                                                          aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            >  4:337317/1‑102 (MQ=255)
                                                                                                                                          aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                            <  3:337317/102‑1 (MQ=255)
                                                                                                                                          aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGa  <  4:1279411/144‑1 (MQ=255)
                                                                                                                                          aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGa  >  3:1279411/1‑144 (MQ=255)
                                                                                                                                            |                                                                                                                                              
ATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGA  >  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: