Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,718,230 T→G 6.1% V291V (GTT→GTG PP_5020 → methyl‑accepting chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,718,2300TG 6.1% 188.4 / 4.2 66V291V (GTT→GTGPP_5020methyl‑accepting chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base G (2/2);  total (30/36)
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.22e-01

CGGCAATGCCCAGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTACAACGCCGCGAC  >  NC_002947/5718083‑5718377
                                                                                                                                                    |                                                                                                                                                   
cGGCAATGCCCAGCTCGCCAGCGAGCAAGAGACAGCACCCATGGTTGAACAAAGCATGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGc                                                                                                                                                   <  2:690786/150‑1 (MQ=255)
    aaTGCCCAGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCAcc                                                                                                                                               >  1:91131/1‑150 (MQ=255)
    aaTGCCCAGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCAcc                                                                                                                                               >  2:517246/1‑150 (MQ=255)
     aTGCCCAGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGACAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCg                                                                                                                                              >  2:808228/1‑150 (MQ=255)
        cccAGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTGGCAACCGTga                                                                                                                                           <  1:153781/150‑2 (MQ=255)
          cAGCGCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTg                                                                                                                                         <  2:301832/150‑1 (MQ=255)
           aGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTgg                                                                                                                                        >  1:332511/1‑150 (MQ=255)
                gCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTg                                                                                                                                   <  1:888535/150‑1 (MQ=255)
                gCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTg                                                                                                                                   <  2:1406387/150‑1 (MQ=255)
                  cAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCACCAAGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGc                                                                                                                                 <  1:100734/150‑1 (MQ=255)
                  cAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGc                                                                                                                                 >  2:1081410/1‑150 (MQ=255)
                  cAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGc                                                                                                                                 >  2:1181642/1‑150 (MQ=255)
                       aGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTg                                                                                                                            <  1:1003579/150‑1 (MQ=255)
                         cAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGc                                                                                                                          <  1:804338/150‑1 (MQ=255)
                                  gCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGgc                                                                                                                  >  2:1306950/1‑150 (MQ=255)
                                     aCGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGgccgc                                                                                                               <  1:99534/150‑1 (MQ=255)
                                      cGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGgccgcc                                                                                                              >  2:632004/1‑150 (MQ=255)
                                            gTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGcc                                                                                                        <  2:360458/150‑1 (MQ=255)
                                             ttGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTg                                                                                                                            <  1:1338576/129‑1 (MQ=255)
                                             ttGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTg                                                                                                                            >  2:1338576/1‑129 (MQ=255)
                                               gAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCAtt                                                                                                     <  1:517246/150‑1 (MQ=255)
                                                  cAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTAcc                                                                                                  <  2:435875/150‑1 (MQ=255)
                                                   aaGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGt                                                                                                                                             <  1:382540/106‑1 (MQ=255)
                                                   aaGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGt                                                                                                                                             >  2:382540/1‑106 (MQ=255)
                                                     ggCACGGTACTGCTGAACCGTGGCCAGATGATGACCGCTTCGCAAACCGAGCTGCGCGATGCCGCTGCCGCGCCGGCAATACCACTGCTCAGGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCgg                                                                                               <  1:434412/150‑1 (MQ=255)
                                                      gCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGc                                                                                              >  1:1147149/1‑150 (MQ=255)
                                                        aCGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAg                                                                                            >  2:816563/1‑150 (MQ=255)
                                                          ggTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGat                                                                                          <  2:627755/150‑1 (MQ=255)
                                                               tgctgcAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGatca                                                                                        >  1:937071/1‑147 (MQ=255)
                                                               tgctgcAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGatca                                                                                        <  2:937071/147‑1 (MQ=255)
                                                                         cagccagATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTgcg                                                                           >  2:710796/1‑150 (MQ=255)
                                                                                atgatgACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAaca                                                                    >  2:220952/1‑150 (MQ=255)
                                                                                   atgaCCGCTTCGCAAACCGAGGGGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTg                                                                 <  1:384088/150‑1 (MQ=255)
                                                                                     gaCCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCt                                                                              >  1:148627/1‑135 (MQ=255)
                                                                                     gaCCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCt                                                                              <  2:148627/135‑1 (MQ=255)
                                                                                        cGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCctgct                                                                                                                     >  1:1191773/1‑93 (MQ=255)
                                                                                        cGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCctgct                                                                                                                     <  2:1191773/93‑1 (MQ=255)
                                                                                          cTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGcgccg                                                          <  2:1147149/150‑1 (MQ=255)
                                                                                           ttCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACAATACGcgccga                                                         >  2:153934/1‑149 (MQ=255)
                                                                                            tCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACAGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGCGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCa                                                        <  2:412908/150‑1 (MQ=255)
                                                                                              gCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCaga                                                      <  1:808228/150‑1 (MQ=255)
                                                                                                aaaCCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTTGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGc                                                    <  1:456669/150‑1 (MQ=255)
                                                                                                    cGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGNCAGATCNTCATTCCGCTGCNCCAAACACTGCGCGCCGCAGAGCgtgt                                                <  1:1241800/150‑1 (MQ=255)
                                                                                                    cGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCgtgt                                                <  2:399902/150‑1 (MQ=255)
                                                                                                          gcgcGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAg                                          >  1:113660/1‑150 (MQ=255)
                                                                                                              gATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGga                                      <  1:613605/150‑1 (MQ=255)
                                                                                                                 gccgctgccgcGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGAc                                   <  1:632004/150‑1 (MQ=255)
                                                                                                                   cgctgccgcGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCt                                 >  1:160413/1‑150 (MQ=255)
                                                                                                                    gctgccgcGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTg                                >  1:994932/1‑150 (MQ=255)
                                                                                                                    gctgccgcGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTg                                >  1:795967/1‑150 (MQ=255)
                                                                                                                     ctgccgcGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGa                               <  2:160413/150‑1 (MQ=255)
                                                                                                                      tgccgcGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGAc                              <  1:1306950/150‑1 (MQ=255)
                                                                                                                                 ccAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGc                   >  1:149598/1‑150 (MQ=255)
                                                                                                                                 ccAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTAGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGc                   <  1:816563/150‑1 (MQ=255)
                                                                                                                                  cAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTgc                                                                >  1:85488/1‑104 (MQ=255)
                                                                                                                                  cAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTgc                                                                <  2:85488/104‑1 (MQ=255)
                                                                                                                                  cAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCa                  >  2:612976/1‑150 (MQ=255)
                                                                                                                                      acacTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTa              >  1:756982/1‑150 (MQ=255)
                                                                                                                                           cagcACGGTGGCAACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGa                               >  2:142830/4‑128 (MQ=255)
                                                                                                                                           cagcACGGTGGCAACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGa                               <  1:142830/125‑1 (MQ=255)
                                                                                                                                           cagcACGGTGGCAACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTACAAcg         >  2:421839/4‑150 (MQ=255)
                                                                                                                                             tCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCgggg                                                                                                          <  2:1022575/51‑1 (MQ=255)
                                                                                                                                             tCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTggg                                                                                                          >  1:1022575/1‑51 (MQ=255)
                                                                                                                                               aCGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTACAAcgccgc     >  1:764581/1‑150 (MQ=255)
                                                                                                                                                 ggTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTACAACGCCGcga   <  1:710796/150‑1 (MQ=255)
                                                                                                                                                  gTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTACAACGCCGcgac  <  1:585632/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGCAATGCCCAGCTCGCCAGCGAGCAAGCGCTGCAACGCATGGTTGAACAAGGCACGGTACTGCTGCAGGCCAGCCAGATGATGACCGCTTCGCAAACCGAGGTGCGCGATGCCGCTGCCGCGCAGGCCAAAACACTGCTCACGGTTGCCACCGTGCTGGCACTGGCCCTTGGCCTGCTGGCCGCCTGGGCCATTACCCGGCAGATCATCATTCCGCTGCGCCAAACACTGCGCGCCGCAGAGCGTGTGGCCAGTGGAGACCTGACCCAGAGCCTGCAGGTACAACGCCGCGAC  >  NC_002947/5718083‑5718377

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