Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 5,386,506 1 bp→GT 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.G90.5% 61.5 / 0.3 22intergenic (‑514/+1137)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (13/6);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*ACB1435,386,5060CT85.7% 47.4 / 0.6 22intergenic (‑515/+1136)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (12/6);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑‑‑CA‑‑GGATC‑‑AT‑‑CCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCT  >  ACB143/5386410‑5386607
                                                                                                     ||                                                                                                            
gcctttttATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGcc                                                                   <  5:75039/134‑1 (MQ=255)
   tttttATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGccccc                                                                <  2:21790/137‑1 (MQ=255)
             aaaCAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCAc                                                      <  2:8392/136‑1 (MQ=255)
                           cGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttt                                        >  4:361102/1‑137 (MQ=255)
                            gAATACTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTGGGCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttt                                       <  7:245488/137‑1 (MQ=255)
                                 ttttCATTTAACACCCTAATATCACATCTCGCCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGT‑‑‑GAACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                  <  1:170691/136‑1 (MQ=255)
                                      aTTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCt                             <  1:325080/137‑1 (MQ=255)
                                       tttAATACCGTAATATGACATCTCGGGAGGCTTGTCGACGACGGCGGTCTCAGTGGTGT‑‑‑CA‑‑GGTAC‑‑ATACCCCGACGCTAGCAAAGCATGAGAAAGCCACCGGAATATCTCGTTCCGGGGGCTTTTTTATTGCGCTg                            <  8:4602/137‑1 (MQ=255)
                                                  aaTATGACATCTCGGCCGGCTTTGCGACGACGGGGGTCTCAGTGGGGTA‑‑CT‑‑GTACAA‑‑AC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg                 <  5:360596/137‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAggg                >  8:161605/1‑137 (MQ=255)
                                                    tatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGc               >  4:296994/1‑137 (MQ=255)
                                                    tatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGc               >  6:388511/1‑137 (MQ=255)
                                                     atGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCg              >  5:155029/1‑137 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTTCGATGCGGGTGCCAGGGCGTGCCCtt       >  6:156043/1‑137 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt       >  8:189773/1‑137 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGAACACCTTCCGGGGGCGTTTTTATTTGGCAGCGGGTGCCAGGGCGGGCCCtt       >  1:194952/1‑137 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTAAGT‑‑GGTGTACG‑‑GTACA‑‑AA‑‑CACCTACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt       >  6:374407/1‑137 (MQ=255)
                                                            ctcGGCCGGCTTGGCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt       <  5:355881/137‑1 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCCCCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTggg    >  5:148875/1‑137 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTggg    >  1:305794/1‑137 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTggg    >  1:117363/1‑137 (MQ=255)
                                                                 ccggcTTGTCGACGACGGCGGTCTCAGT‑‑GGTG‑‑TACGGTACAAAC‑‑CCCGACGCTAGCAACGCATGAGAAAGCCCCAGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCt  >  1:243328/1‑137 (MQ=255)
                                                                                                     ||                                                                                                            
CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CGT‑‑‑CA‑‑GGATC‑‑AT‑‑CCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGGGCT  >  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: