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.C95.2% 73.2 / ‑2.5 21intergenic (‑516/+1135)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (10/10);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGTCAGGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCTCC  >  ACB143/5386410‑5386609
                                                                                                     |                                                                                                        
ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGTTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGc                                                                      >  2:79784/1‑136 (MQ=255)
ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGc                                                                      >  8:380635/1‑136 (MQ=255)
     tttATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCg                                                                 >  7:225296/1‑136 (MQ=255)
            aTAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGAt                                                           <  5:383010/135‑1 (MQ=255)
                 aaaaaTTGATCGAATTCTTTCATTAAAGACCCTAATATGACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCtt                                                     <  3:311992/136‑1 (MQ=255)
                     aTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                 <  7:290729/136‑1 (MQ=255)
                                      aTTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                                <  4:82967/136‑1 (MQ=255)
                                      aTTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                                <  5:352848/136‑1 (MQ=255)
                                             aCCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGt                         <  7:195165/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                   <  5:381243/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCACGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                   <  3:267386/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                   <  8:309452/136‑1 (MQ=255)
                                                   atatGACATCGCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTCCGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                   <  3:36109/136‑1 (MQ=255)
                                                         cATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGc             >  4:178863/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt          >  5:12580/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt          >  8:318050/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt          >  7:317368/1‑136 (MQ=255)
                                                              cggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCCCCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTg        >  6:25926/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg       >  1:36307/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg       >  3:188093/1‑136 (MQ=255)
                                                                    gcTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGAACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCTcc  <  2:86595/136‑1 (MQ=255)
                                                                                                     |                                                                                                        
CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGTCAGGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCTCC  >  ACB143/5386410‑5386609

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: