Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913730,3890AT33.3% 20.6 / 13.7 18E269D (GAA→GATrhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base T (1/5);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTCCT  >  NC_000913/730291‑730404
                                                                                                  |               
gTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGaa                 <  1:1076647/99‑1 (MQ=33)
gTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGaa                 <  2:1085275/99‑1 (MQ=33)
gTGACGGATGGTGCCGGGCGTCACTGCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGaa                 <  2:2133933/99‑1 (MQ=21)
  gACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGATGGAAcc               <  1:2980798/99‑1 (MQ=21)
  gACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAAcc               >  1:2287921/1‑99 (MQ=34)
     ggATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTc            >  1:296413/1‑99 (MQ=255)
      gATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTcc           <  1:1102965/99‑1 (MQ=33)
        tGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGc         <  2:2485441/99‑1 (MQ=34)
         ggTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCTGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgc        <  1:1721413/99‑4 (MQ=11)
         ggTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgc        <  2:292445/99‑4 (MQ=11)
          gTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgcc       <  1:1552360/99‑5 (MQ=11)
          gTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgcc       <  2:1289217/99‑5 (MQ=11)
          gTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgcc       <  2:2069061/99‑5 (MQ=11)
          gTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCtt       >  1:1774689/1‑99 (MQ=34)
           tGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTTCGCttt      >  1:1514279/1‑99 (MQ=21)
           tGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCttt      >  2:1898109/1‑99 (MQ=34)
            gCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCtttt     >  1:1466911/1‑99 (MQ=33)
             ccGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGACGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTc    <  1:2829278/99‑1 (MQ=18)
             ccGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTc    <  1:311306/99‑1 (MQ=33)
             ccGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTc    <  2:782002/99‑1 (MQ=33)
               gggCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTCCt  <  2:559112/99‑1 (MQ=31)
                                                                                                  |               
GTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTCCT  >  NC_000913/730291‑730404

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