Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,234,7150GT20.0% 64.8 / 10.4 30intergenic (‑260/+24)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 G (11/13);  new base T (3/3);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTGGC  >  NC_000913/3234580‑3234858
                                                                                                                                       |                                                                                                                                               
gATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGt                                                                                                                                     >  5:147421/1‑148 (MQ=255)
                ttttGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc                                                                                                                     <  1:167781/148‑1 (MQ=255)
                ttttGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc                                                                                                                     >  2:39693/1‑148 (MQ=255)
                   tGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACcgac                                                                                                                  <  1:39693/148‑1 (MQ=255)
                     ttttATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTAAACGCCTTATCCGGc                                                                                                                                         >  1:137569/1‑123 (MQ=255)
                     ttttATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTAAACGCCTTATCCGGc                                                                                                                                         <  2:137569/123‑1 (MQ=255)
                          tGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                           >  6:103994/1‑116 (MQ=38)
                          tGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                           <  5:103994/116‑1 (MQ=38)
                                                       aTTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCa                                                                                          <  8:29412/136‑1 (MQ=255)
                                                       aTTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCa                                                                                          >  7:29412/1‑136 (MQ=255)
                                                            cGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                               <  1:136781/78‑1 (MQ=255)
                                                            cGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                               >  2:136781/1‑78 (MQ=255)
                                                                                      ccGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc                                                                                                                     >  8:84410/1‑78 (MQ=255)
                                                                                      ccGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc                                                                                                                     <  7:84410/78‑1 (MQ=255)
                                                                                                        cTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCCTTCAGTACCACAAATTTATTATTCGTCGCAATAgtg                                                          <  2:67/119‑1 (MQ=255)
                                                                                                        cTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCCTTCAGTACCACAAATTTATTATTCGTCGCAATAgtg                                                          >  1:67/1‑119 (MQ=255)
                                                                                                         tGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTtattat                                                                        >  5:7131/1‑104 (MQ=255)
                                                                                                         tGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTtattat                                                                        <  6:7131/104‑1 (MQ=255)
                                                                                                                aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAATTTTATTATtcgt                                                                    <  2:53452/101‑1 (MQ=255)
                                                                                                                aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAATTTTATTATtcgt                                                                    >  1:53452/1‑101 (MQ=255)
                                                                                                                  gCCGGATGCGGCGTAAACGCCTTATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGc                                                      >  8:22937/1‑113 (MQ=38)
                                                                                                                  gCCGGATGCGGCGTAAACGCCTTATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGc                                                      <  7:22937/113‑1 (MQ=38)
                                                                                                                     ggATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCa               <  3:126179/149‑1 (MQ=255)
                                                                                                                       aTAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGt                                                                                        >  8:132807/1‑74 (MQ=255)
                                                                                                                       aTAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGt                                                                                        <  7:132807/74‑1 (MQ=255)
                                                                                                                              gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTa                    >  1:126191/1‑135 (MQ=255)
                                                                                                                              gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTa                    <  2:126191/135‑1 (MQ=255)
                                                                                                                               tttACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACgg       >  4:111832/1‑147 (MQ=255)
                                                                                                                               tttACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACgg       <  3:111832/147‑1 (MQ=255)
                                                                                                                                 tACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTg    <  6:147421/148‑1 (MQ=255)
                                                                                                                                   cgccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTGGc  >  1:90293/1‑148 (MQ=255)
                                                                                                                                    gccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATcttctt                                   >  4:11382/1‑114 (MQ=255)
                                                                                                                                    gccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATcttctt                                   <  3:11382/114‑1 (MQ=255)
                                                                                                                                     ccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGaa                       >  2:137796/1‑125 (MQ=255)
                                                                                                                                     ccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGaa                       <  1:137796/125‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                               
GATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTGGC  >  NC_000913/3234580‑3234858

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