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.G95.7% 74.2 / ‑2.3 23intergenic (‑514/+1137)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (8/14);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-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.
*ACB1435,386,5060CT95.6% 52.5 / ‑3.9 23intergenic (‑515/+1136)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/14);  minor base A (0/1);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-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.

TTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑CA‑GGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG  >  ACB143/5386415‑5386604
                                                                                              ||                                                                                                    
tttATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGccccc                                                             <  7:86355/135‑1 (MQ=255)
       aTAACAAATATTGATCGACTTGTTTCATTTAAGCCCATAATATGTCATCTCGGCCGCCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGAt                                                      <  7:265972/135‑1 (MQ=255)
         aaCAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCAc                                                   >  4:53616/1‑136 (MQ=255)
            aaaaaTTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCtt                                                >  6:132342/1‑136 (MQ=255)
                aTTGATCGAATTCTTTCATTTAAGACCCTAATATGCCATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGG‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                            <  6:113882/136‑1 (MQ=255)
                               tCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                               <  6:77952/134‑1 (MQ=255)
                                 aTTTCAGACCCTAATATCACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                           <  2:204322/136‑1 (MQ=255)
                                 aTTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                           <  3:65698/136‑1 (MQ=255)
                                 aTTTAAGACCAAAATATCACATCTCGGCCGGCTTGTCGACGACGGCGGGCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                           <  3:138308/136‑1 (MQ=255)
                                    tAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTg                         <  1:56774/135‑1 (MQ=255)
                                     aaGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCg                       >  6:274834/1‑136 (MQ=255)
                                       gACCCTATTATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGGGT‑AC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCggg                     <  7:188187/136‑1 (MQ=255)
                                        aCCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGt                    <  1:34559/136‑1 (MQ=255)
                                              atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  6:208046/136‑1 (MQ=255)
                                              atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCAACTTCCGGGGGCTTTATTATTGCGCTGCGGGTGCAAgg              >  2:57723/1‑136 (MQ=255)
                                              atatGACATCTCGGCCGGCTTGGCGACGCCGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCACGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  6:55363/136‑1 (MQ=255)
                                              atatGACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  7:185869/136‑1 (MQ=255)
                                                atGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGc            >  2:226982/1‑136 (MQ=255)
                                                 tGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGc            <  1:199562/135‑1 (MQ=255)
                                                      tctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGccc      >  8:96811/1‑136 (MQ=255)
                                                       ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑TGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCACCCGGAAGATCACCGTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCTTTCCCt     >  2:118239/1‑136 (MQ=255)
                                                          ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGccc      <  3:112646/132‑1 (MQ=255)
                                                          ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTGTACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg  >  4:191849/1‑136 (MQ=255)
                                                                                              ||                                                                                                    
TTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑CA‑GGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG  >  ACB143/5386415‑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: ‑