Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,966,502 Δ1 bp coding (447/1662 nt) tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,966,5000A.95.0% 70.8 / ‑0.8 20coding (449/1662 nt)tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (7/12);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAA  >  NZ_CP009273/1966418‑1966581
                                                                                  |                                                                                 
gcgATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCt                                                                           >  2:206555/1‑90 (MQ=255)
  gATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                         <  2:436932/90‑1 (MQ=255)
       aGTTTTTCTCTGCTGAGGGCGTACTGAGCAAAGGCTTCGCCCCTTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTAt                                                                    <  1:304723/90‑1 (MQ=255)
        gTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTCCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTAt                                                                    <  1:3766/90‑1 (MQ=255)
           tttCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCCCCAGTATTGCCa                                                               <  1:169549/90‑1 (MQ=255)
            ttCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCa                                                               <  1:111780/90‑1 (MQ=255)
             tCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCAt                                                              >  2:152966/1‑90 (MQ=255)
                     gAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGa                                                      >  2:327307/1‑90 (MQ=255)
                            tACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                               >  1:276383/1‑90 (MQ=255)
                              cTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaatc                                             <  2:266306/90‑1 (MQ=255)
                                 aGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagt                                          >  1:24026/1‑90 (MQ=255)
                                  gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                         <  1:180290/90‑1 (MQ=255)
                                     aaCGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaatc                                             <  1:330652/83‑1 (MQ=255)
                                     aaCGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaatc                                             >  2:330652/1‑83 (MQ=255)
                                          ttCGCCCATTGCATTTTGCATCCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAAc                                 <  1:71830/90‑1 (MQ=255)
                                            cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                               <  2:38485/90‑1 (MQ=255)
                                            cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                               <  2:74835/90‑1 (MQ=255)
                                                                       ttGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatc    <  2:221750/90‑1 (MQ=255)
                                                                        tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca   <  2:201310/90‑1 (MQ=255)
                                                                         ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCaa  >  1:212973/1‑90 (MQ=255)
                                                                                  |                                                                                 
GCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAA  >  NZ_CP009273/1966418‑1966581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: