Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,386,506 1 bp→GT 100% intergenic (‑515/+1136) 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,5051.G86.4% 70.0 / 4.7 22intergenic (‑514/+1137)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (12/7);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 7.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*ACB1435,386,5060CT85.0% 48.0 / ‑0.0 22intergenic (‑515/+1136)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (11/7);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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: