Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,709 G→T 20.4% H273N (CAC→AAC)  ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,7090GT20.4% 432.9 / 102.1 186H273N (CAC→AAC) ftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (77/71);  new base T (18/20);  total (95/91)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 188 total aligned reads displayed.

CACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAA  >  NC_002947/5367560‑5367857
                                                                                                                                                     |                                                                                                                                                    
 accTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGa                                                                                                                                                     <  1:2062629/150‑1 (MQ=255)
        ggTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACcgccgc                                                                                                                                              >  2:1091687/1‑150 (MQ=255)
            gAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCAt                                                                                                                                          >  2:325701/1‑150 (MQ=255)
            gAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCAt                                                                                                                                          >  1:2563905/1‑150 (MQ=255)
            gAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCAt                                                                                                                                          >  1:1726420/1‑150 (MQ=255)
            gAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCAt                                                                                                                                          >  2:2486170/1‑150 (MQ=255)
               gCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGcc                                                                                                                                       <  2:462016/150‑1 (MQ=255)
                cGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCg                                                                                                                                      <  1:325627/150‑1 (MQ=255)
                      gggCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCa                                                                                                                                >  2:3994835/1‑150 (MQ=255)
                           cagctgCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATg                                                                                                                           <  1:2330285/150‑1 (MQ=255)
                           cagcagCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATg                                                                                                                           <  2:4447110/150‑1 (MQ=255)
                            agcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATgg                                                                                                                          <  2:4617105/150‑1 (MQ=255)
                            agcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATgg                                                                                                                          <  2:4765863/150‑1 (MQ=255)
                              cagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCg                                                                                                                        >  1:191058/1‑150 (MQ=255)
                                  gcAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCa                                                                                                                    >  1:3616270/1‑150 (MQ=255)
                                       gTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGc                                                                                                               <  2:2901840/150‑1 (MQ=255)
                                         cGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGt                                                                                                             <  1:3295055/150‑1 (MQ=255)
                                          gAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTc                                                                                                            >  2:223517/1‑150 (MQ=255)
                                            gCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGa                                                                                                          >  2:2676662/1‑150 (MQ=255)
                                             cacaTCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGgc                                                                                                                                    <  2:2640603/123‑1 (MQ=255)
                                                 tCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGTGACCAACGGCGTCGATTtcg                                                                                                     >  2:4580661/1‑150 (MQ=255)
                                                 tCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtcg                                                                                                     <  1:1473948/150‑1 (MQ=255)
                                                   cGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtcgtc                                                                                                   <  1:4237883/150‑1 (MQ=255)
                                                     gACGGTTGGTGGCGGCTATGACGATGATGCCATCATTCATCTCGAAGCCTTCCATCTCCACCAGCAACTTGTTGAGGGTTTTCTCACGCTCGTCGTGACCGCCGCCCATTCCGGCTCCACGATGGCGACCAACGGCGTCGATTTCGTCGa                                                                                                 <  1:1149183/150‑1 (MQ=255)
                                                        ggTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGa                                                                                              >  1:1259437/1‑150 (MQ=255)
                                                          ttGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAg                                                                                            >  2:3249980/1‑150 (MQ=255)
                                                                   gcgATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCGAGTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCAc                                                                                   >  2:1077757/1‑150 (MQ=255)
                                                                       tgacgatgaTGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACgggg                                                                               <  2:2585580/150‑1 (MQ=255)
                                                                             tgatgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCgccgcc                                                                                                                                             >  1:960616/1‑82 (MQ=255)
                                                                              gatgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGctt                                                                        >  2:1348390/1‑150 (MQ=255)
                                                                                 gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                                                     >  1:2122955/1‑150 (MQ=255)
                                                                                 gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                                                     >  1:3895399/1‑150 (MQ=255)
                                                                                  ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                    <  1:140041/150‑1 (MQ=255)
                                                                                    atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                                                  >  1:2751430/1‑150 (MQ=255)
                                                                                       attcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTg                                                               >  1:1575992/1‑150 (MQ=255)
                                                                                          catCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTc                                                            <  1:1659784/150‑1 (MQ=255)
                                                                                           atCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCg                                                           >  2:751045/1‑150 (MQ=255)
                                                                                            tctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGa                                                          >  2:2138295/1‑150 (MQ=255)
                                                                                             ctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGaa                                                         <  1:1590878/150‑1 (MQ=255)
                                                                                             ctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGaa                                                         <  2:4724830/150‑1 (MQ=255)
                                                                                              tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                                        <  1:1865643/150‑1 (MQ=255)
                                                                                                gAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAt                                                      >  2:3519665/1‑150 (MQ=255)
                                                                                                gAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAt                                                      >  1:2484175/1‑150 (MQ=255)
                                                                                                 aaaccatCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATg                                                     >  1:830161/7‑150 (MQ=255)
                                                                                                     catCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcac                                                 <  2:2259156/148‑3 (MQ=255)
                                                                                                     catCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcac                                                 <  1:585669/148‑3 (MQ=255)
                                                                                                     catCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcac                                                 >  2:3676098/3‑148 (MQ=255)
                                                                                                      atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                >  1:5042428/2‑150 (MQ=255)
                                                                                                      atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                >  1:2969049/2‑150 (MQ=255)
                                                                                                       tCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAATATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacgg                                               >  1:3194783/1‑149 (MQ=255)
                                                                                                        ccATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacgga                                              <  1:3095643/150‑3 (MQ=255)
                                                                                                          atttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgacg                                            >  1:3948617/7‑149 (MQ=255)
                                                                                                           tttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgacgc                                           >  1:374883/6‑148 (MQ=255)
                                                                                                            ttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGAcgcg                                          <  1:1116798/146‑1 (MQ=255)
                                                                                                             tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGACCATGTCGCGCACACgg                                         >  2:3383527/1‑150 (MQ=255)
                                                                                                             tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACgg                                         <  1:4222841/150‑1 (MQ=255)
                                                                                                             tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACgg                                         >  2:2940415/1‑150 (MQ=255)
                                                                                                               gaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCTGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTg                                                                    <  1:4939635/120‑1 (MQ=255)
                                                                                                               gaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCTGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTg                                                                    >  2:4939635/2‑121 (MQ=255)
                                                                                                               caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtc                                                                                                      >  2:1270551/1‑87 (MQ=255)
                                                                                                               caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                                       <  1:668045/150‑1 (MQ=255)
                                                                                                               caccaGCAACTGGGTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAGGTCGCGCACACGGCt                                       >  2:366128/1‑150 (MQ=255)
                                                                                                                  cagaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCTAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGc                                    >  1:410381/7‑150 (MQ=255)
                                                                                                                  caGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGc                                    <  2:2417946/150‑1 (MQ=255)
                                                                                                                  caGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGc                                    >  1:4363175/1‑150 (MQ=255)
                                                                                                                   agaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCa                                   <  1:3837570/145‑1 (MQ=255)
                                                                                                                   aGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGACGCGCACACGGCTGGCa                                   <  1:499561/150‑1 (MQ=255)
                                                                                                                    gaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                  >  2:825789/5‑150 (MQ=255)
                                                                                                                     cAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCAcc                                 <  2:4563802/150‑1 (MQ=255)
                                                                                                                       aCTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGa                               >  1:474135/1‑150 (MQ=255)
                                                                                                                           gTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCCAGGCGCATTCTGCTTGGCCTGTTCGAGCATGTCACGGACGCGGCGTGCAcccgagcg                           >  1:5204631/1‑144 (MQ=255)
                                                                                                                           gTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACAcc                           >  1:4868992/1‑150 (MQ=255)
                                                                                                                            ttGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAccc                          >  2:637096/1‑150 (MQ=255)
                                                                                                                              gAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgatgat                                                                                       >  1:1341336/1‑87 (MQ=255)
                                                                                                                              gAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCAc                        >  2:194270/1‑150 (MQ=255)
                                                                                                                                    ttGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACAt                  <  2:3277994/150‑1 (MQ=255)
                                                                                                                                     tGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACAtt                 <  2:116584/150‑1 (MQ=255)
                                                                                                                                      gCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACAttt                >  1:4963186/1‑150 (MQ=255)
                                                                                                                                       ttcgcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTc               <  2:2960970/143‑1 (MQ=255)
                                                                                                                                       cTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                                       <  2:927690/126‑1 (MQ=255)
                                                                                                                                       cTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                                       >  1:927690/1‑126 (MQ=255)
                                                                                                                                        tcgcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTca              <  1:761373/144‑2 (MQ=255)
                                                                                                                                        tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCAc                                                                                   >  2:3352449/1‑81 (MQ=255)
                                                                                                                                        tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCAc                                                                                   <  1:3352449/81‑1 (MQ=255)
                                                                                                                                        tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc              >  2:5380843/1‑150 (MQ=255)
                                                                                                                                         cgcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTcaa             <  2:2828095/145‑3 (MQ=255)
                                                                                                                                         cACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCa             <  2:4190604/150‑1 (MQ=255)
                                                                                                                                         cACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCa             <  1:3682972/150‑1 (MQ=255)
                                                                                                                                          gcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                  >  1:303250/5‑128 (MQ=255)
                                                                                                                                          gcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                  <  2:303250/124‑1 (MQ=255)
                                                                                                                                          aCGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCAc            >  1:4333694/1‑150 (MQ=255)
                                                                                                                                          aCGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCAAc            >  2:4128347/1‑150 (MQ=255)
                                                                                                                                            gCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGa          >  1:4609798/1‑150 (MQ=255)
                                                                                                                                             ttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACaaa         <  1:1531178/149‑1 (MQ=255)
                                                                                                                                             cTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGaa         <  2:1930744/150‑1 (MQ=255)
                                                                                                                                             cTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGaa         <  2:3949210/150‑1 (MQ=255)
                                                                                                                                                gtcgtTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTc      <  1:2783871/150‑1 (MQ=255)
                                                                                                                                                gtcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTc      >  2:1760281/1‑150 (MQ=255)
                                                                                                                                                 tcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATCATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGtaa     <  2:3088046/150‑3 (MQ=255)
                                                                                                                                                    tGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGaa  >  2:4964605/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAA  >  NC_002947/5367560‑5367857

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