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.96.8% 246.2 / 3.4 63coding (449/1662 nt)tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (24/37);  total (25/38)
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.42e-01

ATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATAT  >  NZ_CP009273/1966413‑1966584
                                                                                       |                                                                                    
atatCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaat                                                                                    >  1:896001/1‑89 (MQ=255)
atatCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTGGTCTGATCGaat                                                                                    >  1:736381/1‑89 (MQ=255)
 tatCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaata                                                                                   <  1:1002003/90‑3 (MQ=255)
  atCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaataa                                                                                  <  2:15994/90‑4 (MQ=255)
  atCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaataa                                                                                  <  2:457675/90‑4 (MQ=255)
  atCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaataa                                                                                  <  1:1213365/90‑4 (MQ=255)
   tCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaataag                                                                                 >  1:666584/1‑86 (MQ=255)
       gATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            <  2:355024/90‑1 (MQ=255)
       gATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            <  2:326023/90‑1 (MQ=255)
       gATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            <  1:182367/90‑1 (MQ=255)
       gATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            <  1:656080/90‑1 (MQ=255)
        aTACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCa                                                                           >  1:927503/1‑90 (MQ=255)
        aTACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCa                                                                           >  1:1080335/1‑90 (MQ=255)
            aGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTAt                                                                       <  2:927503/90‑1 (MQ=255)
             gTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTAtt                                                                      <  2:1572560/90‑1 (MQ=255)
                 ttCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCa                                                                  <  2:478450/90‑1 (MQ=255)
                 ttCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCa                                                                  <  2:1080335/90‑1 (MQ=255)
                 ttCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCa                                                                  <  1:534802/90‑1 (MQ=255)
                 ttCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTAATGCCa                                                                  <  1:1204216/90‑1 (MQ=255)
                     ctgctgAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCa                                                                  >  2:422362/1‑86 (MQ=255)
                     ctgctgAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCa                                                                  <  1:422362/86‑1 (MQ=255)
                      tgctgAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATc                                                             <  1:411284/90‑1 (MQ=255)
                      tgctgAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATc                                                             <  1:1010257/90‑1 (MQ=255)
                           aGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGAt                                                        <  2:1148826/90‑1 (MQ=255)
                                 tACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                                  >  1:1075962/1‑90 (MQ=255)
                                 tACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                                  >  2:1019163/1‑90 (MQ=255)
                                 tACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                                  >  1:1542385/1‑90 (MQ=255)
                                 tACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATaatcaa                                                  >  1:1311246/1‑90 (MQ=255)
                                    tGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            >  1:957559/1‑61 (MQ=255)
                                    tGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            <  2:957559/61‑1 (MQ=255)
                                    tGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAatcgatca                                               <  1:667355/90‑1 (MQ=255)
                                      aGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAg                                                          <  2:954318/77‑1 (MQ=255)
                                      aGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAg                                                          >  1:954318/1‑77 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            <  1:91792/90‑1 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            <  2:305486/90‑1 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            <  2:969685/90‑1 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            <  2:834223/90‑1 (MQ=255)
                                       gCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            <  1:575516/90‑1 (MQ=255)
                                          aaCGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            >  2:177172/1‑55 (MQ=255)
                                          aaCGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTcc                                                                            <  1:177172/55‑1 (MQ=255)
                                               ttCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAAc                                    >  1:1388500/1‑90 (MQ=255)
                                               ttCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAAc                                    >  2:1168151/1‑90 (MQ=255)
                                                 cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  1:1370929/90‑1 (MQ=255)
                                                 cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  1:1019163/90‑1 (MQ=255)
                                                 cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  2:896001/90‑1 (MQ=255)
                                                 cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  1:438075/90‑1 (MQ=255)
                                                 cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  1:2103/90‑1 (MQ=255)
                                                 cGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  2:1542385/90‑1 (MQ=255)
                                                  gCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  1:384146/90‑1 (MQ=255)
                                                   cccATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCtg                                <  1:1203685/90‑1 (MQ=255)
                                                      aTTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTgtgt                             >  2:631649/1‑90 (MQ=255)
                                                      aTTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTgtgt                             >  2:178642/1‑90 (MQ=255)
                                                        tgcattttgcattCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTaa                           <  1:1146155/90‑1 (MQ=255)
                                                              ttgcattCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGttt                     >  1:1249834/1‑90 (MQ=255)
                                                                   ttCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTtata                >  2:187147/1‑90 (MQ=255)
                                                                   ttCCCTGGGTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTtata                >  2:952368/1‑90 (MQ=255)
                                                                          ggTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGc                                  <  2:1256526/65‑1 (MQ=255)
                                                                          ggTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtca         >  1:874794/1‑90 (MQ=255)
                                                                          ggTTGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGGTCAGTTAACGc                                  >  1:1256526/1‑65 (MQ=255)
                                                                             tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca      <  2:1314100/90‑1 (MQ=255)
                                                                             tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca      <  1:926281/90‑1 (MQ=255)
                                                                             tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca      <  1:294838/90‑1 (MQ=255)
                                                                             tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca      <  1:1572244/90‑1 (MQ=255)
                                                                             tGGCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca      >  1:100886/1‑90 (MQ=255)
                                                                              ggCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatca      >  2:677143/1‑90 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            <  2:439380/51‑1 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAAtcagtt                                            >  1:439380/1‑51 (MQ=255)
                                                                              ggCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCaa     >  1:278715/1‑90 (MQ=255)
                                                                               gCTGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAat    >  1:691627/1‑90 (MQ=255)
                                                                                 tGAGCGAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAatat  >  1:758820/1‑90 (MQ=255)
                                                                                       |                                                                                    
ATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATAT  >  NZ_CP009273/1966413‑1966584

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