Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 480,334 C→T *125* (TAG→TAA)  tomB ← Hha toxicity attenuator, conjugation‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913480,3340CT97.7% 130.1 / ‑4.3 43*125* (TAG→TAA) tomBHha toxicity attenuator, conjugation‑related protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/23);  total (19/24)
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.96e-01

GATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAA  >  NC_000913/480202‑480459
                                                                                                                                    |                                                                                                                             
gATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACaa                                                                                                                         >  2:380264/1‑139 (MQ=255)
       cGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAAc                                                                                                                  >  2:54681/1‑139 (MQ=255)
         tatttattTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTGGTTAAAGGTTTTTCGGACATACTGCTACCTATGGCTGTAATAATTTTAACAAGATAAACTc                                                                                                                >  1:512517/1‑139 (MQ=255)
           tttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGc                                                                                                              >  2:418963/1‑139 (MQ=255)
            ttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCa                                                                                                             <  2:704710/139‑1 (MQ=255)
            ttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCa                                                                                                             <  1:421975/139‑1 (MQ=255)
                ttttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGAt                                                                                                         >  1:53223/1‑139 (MQ=255)
                                  ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGa                                                                                       <  1:418963/139‑1 (MQ=255)
                                      cGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACaaa                                                                                   >  1:67126/1‑139 (MQ=255)
                                      cGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGCTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACaaa                                                                                   >  1:143532/1‑139 (MQ=255)
                                         gtCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaaca                                                                                <  2:53223/139‑1 (MQ=255)
                                          tCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaacaa                                                                               <  2:222034/139‑1 (MQ=255)
                                          tCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaacaa                                                                               <  1:99989/139‑1 (MQ=255)
                                                gTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAAt                                                                         <  2:67126/139‑1 (MQ=255)
                                                                     aTTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGc                                                    <  2:143532/139‑1 (MQ=255)
                                                                      ttAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGcc                                                   >  2:101735/1‑139 (MQ=255)
                                                                       tAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCa                                                  <  2:387449/139‑1 (MQ=255)
                                                                            aatCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCt                                             <  1:380264/139‑1 (MQ=255)
                                                                             atCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTg                                            >  1:341576/1‑139 (MQ=255)
                                                                                ggTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAg                                         <  1:101735/139‑1 (MQ=255)
                                                                                ggTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATATTAACAAGATAAACTCGCAGGATACTCTTTCGTCGCATTGACAAAACAACTGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAg                                         <  2:512517/139‑1 (MQ=255)
                                                                                      cGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                   <  1:780070/139‑1 (MQ=255)
                                                                                      cGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                   <  1:684841/139‑1 (MQ=255)
                                                                                            aGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATAtt                             <  2:688103/139‑1 (MQ=255)
                                                                                                 tttCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCt                                             >  1:470308/1‑118 (MQ=255)
                                                                                                 tttCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCt                                             <  2:470308/118‑1 (MQ=255)
                                                                                                  ttCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATAcc                       >  2:706/1‑139 (MQ=255)
                                                                                                    cGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCAt                     <  1:706/139‑1 (MQ=255)
                                                                                                              tCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaa           <  2:341576/139‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGc                                                    >  1:142266/1‑91 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGc                                                    <  2:142266/91‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATTCTg                                   >  2:142852/1‑108 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                   <  1:142852/108‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGc                                  >  2:204553/1‑109 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGc                                  <  1:204553/109‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCAt                     <  2:171436/122‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCAt                     >  1:171436/1‑122 (MQ=255)
                                                                                                                      tGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTtc                                                                                                 <  2:445636/45‑1 (MQ=255)
                                                                                                                      tGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTtc                                                                                                 >  1:445636/1‑45 (MQ=255)
                                                                                                                      tGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATaa   >  2:8399/1‑139 (MQ=255)
                                                                                                                       ggTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATaaa  >  2:405773/1‑139 (MQ=255)
                                                                                                                             aataatTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAAt                          >  2:38661/1‑109 (MQ=255)
                                                                                                                             aataatTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAAt                          <  1:38661/109‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
GATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAA  >  NC_000913/480202‑480459

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