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.100.0% 85.7 / NA 20coding (449/1662 nt)tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (6/14);  total (6/14)

ATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCA  >  NZ_CP009273/1966413‑1966580
                                                                                       |                                                                                
atatCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaat                                                                                >  2:167332/1‑89 (MQ=255)
atatCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaat                                                                                >  2:88577/1‑89 (MQ=255)
 tatCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaata                                                                               <  2:27188/90‑3 (MQ=255)
            aGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTAt                                                                   <  1:167332/90‑1 (MQ=255)
                   cACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATa                                                            <  2:232981/90‑1 (MQ=255)
                   cACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATa                                                            <  2:224502/90‑1 (MQ=255)
                      tgctgAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATc                                                         <  1:43282/90‑1 (MQ=255)
                         tgAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAg                                                      <  1:122997/90‑1 (MQ=255)
                                   cTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaatc                                            <  1:333796/90‑1 (MQ=255)
                                   cTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaatc                                            <  2:294746/90‑1 (MQ=255)
                                    tGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAatcaatca                                           >  1:242620/1‑90 (MQ=255)
                                      aGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagt                                         >  2:41791/1‑90 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                              <  2:183349/84‑1 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                              >  1:183349/1‑84 (MQ=255)
                                               ttCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGAt                                                    <  1:190002/70‑1 (MQ=255)
                                               ttCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGAt                                                    >  2:190002/1‑70 (MQ=255)
                                               ttCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAAc                                <  1:325946/90‑1 (MQ=255)
                                               ttCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAAc                                <  1:163536/90‑1 (MQ=255)
                                                                             tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca  <  1:153901/90‑1 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                        <  2:246034/51‑1 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                        >  1:74113/1‑51 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                        <  2:74113/51‑1 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                        >  1:246034/1‑51 (MQ=255)
                                                                                       |                                                                                
ATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCA  >  NZ_CP009273/1966413‑1966580

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