Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,420,641 Δ1 bp 7.6% intergenic (+81/+26) ctpH → / ← PP_t62 methyl‑accepting chemotaxis protein CtpH/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,420,6410A. 7.6% 202.9 / 10.9 66intergenic (+81/+26)ctpH/PP_t62methyl‑accepting chemotaxis protein CtpH/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base . (3/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

GCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCTTGTT  >  NC_002947/2420505‑2420790
                                                                                                                                        |                                                                                                                                                      
gcCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACCGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCttt                                                                                                                                          <  1:1304235/150‑1 (MQ=255)
gcCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACCGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCttt                                                                                                                                          >  2:1304235/1‑150 (MQ=255)
  cTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTtgt                                                                                                                                        >  1:666126/1‑150 (MQ=255)
  cTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTtgt                                                                                                                                        <  2:485879/150‑1 (MQ=255)
    gtaTGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGGCTCGTTCGGGAGGAGTTGCAGTTGCTATAGAGGCTGAAGCCGCTCACGGACGTGTTTTTCTCACACAGCTCAATGCCCGCACAGTGCGGGTTTTTTGa                                                                                                                                      <  2:18963/148‑1 (MQ=255)
       tGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACaa                                                                                                                                   <  2:774411/150‑1 (MQ=255)
           ttGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGAt                                                                                                                               <  1:1436261/150‑1 (MQ=255)
           ttGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGAt                                                                                                                               <  2:1579000/150‑1 (MQ=255)
             gTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGATTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGAttt                                                                                                                             >  2:1770991/1‑150 (MQ=255)
             gTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGAttt                                                                                                                             >  2:1871379/1‑150 (MQ=255)
                 gggCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTAGTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTc                                                                                                                         <  2:1512386/150‑1 (MQ=255)
                    cAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGga                                                                                                                      <  2:2313590/150‑1 (MQ=255)
                      gTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGgaga                                                                                                                    >  2:857226/1‑150 (MQ=255)
                             gAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGcac                                                                                                                                                       >  2:2029354/1‑109 (MQ=255)
                             gAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGcac                                                                                                                                                       <  1:2029354/109‑1 (MQ=255)
                                gTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCTCTCACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTc                                                                                                          >  2:645597/1‑150 (MQ=255)
                                gTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTc                                                                                                          <  2:856302/150‑1 (MQ=255)
                                   aCCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGaa                                                                                                       >  1:331502/1‑150 (MQ=255)
                                     cTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAAcc                                                                                                     >  1:183681/1‑150 (MQ=255)
                                        ccGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGc                                                                                                  <  1:857226/150‑1 (MQ=255)
                                          gCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCgg                                                                                                <  1:1954025/150‑1 (MQ=255)
                                           ccTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCggg                                                                                                                                              <  2:260790/103‑1 (MQ=255)
                                           ccTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCggg                                                                                                                                              >  1:260790/1‑103 (MQ=255)
                                                 tGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTtctc                                                                                         <  2:666126/150‑1 (MQ=255)
                                                    cTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGAc                                                                                                                 >  2:1218370/1‑123 (MQ=255)
                                                    cTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGAc                                                                                                                 <  1:1218370/123‑1 (MQ=255)
                                                     tGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTcc                                                                                     <  1:1011519/150‑1 (MQ=255)
                                                         tCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGaa                                                                                 >  2:1932731/1‑150 (MQ=255)
                                                           ggTCTGAAGGAGTTGCAGGTGCCATAGAGGATAACAAGGATGGCTGGTATACGTTTTTCACACAACACACAGCCTGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACg                                                                               <  1:522678/150‑1 (MQ=255)
                                                                gCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAg                                                                                                               >  1:39002/1‑113 (MQ=255)
                                                                gCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAg                                                                                                               <  2:39002/113‑1 (MQ=255)
                                                                 cAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGgagaga                                                                                                                  >  2:603325/1‑109 (MQ=255)
                                                                 cAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGgagaga                                                                                                                  <  1:603325/109‑1 (MQ=255)
                                                                 cAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGAcgc                                                                         <  1:1932731/150‑1 (MQ=255)
                                                                   ggAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACgcgc                                                                       >  2:1317399/1‑150 (MQ=255)
                                                                       ttGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGga                                                                                                                      >  1:800362/1‑99 (MQ=255)
                                                                       ttGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGga                                                                                                                      <  2:800362/99‑1 (MQ=255)
                                                                       ttGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGCTTtgt                                                                                                                                        <  2:2111503/81‑1 (MQ=38)
                                                                       ttGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGCTTtgt                                                                                                                                        >  1:2111503/1‑81 (MQ=38)
                                                                                ccATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGTTGCCTGGCTgcgc                                                          >  1:479551/1‑150 (MQ=255)
                                                                                  aTAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTa                                                        <  1:1317399/150‑1 (MQ=255)
                                                                                     gagGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCg                                                                                        <  1:311488/115‑1 (MQ=255)
                                                                                     gagGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCg                                                                                        >  2:311488/1‑115 (MQ=255)
                                                                                        gCTGAAGCCGCGTACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGATACGCGCTACCTGGCTGCGCTACACTcc                                                  >  1:697281/1‑150 (MQ=255)
                                                                                                   gcACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGCCAGGACTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACaaaa                                       >  1:1791329/1‑150 (MQ=255)
                                                                                                    cACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGa                                                                                                             >  1:1504148/1‑79 (MQ=255)
                                                                                                    cACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGGCGGAGAGACAGGa                                                                                                             <  2:1504148/79‑1 (MQ=255)
                                                                                                       tGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACga                                                                              <  1:1993452/107‑1 (MQ=255)
                                                                                                       tGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACga                                                                              >  2:1993452/1‑107 (MQ=255)
                                                                                                          cGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACt                                <  2:331502/150‑1 (MQ=255)
                                                                                                                    gCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACCCCGACGAACAAACTCCTACTGCATCCCGtt                      >  2:137209/1‑150 (MQ=255)
                                                                                                                     cacacaACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACGAAATCCTACTGCATCCCGttt                     >  2:1446590/1‑150 (MQ=255)
                                                                                                                        acaacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATc                                     >  1:1319181/1‑131 (MQ=255)
                                                                                                                        acaacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATc                                     <  2:1319181/131‑1 (MQ=255)
                                                                                                                         caacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGc                 >  2:1821570/1‑150 (MQ=255)
                                                                                                                          aacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGt                                                                                                                          >  1:633724/1‑44 (MQ=255)
                                                                                                                          aacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGt                                                                                                                          <  2:633724/44‑1 (MQ=255)
                                                                                                                           acacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTc               >  1:800465/1‑150 (MQ=255)
                                                                                                                            cacaaAGCCCGCCATGTGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTgcgc                                                          >  2:2115301/1‑106 (MQ=39)
                                                                                                                            cacaaAGCCCGCCATGTGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTgcgc                                                          <  1:2115301/106‑1 (MQ=39)
                                                                                                                             acaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGc             >  2:1914495/1‑150 (MQ=255)
                                                                                                                             acaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGc             >  2:1198142/1‑150 (MQ=255)
                                                                                                                             acaaAGCCCGCCATGTGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGc             >  1:137797/1‑150 (MQ=255)
                                                                                                                              caaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATg                                              >  2:1948218/1‑116 (MQ=255)
                                                                                                                              caaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATg                                              <  1:1948218/116‑1 (MQ=255)
                                                                                                                               aaaGCCCGCACATGGCGGGCTTTGTGACAGTGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCaa           >  1:1593067/1‑150 (MQ=255)
                                                                                                                                     cGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCtt     >  1:1540362/1‑150 (MQ=255)
                                                                                                                                      gCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCttg    >  1:1500318/1‑150 (MQ=255)
                                                                                                                                       cacaTGGCGGGCTTTGTGACAACGATTTCGTCGGAGACACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCCAGCCCttgt   <  2:183681/150‑1 (MQ=255)
                                                                                                                                        acaTGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCttgtt  <  1:1914495/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                      
GCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCTTGTT  >  NC_002947/2420505‑2420790

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: