Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,872,021 T→A 12.7% D211V (GAC→GTC)  ptsP ← phosphoenolpyruvate‑dependent regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,872,0210TA12.7% 162.4 / 14.3 63D211V (GAC→GTC) ptsPphosphoenolpyruvate‑dependent regulator
Reads supporting (aligned to +/- strand):  ref base T (32/22);  new base A (5/3);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATCGA  >  NC_002947/5871890‑5872169
                                                                                                                                   |                                                                                                                                                     
tcatcaGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTcc                                                                                                                                     <  4:115048/149‑1 (MQ=255)
  atcaGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCgg                                                                                                                                   >  4:185318/1‑149 (MQ=255)
     aGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 <  3:185318/148‑1 (MQ=255)
     aGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 <  2:187789/148‑1 (MQ=255)
      gATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCgg                                                                                                                                   <  4:279550/145‑1 (MQ=255)
      gATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCgg                                                                                                                                   >  3:279550/1‑145 (MQ=255)
         aGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACt                                                                                                                            >  7:240337/1‑149 (MQ=255)
             gTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 >  1:187604/1‑140 (MQ=255)
             gTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGa                                                                                                                                 <  2:187604/140‑1 (MQ=255)
                 aaCAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTc                                                                                                                    >  6:348164/1‑149 (MQ=255)
                                       gcAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCAGGAGGCAACATCaccacc                                                                                              <  8:273558/149‑1 (MQ=255)
                                         aGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGAACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCg                                                                                             >  3:242800/1‑148 (MQ=255)
                                              ggTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCAt                                                                                                                                                   <  4:103551/89‑1 (MQ=255)
                                              ggTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCAt                                                                                                                                                   >  3:103551/1‑89 (MQ=255)
                                              ggTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGcc                                                                                       >  8:317680/1‑149 (MQ=255)
                                              ggTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGc                                                                                        <  4:402851/148‑1 (MQ=255)
                                                   ccAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCTGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCAcc                                                                                  >  3:165143/1‑149 (MQ=255)
                                                      gTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACcgcg                                                                               <  4:156597/149‑1 (MQ=255)
                                                            ccGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCAc                                                                           >  7:124983/1‑147 (MQ=255)
                                                                   cTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAc                                                                  >  2:9755/1‑149 (MQ=255)
                                                                     tGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAcc                                                                 >  8:35060/1‑148 (MQ=255)
                                                                     tGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAcc                                                                 <  7:35060/148‑1 (MQ=255)
                                                                      gcgcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCg                                                               >  3:366609/1‑149 (MQ=255)
                                                                       cgcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCgg                                                              >  5:338903/1‑149 (MQ=255)
                                                                        gcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCTGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGc                                                             <  4:165143/149‑1 (MQ=255)
                                                                        gcATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGc                                                             >  4:229979/1‑149 (MQ=255)
                                                                         cATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGc                                                             >  1:117507/1‑149 (MQ=255)
                                                                         cATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAg                                                                        >  5:406721/1‑137 (MQ=255)
                                                                         cATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAg                                                                        <  6:406721/137‑1 (MQ=255)
                                                                           tGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGcaca                                                          <  7:374206/149‑1 (MQ=255)
                                                                            gtcgtcGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTcc                                                                                                                         >  7:176876/1‑85 (MQ=255)
                                                                            gtcgtcGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTcc                                                                                                                         <  8:176876/85‑1 (MQ=255)
                                                                                tcgCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCAc                                                     >  6:60536/1‑149 (MQ=255)
                                                                                   cgcACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGaa                                                  <  8:240337/149‑1 (MQ=255)
                                                                                      acacCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACAcc                                                        >  5:421402/1‑140 (MQ=255)
                                                                                      acacCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACAcc                                                        <  6:421402/140‑1 (MQ=255)
                                                                                          cTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAc                                                                  >  7:289271/1‑126 (MQ=255)
                                                                                          cTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAAc                                                                  <  8:289271/126‑1 (MQ=255)
                                                                                               aGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTg                                      <  3:229979/149‑1 (MQ=255)
                                                                                                 ggCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGaa                                    >  8:93476/1‑149 (MQ=255)
                                                                                                 ggCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGGCATCGACGGTCTTGTCCGGCACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGaa                                    <  4:242800/149‑1 (MQ=255)
                                                                                                         ggAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCCCCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTg                            >  4:53582/1‑149 (MQ=255)
                                                                                                         ggAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCa                                                      >  3:114348/1‑123 (MQ=255)
                                                                                                         ggAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCa                                                      <  4:114348/123‑1 (MQ=255)
                                                                                                          gAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACAACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGc                           >  1:130826/1‑149 (MQ=255)
                                                                                                           aaCAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCCGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGcc                          >  6:75115/1‑149 (MQ=255)
                                                                                                                cTTGAGTTCGGCGTCGATGTCATCGAGGGTGTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACGACCGCACGCCCCACGGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGGGCTTCCTGAATACCTTTGCCCTGa                      <  1:303699/148‑1 (MQ=255)
                                                                                                                cTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACGACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCCCGAAGCACGCTTCCTGAATACCTTTGCCCTGa                      >  2:303699/1‑148 (MQ=255)
                                                                                                                 ttGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACg                    >  5:77747/1‑149 (MQ=255)
                                                                                                                 ttGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACg                    >  5:289203/1‑149 (MQ=255)
                                                                                                                     gTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGccc                <  7:93476/149‑1 (MQ=255)
                                                                                                                       agaccGTCTATGACATCGACTGTCTTGTCCGGGACCACTTACAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGc                           <  8:316118/132‑1 (MQ=255)
                                                                                                                       agaccGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGACCCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGc                           >  7:316118/5‑136 (MQ=255)
                                                                                                                       agaccGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGc                           <  8:310394/132‑1 (MQ=255)
                                                                                                                       agaccGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGc                           >  7:310394/5‑136 (MQ=255)
                                                                                                                        cGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAgg             >  4:353869/1‑149 (MQ=255)
                                                                                                                         ggCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGc            >  1:60262/1‑149 (MQ=255)
                                                                                                                         accGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCttt                             >  7:388975/3‑132 (MQ=255)
                                                                                                                         accGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCttt                             <  8:388975/130‑1 (MQ=255)
                                                                                                                          ccGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGcc           >  2:234521/2‑149 (MQ=255)
                                                                                                                           cGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCAAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGCACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCACGCCa          >  7:290630/1‑149 (MQ=255)
                                                                                                                            gTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCAc         <  5:60536/149‑1 (MQ=255)
                                                                                                                                 tGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATc    >  1:54868/1‑149 (MQ=255)
                                                                                                                                   tCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATCGa  <  8:124983/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                     
TCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACGGATCGA  >  NC_002947/5871890‑5872169

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