Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,927,503 G→A 9.5% H127Y (CAC→TAC) ‡ cheBA ← chemotaxis response regulator protein‑glutamate methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,927,5030GA 9.5% 179.7 / 6.6 63H127Y (CAC→TAC) ‡cheBAchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base A (3/3);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GCCGCCGGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATCAACACG  >  NC_002947/4927358‑4927651
                                                                                                                                                 |                                                                                                                                                    
gccgccGGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCAc                                                                                                                                                   >  3:374154/1‑149 (MQ=255)
     cgGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtct                                                                                                                                              >  5:80754/1‑149 (MQ=255)
      ggCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCACGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtctc                                                                                                                                             <  7:22607/149‑1 (MQ=255)
             gggCCGCGGCACGGGCAGGCGAGGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGTCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCTTAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACGTTCTCGCACAGc                                                                                                                                      <  8:275797/149‑1 (MQ=255)
               gCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCaa                                                                                                                                    >  6:84223/1‑149 (MQ=255)
               gCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCaa                                                                                                                                    >  4:151338/1‑149 (MQ=255)
                 cGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGTGCGAGAAGGTGCACACCATTTCGCGCAGCAAcc                                                                                                                                  <  6:59346/149‑2 (MQ=38)
                 cGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACt                                                                                                                                  >  2:185229/1‑149 (MQ=255)
                      cACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTg                                                                                                                             <  5:144578/149‑1 (MQ=255)
                      cACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTg                                                                                                                             <  6:112989/149‑1 (MQ=255)
                      cACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGGTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTg                                                                                                                             <  5:179765/149‑1 (MQ=255)
                        cGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGAc                                                                                                                           <  3:151338/149‑1 (MQ=255)
                           gCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCtt                                                                                                                        >  3:233320/1‑149 (MQ=255)
                             aGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGt                                                                                                                      <  1:296041/149‑1 (MQ=255)
                             aGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGt                                                                                                                      <  4:115380/149‑1 (MQ=255)
                              ggCGCAGACGCGCCCGGCTGCGCCTGGCTGGTAAAGGTGCTGGTCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGt                                                                                                                      >  4:233037/1‑148 (MQ=255)
                               gcgcTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCa                                                                                                                    <  4:374154/149‑1 (MQ=255)
                                  cTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAggg                                                                                                                 <  1:191125/149‑1 (MQ=255)
                                  cTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAggg                                                                                                                 <  6:80754/149‑1 (MQ=255)
                                   tGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtctc                                                                                                                                             <  8:144936/120‑1 (MQ=255)
                                   tGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtctc                                                                                                                                             >  7:144936/1‑120 (MQ=255)
                                        ggCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTgcgc                                                                                                           <  4:370025/149‑1 (MQ=255)
                                             gTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGAt                                                                                                      >  6:6687/1‑149 (MQ=255)
                                                  cgcCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTc                                                                                                   <  1:177343/147‑1 (MQ=255)
                                                        gctggTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGt                                                                                           >  5:351087/1‑149 (MQ=255)
                                                              aaaGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCg                                                                                     >  8:49225/1‑149 (MQ=255)
                                                                 ggTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAAc                                                                                                                                   >  2:21126/1‑100 (MQ=255)
                                                                 ggTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAAc                                                                                                                                   <  1:21126/100‑1 (MQ=255)
                                                                   tGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAgg                                                                                <  5:6687/149‑1 (MQ=255)
                                                                     cTGGCCGGCGCTGGCGCGGCGACCGTCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACNTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGt                                                                               <  5:100231/148‑1 (MQ=255)
                                                                     cTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAgg                                                                                >  4:150677/1‑147 (MQ=255)
                                                                        gCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTc                                                                           <  4:129821/149‑1 (MQ=255)
                                                                                tGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCga                                                                                                         <  5:356614/111‑1 (MQ=255)
                                                                                tGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCga                                                                                                         >  6:356614/1‑111 (MQ=255)
                                                                                  gcgcGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGcc                                                                  <  1:185229/148‑1 (MQ=255)
                                                                                  gcgcGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCg                                                                 <  5:84223/149‑1 (MQ=255)
                                                                                    gcggcgACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTgc                                                                                                             >  4:133294/1‑103 (MQ=255)
                                                                                       gcgACCGGCTCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATTGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTTCGCGAGATGTCATCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATc                                                            <  8:330642/149‑1 (MQ=255)
                                                                                            cGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATCTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAgcg                                                       >  2:235115/1‑149 (MQ=255)
                                                                                                      cTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAggg                                              <  4:233320/148‑1 (MQ=255)
                                                                                                        ggCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTACTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGa                                           >  4:203848/1‑149 (MQ=255)
                                                                                                         gCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGAc                                          >  3:335847/1‑149 (MQ=255)
                                                                                                         gCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGAc                                          >  4:236483/1‑149 (MQ=255)
                                                                                                           aTAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAgc                                                        <  1:310222/133‑1 (MQ=255)
                                                                                                           aTAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAgc                                                        >  2:310222/1‑133 (MQ=255)
                                                                                                               cTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGgc                                    <  3:157288/149‑1 (MQ=255)
                                                                                                                     aTCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTc                              <  3:236483/149‑1 (MQ=255)
                                                                                                                     aTCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTc                              <  6:351087/149‑1 (MQ=255)
                                                                                                                        cggcggTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAATTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCAt                            <  3:233037/148‑1 (MQ=255)
                                                                                                                          gcggTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATggg                         >  7:236022/1‑149 (MQ=255)
                                                                                                                             gTTGCTGTGCGAGAAGGTGCACACCATTTCGCGCAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGt                        >  2:94306/1‑147 (MQ=38)
                                                                                                                             gTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGt                        >  4:332626/1‑147 (MQ=255)
                                                                                                                                gctgcGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGc                    >  7:101843/1‑148 (MQ=255)
                                                                                                                                gctgcGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGAGGGCGCCTTCATGGGTCAGcg                   >  4:32011/1‑149 (MQ=255)
                                                                                                                                  tgtgcgAGAAGGTGCACACCATTTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGttttt                                                                                       <  2:25393/76‑1 (MQ=21)
                                                                                                                                  tgtgcgAGAAGGTGCACACCATTTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGttttt                                                                                       <  2:19543/76‑1 (MQ=21)
                                                                                                                                  tgtgcgAGAAGGTGCACACCATTTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGttttt                                                                                       >  1:25393/4‑79 (MQ=21)
                                                                                                                                  tgtgcgAGAAGGTGCACACCATTTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGttttt                                                                                       >  1:19543/4‑79 (MQ=21)
                                                                                                                                  tgcGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGAGCGTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATGCAGCGCGTCGAGGGTGTCGAGGGCGACTTCATGGGTCAGcgac                 >  6:282908/1‑149 (MQ=255)
                                                                                                                                      cgAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACgaga             <  4:293927/149‑1 (MQ=255)
                                                                                                                                           tGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACAt         >  2:225284/1‑148 (MQ=255)
                                                                                                                                           tGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATc        <  3:332626/149‑1 (MQ=255)
                                                                                                                                           tGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATc        <  8:236022/149‑1 (MQ=255)
                                                                                                                                                 gCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATCAACACg  >  4:161406/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCGCCGGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATCAACACG  >  NC_002947/4927358‑4927651

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