Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,680,915 T→C 7.8% E123G (GAG→GGG) ‡ 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,9150TC 7.8% 183.5 / 5.1 64E123G (GAG→GGG) ‡PP_4986hemolysin III family channel protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base C (2/3);  total (29/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 = 8.78e-01

GAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACA  >  NC_002947/5680767‑5681056
                                                                                                                                                    |                                                                                                                                             
gaagaaAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCt                                                                                                                                               <  3:130677/149‑1 (MQ=255)
 aagaaAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTc                                                                                                                                              >  5:299473/1‑149 (MQ=255)
      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                                                                                                                      >  8:200221/1‑127 (MQ=255)
                                               agGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCaa                                                                                                                      <  7:200221/127‑1 (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                                    <  2:361886/149‑1 (MQ=255)
                                                                                                           aGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGc                                    <  5:49947/149‑1 (MQ=255)
                                                                                                           aGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGc                                    >  3:103751/1‑149 (MQ=255)
                                                                                                               cATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCgtg                                 >  1:298549/1‑148 (MQ=255)
                                                                                                                     cGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCa                                               >  4:123443/1‑128 (MQ=255)
                                                                                                                     cGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCa                                               <  3:123443/128‑1 (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                                                                 >  3:199117‑M1/8‑93 (MQ=255)
                                                                                                                                      gatcaagCCGCGCTCCGAGGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAg                                                                 <  4:199117‑M1/86‑1 (MQ=255)
                                                                                                                                         ggATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATc      <  3:220932/149‑1 (MQ=255)
                                                                                                                                         ggATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATACCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATc      <  5:349113/149‑1 (MQ=255)
                                                                                                                                          gATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATc      >  3:315693/1‑148 (MQ=255)
                                                                                                                                            tCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCa                     >  5:322383/1‑131 (MQ=255)
                                                                                                                                            tCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCa                     <  6:322383/131‑1 (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)
                                                                                                                                                    |                                                                                                                                             
GAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGCGCGGCTTGATCTCTTGCAACATGCCGATGATCGCCAGCCCCCAGACCACGCCGAACAGGCTCCATCCCCAAGGGCCGCGCAGGCTGACCAGGCAAAACGGCGTGTAGCTGCCGGCAATCAGCAGATAGATCGACA  >  NC_002947/5680767‑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: ‑