Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,420,641 Δ1 bp 8.3% 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. 8.3% 179.8 / 10.3 60intergenic (+81/+26)ctpH/PP_t62methyl‑accepting chemotaxis protein CtpH/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (35/19);  new base . (2/3);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CAGGCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCTTG  >  NC_002947/2420492‑2420788
                                                                                                                                                     |                                                                                                                                                    
cAGGCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCGCACAACACAAAGCCCGca                                                                                                                                                      <  2:295685/150‑1 (MQ=255)
   gCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTAGGGGGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGTACGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACAt                                                                                                                                                   <  2:101727/150‑1 (MQ=255)
           gcgcCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGC‑CATGTGCGGGCt                                                                                                                                          <  1:1014117/150‑1 (MQ=255)
           gcgcCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGC‑CATGTGCGGGCt                                                                                                                                          >  1:101028/1‑150 (MQ=255)
             gcCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCttt                                                                                                                                        <  2:293057/150‑1 (MQ=255)
               cTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTtgt                                                                                                                                      <  1:893383/150‑1 (MQ=255)
               cTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGC‑CATGTGCGGGCTTtgt                                                                                                                                      <  1:954348/150‑1 (MQ=255)
                    tGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAACGCCCGCTCAGG‑GCGGGCTTTGTGACaa                                                                                                                                 <  2:114680/150‑1 (MQ=255)
                    tGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGC‑CATGTGCGGGCTTtgtgttgt                                                                                                                                 <  2:101028/150‑5 (MQ=255)
                     ggTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAg                                                                                                                                >  2:296719/1‑150 (MQ=255)
                       tttGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGAGTTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGa                                                                                                                              >  1:281526/1‑150 (MQ=255)
                           tCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTtgt                                                                                                                                      <  2:779293/138‑1 (MQ=255)
                           tCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTtgt                                                                                                                                      >  1:779293/1‑138 (MQ=255)
                              gggCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGt                                                                                                                        >  2:185710/1‑149 (MQ=255)
                              gggCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGt                                                                                                                        <  1:185710/149‑1 (MQ=255)
                               ggCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTc                                                                                                                       >  2:657597/1‑149 (MQ=255)
                               ggCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTc                                                                                                                       <  1:657597/149‑1 (MQ=255)
                                           aaGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGAt                                                                                                          <  1:231178/150‑1 (MQ=255)
                                              tCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCg                                                                                                       >  1:362169/1‑150 (MQ=255)
                                                aCCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACaa                                                                                                                                 <  1:232861/122‑1 (MQ=255)
                                                aCCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACaa                                                                                                                                 >  2:232861/1‑122 (MQ=255)
                                                           gAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCtt                                                                                          <  1:296719/150‑1 (MQ=255)
                                                            aGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTtc                                                                                         <  2:362169/150‑1 (MQ=255)
                                                                 cTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGAt                                                                                                          <  1:561418/128‑1 (MQ=255)
                                                                 cTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGAt                                                                                                          >  2:561418/1‑128 (MQ=255)
                                                                 cTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTc                                                                                    >  1:761397/1‑150 (MQ=255)
                                                                  tGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTcc                                                                                   >  1:126508/1‑150 (MQ=255)
                                                                        ggTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACg                                                                             >  1:630466/1‑150 (MQ=255)
                                                                         gTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTAGACCCGAACgg                                                                            >  2:31911/1‑149 (MQ=255)
                                                                         gTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATAAATCGAACACAACACAAAGCTGGCACATT‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACga                                                                            >  1:157626/1‑150 (MQ=255)
                                                                            tGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGAc                                                                         >  2:232739/1‑150 (MQ=255)
                                                                              cAGGAGTTGCAGGTGCCATGGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTgg                                                                                                                         <  1:546667/100‑1 (MQ=255)
                                                                              cAGGAGTTGCAGGTGCCATGGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTgg                                                                                                                         >  2:546667/1‑100 (MQ=255)
                                                                              cAGGAGTTGCAGGTGCCATAGAGGTTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGAcgc                                                                       >  1:1185256/1‑150 (MQ=255)
                                                                              cAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGAcgc                                                                       >  1:97455/1‑150 (MQ=255)
                                                                              cAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCAGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGAcgc                                                                       >  2:813067/1‑150 (MQ=255)
                                                                                      gCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGATTTGTGACAAGGATTTGGTCGACGCGACCCGACGTTTACCGTCTGCTTTCTTGTCCCGAACGAGACGCGCTACCTg                                                               >  1:1248862/1‑150 (MQ=255)
                                                                                        aGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGc                                                             >  2:909250/1‑150 (MQ=255)
                                                                                          gTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTg                                                           >  2:827229/1‑150 (MQ=255)
                                                                                          gTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTg                                                           >  1:991566/1‑150 (MQ=255)
                                                                                            gCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTgcg                                                         >  1:1069055/1‑150 (MQ=255)
                                                                                             ccATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTgcgc                                                        <  2:126508/150‑1 (MQ=255)
                                                                                               aTAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTa                                                      >  2:1015467/1‑150 (MQ=255)
                                                                                                 agagGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTaca                                                    <  2:630466/150‑1 (MQ=255)
                                                                                                          aGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGa                                           >  1:490967/1‑150 (MQ=255)
                                                                                                                  aCTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATc                                   >  1:209455/1‑150 (MQ=255)
                                                                                                                      aCGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGAt                                             >  2:42586/1‑136 (MQ=255)
                                                                                                                      aCGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGAt                                             <  1:42586/136‑1 (MQ=255)
                                                                                                                      aCGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGCCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTAc                               >  2:20482/1‑150 (MQ=255)
                                                                                                                       cGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACt                              <  2:1185256/150‑1 (MQ=255)
                                                                                                                        gTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTg                             >  2:164144/1‑150 (MQ=255)
                                                                                                                          aTTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGCCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCa                           <  1:20482/150‑1 (MQ=255)
                                                                                                                            ttttCGCACACAACACAAAGGCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATc                         >  2:1230094/1‑150 (MQ=255)
                                                                                                                                     ttgtcacaaAGCCCGC‑CATGTGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTg                >  1:446473/5‑150 (MQ=255)
                                                                                                                                         cacaaAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCg            >  2:1109021/1‑150 (MQ=255)
                                                                                                                                          acaaAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCGCGTTTTTGCTCGc           >  2:257986/1‑150 (MQ=255)
                                                                                                                                          acaaAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCATACTGAATCCCGTTTTTGCTccc           >  2:713277/1‑148 (MQ=255)
                                                                                                                                           caaAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCa          >  2:796016/1‑150 (MQ=255)
                                                                                                                                            aaaGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCaa         <  1:1109021/150‑1 (MQ=255)
                                                                                                                                                   gCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCttg  <  2:490967/150‑1 (MQ=255)
                                                                                                                                                   gCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCttg  >  2:1156239/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATG‑GCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCTTG  >  NC_002947/2420492‑2420788

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