Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,680,919 A→G 7.9% S122P (TCC→CCC) ‡ PP_4986 ← hemolysin III family channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,680,9190AG 7.9% 178.0 / 4.6 63S122P (TCC→CCC) ‡PP_4986hemolysin III family channel protein
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base G (2/3);  total (28/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 = 1.00e+00

AATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACA  >  NC_002947/5680773‑5681056
                                                                                                                                                  |                                                                                                                                         
aaTGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAgc                                                                                                                                         >  3:133333/1‑149 (MQ=255)
 atgatgCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAgcg                                                                                                                                        >  3:217578/1‑149 (MQ=255)
   gatgCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGcgcg                                                                                                                                      >  5:370128/1‑149 (MQ=255)
    atgCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGcgcg                                                                                                                                      <  4:133333/148‑1 (MQ=255)
     tgCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGc                                                                                                                                    <  6:195496/149‑1 (MQ=255)
      gCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCt                                                                                                                                   >  3:37664/1‑149 (MQ=255)
       ccGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCtt                                                                                                                                  >  8:92687/1‑149 (MQ=255)
         gACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCGAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCTCCGAGGCGCggatcct                                                                                                                                <  8:239871/149‑6 (MQ=255)
              tgtAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATctct                                                                                                                           >  2:280261/1‑149 (MQ=255)
                tAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCtt                                                                                                                          <  5:45606/148‑1 (MQ=255)
                tAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCtt                                                                                                                          <  7:85824/148‑1 (MQ=255)
                  gAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGc                                                                                                                        <  1:13641/148‑1 (MQ=255)
                    acacGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                      >  5:187440/1‑148 (MQ=255)
                    acacGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                      >  4:17625/1‑148 (MQ=255)
                        gccgccGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCCGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGc                                                                                                                 >  2:320392/1‑149 (MQ=255)
                        gccgccGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTg                                                                                                                         >  2:240809/1‑141 (MQ=255)
                        gccgccGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTg                                                                                                                         <  1:240809/141‑1 (MQ=255)
                        gccgccGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGc                                                                                                                 <  8:237864/149‑1 (MQ=255)
                        gccgccGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGc                                                                                                                 <  4:217578/149‑1 (MQ=255)
                        gccgccGGCAGCCAGCCAGGCGAAGACGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCCCGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGc                                                                                                                 >  7:167932/1‑149 (MQ=255)
                           gccgGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCga                                                                                                              <  7:26970/149‑1 (MQ=255)
                                       ccagGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGcc                                                                                                  <  4:37664/149‑1 (MQ=255)
                                         agGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                      <  7:200221/127‑1 (MQ=255)
                                         agGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                      >  8:200221/1‑127 (MQ=255)
                                           gCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCa                                                                                              <  6:187440/149‑1 (MQ=255)
                                                  cggccgTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGTCCACg                                                                                       <  6:299473/149‑1 (MQ=255)
                                                          cccAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATCATCGCCAGCCCCCAGACCACGCCGAACa                                                                                <  1:208646/148‑1 (MQ=255)
                                                           ccAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAgg                                                                              <  7:92687/149‑1 (MQ=255)
                                                           ccAGCGAGTTCAGCAGCGGCTTGACTGCCAACAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACATGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAgg                                                                              <  8:167932/149‑1 (MQ=255)
                                                             aGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCGATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCt                                                                            <  5:61848/149‑1 (MQ=255)
                                                                  gTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATc                                                                       >  6:49947/1‑149 (MQ=255)
                                                                     cagcagCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATccc                                                                     >  4:220932/1‑148 (MQ=255)
                                                                         agcGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAgg                                                                <  3:17625/149‑1 (MQ=255)
                                                                          gcGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAggg                                                               <  1:280261/149‑1 (MQ=255)
                                                                           cGGCTTGACTGCCACCAGGACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGc                                                              <  5:221837/149‑1 (MQ=255)
                                                                               ttGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCcgcg                                                          >  6:2978/1‑149 (MQ=255)
                                                                                      ccaccaGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCt                                                    <  6:370128/148‑1 (MQ=255)
                                                                                          caGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCa                                               <  6:301355/149‑1 (MQ=255)
                                                                                              aCGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATACCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGc                                            >  6:349113/1‑148 (MQ=255)
                                                                                                gATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCa                                                                                                     >  2:318256/1‑89 (MQ=255)
                                                                                                gATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCa                                                                                                     <  1:318256/89‑1 (MQ=255)
                                                                                                     aGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGc                                    <  5:49947/149‑1 (MQ=255)
                                                                                                     aGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGc                                    <  2:361886/149‑1 (MQ=255)
                                                                                                     aGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGc                                    >  3:103751/1‑149 (MQ=255)
                                                                                                         cATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCgtg                                 >  1:298549/1‑148 (MQ=255)
                                                                                                               cGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCa                                               <  3:123443/128‑1 (MQ=255)
                                                                                                               cGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCa                                               >  4:123443/1‑128 (MQ=255)
                                                                                                                           caagagatcaagCCGCGCTCCGAGGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAgg                                                                >  1:137241‑M1/13‑99 (MQ=255)
                                                                                                                           caagagatcaagCCGCGCTCCGAGGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAgg                                                                <  2:137241‑M1/87‑1 (MQ=255)
                                                                                                                             tCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCaga            <  5:2978/149‑1 (MQ=255)
                                                                                                                              cGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCaga            >  3:80137/1‑148 (MQ=255)
                                                                                                                               gACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCagat           <  1:320392/148‑1 (MQ=255)
                                                                                                                               gACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCaga            >  5:249321/1‑147 (MQ=255)
                                                                                                                                gatcaagCCGCGCTCCGAGGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAg                                                                 <  4:199117‑M1/86‑1 (MQ=255)
                                                                                                                                gatcaagCCGCGCTCCGAGGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAg                                                                 >  3:199117‑M1/8‑93 (MQ=255)
                                                                                                                                   ggATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATc      <  3:220932/149‑1 (MQ=255)
                                                                                                                                   ggATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATACCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATc      <  5:349113/149‑1 (MQ=255)
                                                                                                                                    gATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATc      >  3:315693/1‑148 (MQ=255)
                                                                                                                                      tCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCa                     <  6:322383/131‑1 (MQ=255)
                                                                                                                                      tCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCa                     >  5:322383/1‑131 (MQ=255)
                                                                                                                                        CGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACa  >  4:315693‑M1/1‑148 (MQ=255)
                                                                                                                                         GCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATc                  <  7:128996‑M1/131‑1 (MQ=255)
                                                                                                                                         GCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATc                  >  8:128996‑M1/1‑131 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
AATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACA  >  NC_002947/5680773‑5681056

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