Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,927,508 T→A 9.2% K125M (AAG→ATG)  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,5080TA 9.2% 189.6 / 7.5 65K125M (AAG→ATG) cheBAchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (3/3);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATCAACACGGGCG  >  NC_002947/4927363‑4927655
                                                                                                                                                 |                                                                                                                                                   
cgGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtct                                                                                                                                                  >  5:80754/1‑149 (MQ=255)
 ggCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCACGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtctc                                                                                                                                                 <  7:22607/149‑1 (MQ=255)
        gggCCGCGGCACGGGCAGGCGAGGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGTCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCTTAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACGTTCTCGCACAGc                                                                                                                                          <  8:275797/149‑1 (MQ=255)
          gCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCaa                                                                                                                                        >  4:151338/1‑149 (MQ=255)
          gCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCaa                                                                                                                                        >  6:84223/1‑149 (MQ=255)
            cGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGTGCGAGAAGGTGCACACCATTTCGCGCAGCAAcc                                                                                                                                      <  6:59346/149‑2 (MQ=38)
            cGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACt                                                                                                                                      >  2:185229/1‑149 (MQ=255)
                 cACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTg                                                                                                                                 <  6:112989/149‑1 (MQ=255)
                 cACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTg                                                                                                                                 <  5:144578/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                                                                                                                     <  6:80754/149‑1 (MQ=255)
                             cTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAggg                                                                                                                     <  1:191125/149‑1 (MQ=255)
                              tGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtctc                                                                                                                                                 >  7:144936/1‑120 (MQ=255)
                              tGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTtctc                                                                                                                                                 <  8:144936/120‑1 (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                                              >  4:236483/1‑149 (MQ=255)
                                                                                                    gCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGAc                                              >  3:335847/1‑149 (MQ=255)
                                                                                                      aTAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAgc                                                            >  2:310222/1‑133 (MQ=255)
                                                                                                      aTAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAgc                                                            <  1:310222/133‑1 (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)
                                                                                                                                             cACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATCAACACgg     <  3:150677/149‑1 (MQ=255)
                                                                                                                                                cTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTCCGCCGGCATCCAGCCCGTCGAGGGTGAGGCGGGCGACTTCATGGGTCTGAGACGAGAACATCAACACGGGCg  >  7:302992/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGCGTCGGGGCCGCTGCACGGGCAGGCGCTGCCGGGCGCGTCTGCGCCTGGCTGGTAAAGGTGCTGGCCGGCGCTGGCGCGGCGACCGGCGCCGTTCTGGCATAGCTGCCGATCCGGCGGTTGCTGCGCGAAATGGTGTGCACCTTCTCGCACAGCAACTGCTTGACCTTGTCAGGGTTGCGCGAGATGTCCTCGAAGTTTTTCGGCAGGTAGTCCACTGCGCCGGCATCCAGCGCGTCGAGGGTGACGCGGGCGCCTTCATGGGTCAGCGACGAGAACATCAACACGGGCG  >  NC_002947/4927363‑4927655

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: ‑