Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 516,912 C→T 7.2% R27H (CGC→CAC)  PP_0429 ← short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947516,9120CT 7.2% 202.3 / 4.4 69R27H (CGC→CAC) PP_0429short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/3);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTTTG  >  NC_002947/516766‑517058
                                                                                                                                                   |                                                                                                                                                  
cttcCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTgc                                                                                                                                                    <  3:253919/147‑1 (MQ=255)
   ccAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGCCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCCGCGTATCTCGATGGCGTGCCACCAGGATCAGATTGCGCCcg                                                                                                                                               >  8:358032/1‑149 (MQ=255)
     aGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGATCAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCcgcg                                                                                                                                             <  8:250287/149‑1 (MQ=255)
     aGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCACCcgcg                                                                                                                                             >  7:229150/1‑149 (MQ=255)
        tctAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCg                                                                                                                                          >  1:23925/1‑149 (MQ=255)
          tAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTg                                                                                                                                        >  1:366088/1‑149 (MQ=255)
            gCACGAAGCCTGAAAGGCGCAGCGGCTGACTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGt                                                                                                                                         <  6:305400/146‑1 (MQ=255)
             cACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTCGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCCTTGGCGTGCCACCAGGAGCAGATTGCGCCCGCCCCGTGCCa                                                                                                                                     >  5:83741/1‑149 (MQ=255)
                     cTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAAGAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTcc                                                                                                                             <  2:237103/149‑1 (MQ=255)
                     cTGAAAGGCGCAGCGGCTGGCTGAGGGCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTcc                                                                                                                             <  3:89031/149‑1 (MQ=255)
                      tGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTGTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCg                                                                                                                            >  1:323428/1‑149 (MQ=255)
                      tGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCg                                                                                                                            >  5:205130/1‑149 (MQ=255)
                          aGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCaa                                                                                                                        <  7:291263/149‑1 (MQ=255)
                              gcAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAaggcc                                                                                                                    <  2:230765/149‑1 (MQ=255)
                              gcAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAaggc                                                                                                                     >  6:306548/1‑148 (MQ=255)
                                 gCGGCTGGCTGATGTCGCAGGCGCGGAACAGCACCTCGTCGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGTGCCCGCGCCGTGCCAGCGCTTCCGCCAatgccagg                                                                                                                 <  4:285033/149‑1 (MQ=255)
                                  cGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAaggccaggc                                                                                                                <  3:173188/149‑1 (MQ=255)
                                       ggctgAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGccaggcccagg                                                                                                           >  8:268638/1‑149 (MQ=255)
                                        gctgAGGTCGCATGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGccaggcccag                                                                                                            <  8:196591/147‑1 (MQ=255)
                                          tgAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCCGGCCCAGGCCg                                                                                                        >  2:303387/1‑149 (MQ=255)
                                             ggTCGCAGGCGCTGAAAAGCAACTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTg                                                                                                     <  8:171306/149‑1 (MQ=255)
                                             ggTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTg                                                                                                     <  8:319784/149‑1 (MQ=255)
                                                cGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGaa                                                                                                  <  6:178379/149‑1 (MQ=255)
                                                 gCAGGCGCGGAACAGCACCTCGACGCCGTATCGCTGGGGCAGGTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGTTCAGTTTGCGCCCGCGGCGTGCCAGCGCTTCCGCCAAGGGCAGGGCGAGGCCGCTGGaa                                                                                                  <  6:83741/148‑1 (MQ=255)
                                                 gCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAg                                                                                                 <  5:306548/149‑1 (MQ=255)
                                                       gcgGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAg                                                                                           >  5:84624/1‑149 (MQ=255)
                                                         ggAACAGCACCTCGACGCCGAAACTCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGt                                                                                          >  5:174309/1‑148 (MQ=255)
                                                                cACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATgg                                                                                  >  1:330559/1‑149 (MQ=255)
                                                                 aCCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGTGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCAt                                                                                    <  6:314446/146‑1 (MQ=255)
                                                                 aCCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGTGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCAt                                                                                    <  2:242219/146‑1 (MQ=255)
                                                                 aCCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGTGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCAt                                                                                    >  5:314446/1‑146 (MQ=255)
                                                                 aCCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGTGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCAt                                                                                    >  1:242219/1‑146 (MQ=255)
                                                                     cGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTaa                                                                             >  7:10922/1‑149 (MQ=255)
                                                                        cGCCGAAACTCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACgg                                                                          <  6:174309/149‑1 (MQ=255)
                                                                        cGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACgg                                                                          <  1:226521/149‑1 (MQ=255)
                                                                           cGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGt                                                                        >  2:34097/1‑148 (MQ=255)
                                                                            gAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCa                                                                      <  2:23925/149‑1 (MQ=255)
                                                                             aaaCGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCACCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCAt                                                                     <  8:229150/149‑1 (MQ=255)
                                                                                 gCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCaa                                                                 <  2:246827/149‑1 (MQ=255)
                                                                                 gCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCaa                                                                 <  6:114959/149‑1 (MQ=255)
                                                                                       tCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTcc                                                           >  5:321298/1‑149 (MQ=255)
                                                                                           cTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTggg                                                       <  3:35279/149‑1 (MQ=255)
                                                                                            tCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTgggg                                                      >  3:102750/1‑149 (MQ=255)
                                                                                               aTGGCAACAGGCTCCAGCGTTTCTCGTTGGCGGGCCACCAGTATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGGTGGAATCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCg                                                   <  3:27513/149‑1 (MQ=255)
                                                                                               aTGGCAACAGGCTCCAGCGTTTCTAGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCg                                                   <  5:301819/149‑1 (MQ=255)
                                                                                                 ggCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCgg                                                  >  2:168510/1‑148 (MQ=255)
                                                                                                    aaCAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGaa                                                >  7:16682/1‑147 (MQ=255)
                                                                                                           tCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCAt                                                                     >  1:53058/1‑120 (MQ=255)
                                                                                                             cAGCGTTTTTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTa                                     <  3:285466/149‑1 (MQ=255)
                                                                                                             cAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTGCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTTACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTa                                     <  7:358032/149‑1 (MQ=255)
                                                                                                             cAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTa                                     <  7:106964/149‑1 (MQ=255)
                                                                                                              aGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAg                                    >  6:21667/1‑149 (MQ=255)
                                                                                                              aGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAg                                    <  4:102750/149‑1 (MQ=255)
                                                                                                               gCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAgt                                   >  1:76947/1‑149 (MQ=255)
                                                                                                                  tttCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTa                                <  2:330559/149‑1 (MQ=255)
                                                                                                                      tcGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAAc                            >  5:313516/1‑149 (MQ=255)
                                                                                                                        gTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACtt                          <  6:313516/149‑1 (MQ=255)
                                                                                                                            gCGTGCCACCAGGATCAGCTTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGccc                      <  5:317088/149‑1 (MQ=255)
                                                                                                                            gCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGccc                      <  2:366088/149‑1 (MQ=255)
                                                                                                                            gCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACCGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTTccc                      >  5:327120/1‑149 (MQ=255)
                                                                                                                                 ccaccaGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGc                  >  5:12674/1‑148 (MQ=255)
                                                                                                                                   accaGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAg               <  2:188031/149‑1 (MQ=255)
                                                                                                                                        gATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGtt          >  2:286438/1‑149 (MQ=255)
                                                                                                                                         aTCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTg         >  1:114979/1‑149 (MQ=255)
                                                                                                                                         aTCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTg         >  3:20403/1‑149 (MQ=255)
                                                                                                                                           cAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGcc       >  6:398596/1‑149 (MQ=255)
                                                                                                                                            aGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGcc       <  8:16682/148‑1 (MQ=255)
                                                                                                                                               aTTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCttt    <  6:321298/149‑1 (MQ=255)
                                                                                                                                                ttGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCtttt   <  6:12674/149‑1 (MQ=255)
                                                                                                                                                ttGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCtttt   <  5:398596/149‑1 (MQ=255)
                                                                                                                                                  gcgcCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTTTg  <  5:21667/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGACGCCGAAACGCTGGGTCAGCTCGATGGCAACAGGCTCCAGCGTTTCTCGTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTTTG  >  NC_002947/516766‑517058

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: