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.G94.7% 57.9 / ‑2.9 19intergenic (‑514/+1137)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (8/10);  total (8/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 = 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.
*ACB1435,386,5060CT94.6% 40.6 / ‑2.8 19intergenic (‑515/+1136)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/10);  minor base G (0/1);  total (8/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 = 6.23e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGTCGTCA‑‑GGATC‑‑A‑‑TCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG  >  ACB143/5386410‑5386604
                                                                                                 ||                                                                                                       
ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGc                                                                 <  1:15380/136‑1 (MQ=255)
    ttttATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGccccc                                                             >  5:76582/1‑136 (MQ=255)
            aTAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGAt                                                      <  6:184508/135‑1 (MQ=255)
                 aaaaaTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCtt                                                >  2:175088/1‑136 (MQ=255)
                   aaaTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTcc                                              <  6:235747/136‑1 (MQ=255)
                     aTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                            <  5:82931/136‑1 (MQ=255)
                         aTCGAATTCTTTCATTTAAGCCCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTAC‑‑AACCCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                         <  5:1259/135‑1 (MQ=255)
                                  tttCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                               <  6:7437/136‑1 (MQ=255)
                                    tCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGAACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                               <  7:180546/134‑1 (MQ=255)
                                            gACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCggg                     <  3:107647/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  6:92953/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  3:126192/136‑1 (MQ=255)
                                                    tatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAggg             >  4:216371/1‑136 (MQ=255)
                                                     atGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              >  7:108887/1‑134 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt     >  5:249295/1‑136 (MQ=255)
                                                              cggccggcTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTg   >  5:56563/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg  >  7:177812/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg  >  8:225118/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGGCTCAGTGGTGTACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt    <  1:132415/134‑1 (MQ=255)
                                                                                                 ||                                                                                                       
CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGTCGTCA‑‑GGATC‑‑A‑‑TCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG  >  ACB143/5386410‑5386604

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: