Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 5,386,513 +AA intergenic (‑522/+1129) 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,5121.A88.9% 86.6 / 5.1 27intergenic (‑521/+1130)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (11/13);  total (12/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 = 1.00e+00
*ACB1435,386,5122.A88.9% 76.9 / 5.1 27intergenic (‑521/+1130)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (11/13);  total (12/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 = 8.62e-01

CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGTCGT‑‑‑CAGGATC‑‑A‑‑TCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG  >  ACB143/5386410‑5386604
                                                                                                             ||                                                                                               
ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGc                                                                 >  1:140318/1‑136 (MQ=255)
ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGc                                                                 >  8:196722/1‑136 (MQ=255)
         ttGATAACAAAAATTGATCTAATTCTTTAATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGGCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGaa                                                         <  1:381239/135‑1 (MQ=255)
               aCAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCAc                                                   >  4:93987/1‑135 (MQ=255)
                      ttGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                            <  4:38452/135‑1 (MQ=255)
                                    tCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                               <  8:284278/134‑1 (MQ=255)
                                    tCATTAAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGGT‑‑‑TACGCTAC‑‑AACCCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                               <  6:40934/134‑1 (MQ=255)
                                      aTTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCGCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                           <  7:90429/136‑1 (MQ=255)
                                      aTTTAACACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTAC‑‑AACCCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                           <  7:36010/136‑1 (MQ=255)
                                             aCCCTAATATGACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGt                    <  1:155497/136‑1 (MQ=255)
                                               ccTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGCGGGCTTTTTTATTGCGCTGCGGGTGc                  >  6:291024/1‑136 (MQ=255)
                                                   ttatGACATCGCGGCCGGCTGGTCGACGACGGGGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  2:383321/135‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  6:139809/136‑1 (MQ=255)
                                                   atatGACATCTCGGCCGGCGTGTCGACGACGGCGGTCGCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGAGGTTGCCAgg              <  8:255959/136‑1 (MQ=255)
                                                   atatCACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  1:143764/136‑1 (MQ=255)
                                                    tatGACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg              <  4:337739/135‑1 (MQ=255)
                                                            ctcGGCCGGCTTGTTGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCGTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt     >  5:3657/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt     >  5:69561/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt     >  6:362927/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt     >  1:216374/1‑136 (MQ=255)
                                                            ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTATTATTGCGCTGCGGGAGCCAGGGCGTGccc      >  3:212377/1‑135 (MQ=255)
                                                             tcGGCCGGCTTGTCGACCACGGCGGTCTCAGTGGTGTACGGTACA‑‑A‑‑TCCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGATGCGGGTGCCAGGGCGTGACCtt    >  5:313863/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg  >  5:132480/1‑136 (MQ=255)
                                                               ggccggcTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg  >  7:312245/1‑136 (MQ=255)
                                                                gccggcTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt    <  4:348427/133‑1 (MQ=255)
                                                                gccggcTTGTCGACGACGGCGGTCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt    <  4:192186/133‑1 (MQ=255)
                                                                gccggcTTGTCGACGACGGCGGGCTCAGTGGTG‑‑TACGGTACAAA‑‑CCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCtt    <  4:183160/133‑1 (MQ=255)
                                                                                                             ||                                                                                               
CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGTCGT‑‑‑CAGGATC‑‑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: