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.0% 98.8 / ‑6.8 33*125* (TAG→TAAtomBHha toxicity attenuator, conjugation‑related protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (15/17);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTG  >  NC_000913/480202‑480463
                                                                                                                                    |                                                                                                                                 
gATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACaa                                                                                                                             >  2:341758/1‑139 (MQ=255)
    aTTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATa                                                                                                                         <  1:107295/139‑1 (MQ=255)
                                   cAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGAc                                                                                          >  2:92662/1‑139 (MQ=255)
                                         gtCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaaca                                                                                    >  1:256777/1‑139 (MQ=255)
                                                                      ttAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGc                                                                                                <  1:309272/98‑1 (MQ=255)
                                                                      ttAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGc                                                                                                >  2:309272/1‑98 (MQ=255)
                                                                              tCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCa                                                                                                                 <  2:663218/73‑1 (MQ=255)
                                                                              tCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCa                                                                                                                 >  1:663218/1‑73 (MQ=255)
                                                                                       gTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGc                                      >  2:839842/1‑139 (MQ=255)
                                                                                                  ttCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATAcc                           <  2:415872/139‑1 (MQ=255)
                                                                                                         aTACTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAacta                    <  1:341758/139‑1 (MQ=255)
                                                                                                            cTTCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTg                 <  2:376926/139‑1 (MQ=255)
                                                                                                              tCTACCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaa               <  2:748889/139‑1 (MQ=255)
                                                                                                                 aCCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCa                                                      >  1:676207/1‑97 (MQ=255)
                                                                                                                 aCCTATGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCa                                                      <  2:676207/97‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCAtt                                                    >  1:824378/1‑95 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCAtt                                                    <  2:824378/95‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGa                                               >  1:345411/1‑100 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGa                                               <  2:345411/100‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                       >  1:70190/1‑108 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                       <  2:70190/108‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATAtt                                 <  1:192540/114‑1 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATAtt                                 >  2:192540/1‑114 (MQ=255)
                                                                                                                     aTGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATa        >  1:294588/1‑139 (MQ=255)
                                                                                                                      tGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTtcgtcg                                                                                                 <  2:585325/49‑1 (MQ=255)
                                                                                                                      tGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTtcgtcg                                                                                                 >  1:585325/1‑49 (MQ=255)
                                                                                                                      tGGTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATaa       <  2:294588/139‑1 (MQ=255)
                                                                                                                       ggTTGTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGAATTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATaaa      >  1:580541/1‑139 (MQ=255)
                                                                                                                          tttaattaTTCTAAAAAGATAAACTAGCATTATTCTCTTACGGCGAATTGAAAAACGAACGGAATAGTAGACTTCGCGCCTTCCGCCATTTCTCAAGATCCTCCATATTCATACCATAACTACTGAACAACATAAAGgt   >  2:648284/3‑139 (MQ=255)
                                                                                                                           gTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGgtg  <  2:519805/139‑1 (MQ=255)
                                                                                                                           gTAATAATTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGgtg  <  2:580541/139‑1 (MQ=255)
                                                                                                                            taataaTTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                       >  2:205919/1‑101 (MQ=255)
                                                                                                                            taataaTTTTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                       <  1:205919/101‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
GATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTG  >  NC_000913/480202‑480463

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