Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,241,714:1 +C 7.1% coding (158/633 nt) hmgC ← maleylacetoacetate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,241,7141.C 7.1% 275.7 / 12.8 70coding (158/633 nt)hmgCmaleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base . (28/37);  new base C (2/3);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGT  >  NC_002947/5241566‑5241859
                                                                                                                                                     |                                                                                                                                                 
gCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGGACCAGCACCTCACCACCGTCGGTGCGCAAGGCCg                                                                                                                                                    <  2:267749/149‑1 (MQ=255)
gCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCg                                                                                                                                                    >  1:142683/1‑149 (MQ=255)
      aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTAcgc                                                                                                                                             >  5:351103/1‑149 (MQ=255)
      aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTAcgc                                                                                                                                             >  6:151352/1‑149 (MQ=255)
      aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTAcgc                                                                                                                                             <  3:193174/149‑1 (MQ=255)
      aaTGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCCCACCACCGTCGGTGCGCAAGGCCGGTAcgc                                                                                                                                             >  6:361577/1‑149 (MQ=255)
      aaTGATGGCCGCCACCCCCCGAACCCTGACCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTAcgc                                                                                                                                             <  3:4680/149‑1 (MQ=255)
               cgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCg                                                                                                                                    <  3:347622/149‑1 (MQ=255)
               cgccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCg                                                                                                                                    <  2:142683/149‑1 (MQ=255)
                 ctcccccccGACCCCTGGCCCGTGCCGCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACGGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCggg                                                                                                                                  <  1:318519/146‑1 (MQ=255)
                 ccACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTACCGGCCTTGCgca                                                                                                                                  <  5:362709/149‑3 (MQ=255)
                  cACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGAGCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTACCGGCCTTGCgcac                                                                                                                                 <  5:66586/149‑4 (MQ=255)
                   aCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGtt                                                                                                                                <  7:27638/149‑1 (MQ=255)
                          gAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGtt                                                                                                                                <  4:255645/142‑1 (MQ=255)
                          gAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGtt                                                                                                                                >  3:255645/1‑142 (MQ=255)
                            aCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTAc                                                                                                                       <  5:129273/149‑1 (MQ=255)
                              ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                      >  6:349313/1‑148 (MQ=255)
                              ccTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACa                                                                                                                      >  2:360251/1‑148 (MQ=255)
                                  gCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAg                                                                                                                   <  5:30880/147‑1 (MQ=255)
                                        gCCTCGGCGGTGGCCGGCAGCAACGCAGGCAGTGGATAAACCTCTTCCAGGTACTCGATGCTCGCCCGCGCCTGCACCGGCCCCTCACCCCCGTCGGTGCGCCAGGCCGGTACGCGGCCCGGCGGGGTGACGGCGCCAACGCCCTCGcc                                                                                                           >  8:150721/1‑149 (MQ=255)
                                          cTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCAc                                                                                                         >  1:113944/1‑149 (MQ=255)
                                           tCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                        <  4:180591/133‑1 (MQ=255)
                                           tCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTa                                                                                                                        >  3:180591/1‑133 (MQ=255)
                                                           gcaACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTgcagc                                                                                        <  5:349313/149‑1 (MQ=255)
                                                              aCGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                     >  2:351867/1‑149 (MQ=255)
                                                                gCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACg                                                                                                        <  7:318904/128‑1 (MQ=255)
                                                                gCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACg                                                                                                        >  8:318904/1‑128 (MQ=255)
                                                                    gCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTg                                                                                            <  5:108823/136‑1 (MQ=255)
                                                                    gCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTg                                                                                            >  6:108823/1‑136 (MQ=255)
                                                                      tgtgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGa                                                                                                                              <  1:731/100‑1 (MQ=255)
                                                                      tgtgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGa                                                                                                                              >  2:731/1‑100 (MQ=255)
                                                                        tgGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGc                                                                            <  1:351867/148‑1 (MQ=255)
                                                                          gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACgg                                                                                                                           <  8:334072/99‑1 (MQ=255)
                                                                          gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACgg                                                                                                                           >  7:334072/1‑99 (MQ=255)
                                                                          gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                          <  5:68631/148‑1 (MQ=255)
                                                                          gATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                          <  6:351103/148‑1 (MQ=255)
                                                                             aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                          >  3:289171/1‑145 (MQ=255)
                                                                             aaaCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAg                                                                          <  4:289171/145‑1 (MQ=255)
                                                                               aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                     <  5:367877/132‑1 (MQ=255)
                                                                               aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAgg                                                                                     >  6:367877/1‑132 (MQ=255)
                                                                               aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                    <  4:145580/149‑1 (MQ=255)
                                                                               aCCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTg                                                                    <  5:62503/149‑1 (MQ=255)
                                                                                    ttCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAag                                                               >  4:196318/1‑149 (MQ=255)
                                                                                      ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                <  6:278199/146‑1 (MQ=255)
                                                                                      ccAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCCCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTaa                                                                <  5:361577/146‑1 (MQ=255)
                                                                                         ggTACTCGATGAGCGCCGGCGACTGCACCCGCACCTCACCACCGGCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGGTGACGGCTACATAGCCCTCGCCCCGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagccagc                                                          <  7:150721/149‑1 (MQ=255)
                                                                                          gTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCgg                                                                             <  5:22216/129‑1 (MQ=255)
                                                                                          gTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCgg                                                                             >  6:22216/1‑129 (MQ=255)
                                                                                             cTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCt                                                      >  1:136787/1‑149 (MQ=255)
                                                                                              tCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCtt                                                     >  1:161728/1‑149 (MQ=255)
                                                                                                             gACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAAt                                      <  2:113944/149‑1 (MQ=255)
                                                                                                              aCTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATg                                     >  2:221079/1‑149 (MQ=255)
                                                                                                               cTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                              <  4:232429/123‑1 (MQ=255)
                                                                                                               cTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAagc                                                              >  3:232429/1‑123 (MQ=255)
                                                                                                                         aCCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTa                          <  2:136787/149‑1 (MQ=255)
                                                                                                                         aCCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTa                          <  3:196318/149‑1 (MQ=255)
                                                                                                                            tCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGa                       >  3:34108/1‑149 (MQ=255)
                                                                                                                              accaccGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGAc                      <  2:161728/148‑1 (MQ=255)
                                                                                                                                caccGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAg                   >  1:173195/1‑149 (MQ=255)
                                                                                                                                       ggTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGaa             <  5:151352/148‑1 (MQ=255)
                                                                                                                                         tGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAg                         >  3:107540/1‑134 (MQ=255)
                                                                                                                                         tGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAg                         <  4:107540/134‑1 (MQ=255)
                                                                                                                                          ccgcAGGGCCGCGTACCGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACgg          >  1:310489/2‑148 (MQ=255)
                                                                                                                                           cgcAGGGCCGCGTACCGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTg                >  2:249721/1‑141 (MQ=255)
                                                                                                                                           cgcAGGGCCGCGTACCGGCCTTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTg                <  1:249721/141‑1 (MQ=255)
                                                                                                                                              aaGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTaa       >  7:225361/1‑147 (MQ=255)
                                                                                                                                               aGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAg    <  4:34108/149‑1 (MQ=255)
                                                                                                                                               aGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAg    <  1:360251/149‑1 (MQ=255)
                                                                                                                                                ggCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg   >  5:37346/1‑149 (MQ=255)
                                                                                                                                                ggCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAgg   >  8:154742/1‑149 (MQ=255)
                                                                                                                                                 gCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt  <  6:37346/149‑1 (MQ=255)
                                                                                                                                                 gCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGgt  <  2:173195/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGGGTTGACGGCTACATAGCCCTCGCCACGCTGCTCACCCTGCAGCAGGTTGACCGGCAGGGACTGGTAAGCCAGCCCCTTCAGTGCCAGGGCAATGCGCACCCGGTAGGACGAGGTGGAACGGTAATAGGT  >  NC_002947/5241566‑5241859

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: