Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,386,507:1 +C 100% intergenic (‑516/+1135) CJ019_05070 ← / ← nptII hypothetical protein/neomycin phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,386,5071.C86.4% 70.2 / 5.6 22intergenic (‑516/+1135)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (11/8);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑‑‑CAGGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCT  >  ACB143/5386410‑5386607
                                                                                                        |                                                                                                       
ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTGGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGc                                                                    <  5:121379/136‑1 (MQ=255)
         ttGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGaa                                                            <  4:33657/135‑1 (MQ=255)
                   aaaTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGGGT‑‑‑ACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTcc                                                 <  1:143191/136‑1 (MQ=255)
                   aaaTTGATCGAATTCTTTCATTTAACACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGT‑‑‑CCGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTcc                                                 <  7:362509/136‑1 (MQ=255)
                         aTCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                            <  7:212525/135‑1 (MQ=255)
                                      aTTAAAGCCCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGT‑‑‑CTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                              <  3:286513/136‑1 (MQ=255)
                                             aCCCTAATATGACATCTCGGCCGGCTTCGCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGAACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGt                       <  4:288234/136‑1 (MQ=255)
                                               ccTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACA‑‑AACCACGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTg                      >  8:31322/1‑135 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                 <  4:17800/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                 <  5:299482/136‑1 (MQ=255)
                                                    tatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAggg                >  7:249686/1‑136 (MQ=255)
                                                    tatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAggg                >  6:90161/1‑136 (MQ=255)
                                                     atGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATAACCTTACGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGc               >  3:210869/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTTCGCTGCGGGTGCCAGGGCGTGCCCt        >  6:82164/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt        >  1:75963/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt        >  7:171506/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt        >  4:124443/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt        >  3:266685/1‑136 (MQ=255)
                                                              cggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTg      >  7:363491/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGc                                                                    >  4:354843/1‑73 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGc                                                                    <  3:354843/73‑1 (MQ=255)
                                                                  cggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCATTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCt  >  2:145798/1‑136 (MQ=255)
                                                                                                        |                                                                                                       
CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑‑‑CAGGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCT  >  ACB143/5386410‑5386607

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