Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,253,887 G→C 12.5% F249L (TTC→TTG PP_4630 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,253,8870GC12.5% 172.5 / 13.5 63F249L (TTC→TTGPP_4630MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base C (4/4);  total (30/33)
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.65e-01

GCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGAT  >  NC_002947/5253741‑5254033
                                                                                                                                                  |                                                                                                                                                  
gCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGa                                                                                                                                                   >  4:118378/1‑148 (MQ=255)
  cAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGt                                                                                                                                                >  8:127127/1‑149 (MQ=255)
           ccGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGa                                                                                                                                       <  7:127127/149‑1 (MQ=255)
           ccGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGa                                                                                                                                       <  6:338761/149‑1 (MQ=255)
           ccGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGa                                                                                                                                       <  3:159507/149‑1 (MQ=255)
            cGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAg                                                                                                                                      >  2:318251/1‑149 (MQ=255)
             gCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGt                                                                                                                                     >  1:164253/1‑149 (MQ=255)
              cAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGtt                                                                                                                                    <  2:288148/149‑1 (MQ=255)
                   gcTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTctgct                                                                                                                               <  4:267311/149‑1 (MQ=255)
                      ggCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCtgtg                                                                                                                            >  4:33305/1‑149 (MQ=255)
                         gcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGc                                                                                                                         <  8:76300/149‑1 (MQ=255)
                                    cccTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATgg                                                                                                                <  1:318251/147‑1 (MQ=255)
                                     ccTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTTAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGTGaa                                                                                                             <  1:196304/149‑1 (MQ=255)
                                         ggcgcCGCGGCAGGCTACAGGGCTCCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCACTAACACCGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAAc                                                                                                          <  4:23432/147‑1 (MQ=255)
                                         cgcgcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACa                                                                                                         <  2:164253/149‑1 (MQ=255)
                                            gcCCAGTCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCa                                                                                                      <  1:202843/149‑1 (MQ=255)
                                            gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCa                                                                                                      >  5:289908/1‑149 (MQ=255)
                                            gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCa                                                                                                      >  3:143587/1‑149 (MQ=255)
                                               caggcaggCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCa                                                                                                   <  7:242494/149‑1 (MQ=255)
                                                       cTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTcc                                                                                           <  6:289908/149‑1 (MQ=255)
                                                          cAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGc                                                                                        >  4:34331/1‑149 (MQ=255)
                                                           aGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCACTTCAAGCGCCTGGc                                                                                                                                         <  8:274570/99‑1 (MQ=255)
                                                           aGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCACTTCAAGCGCCTGGc                                                                                                                                         >  7:274570/1‑99 (MQ=255)
                                                             ggCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTc                                                                                     <  3:34331/149‑1 (MQ=255)
                                                                tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGc                                                                                  >  5:149197/1‑149 (MQ=255)
                                                                    agccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCACCCAGTAACAACGGCCATTCCAGCATCGACGCCAGGCGCTTCAAGTGCCTGGCCACTACTGCTGACGGCGACTGAAGGGAAAACATCACCAGGGCGTCCGGCTTCATCCGt                                                                               >  4:268491/1‑148 (MQ=255)
                                                                     gccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCAGCTTCAGCCGTTc                                                                             <  1:59631/149‑1 (MQ=255)
                                                                      ccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCg                                                                            >  1:250827/1‑149 (MQ=255)
                                                                         ggCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAATAACAACGGGCATTCCAGCATCGACGCCAGGCCCTTGAAGTGCCTGGCGAGTTCTGCTGTggg                                                                                                                          >  7:235359/1‑100 (MQ=255)
                                                                               ccTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGTACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                                                   >  6:214429/1‑149 (MQ=255)
                                                                                cTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGt                                                                               >  6:187335/1‑136 (MQ=255)
                                                                                cTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGt                                                                               <  5:187335/136‑1 (MQ=255)
                                                                                   tACCGGTTCCGCTGCTTCACCTGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgcgc                                                                 <  3:118378/147‑1 (MQ=255)
                                                                                     ccGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGc                                                             >  7:205559/1‑149 (MQ=255)
                                                                                     ccGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTACAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGATCTGCTGAGGGCGACTGAAGGGAAAACACCACCAAGCCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAg                                                              >  5:50254/1‑148 (MQ=255)
                                                                                             gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAg                                                     >  2:163282/1‑149 (MQ=255)
                                                                                                gctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCgg                                                  >  4:194883/1‑149 (MQ=255)
                                                                                                     aCCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGcc                                             >  4:316964/1‑149 (MQ=255)
                                                                                                     aCCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGcc                                             >  2:110097/1‑149 (MQ=255)
                                                                                                     aCCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGc                                              <  3:33305/148‑1 (MQ=255)
                                                                                                        agccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGt                                          <  2:250827/149‑1 (MQ=255)
                                                                                                            agTAACAAGGGGCATTCCAGCATCGTCGGATGGGGCTTGATGTGGCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACTGCACGAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCTAGCTCGTGCAGCTGCTGCCCGATctc                                      <  6:50254/149‑3 (MQ=255)
                                                                                                               aacaacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCCCCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAg                                   >  2:258108/1‑149 (MQ=255)
                                                                                                               aacaacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGc                                                    <  5:250566/132‑1 (MQ=255)
                                                                                                               aacaacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGc                                                    >  6:250566/1‑132 (MQ=255)
                                                                                                                 caacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGAt                                                                                                                  >  8:158018/1‑68 (MQ=255)
                                                                                                                 caacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGAt                                                                                                                  <  7:158018/68‑1 (MQ=255)
                                                                                                                      ggCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATc                            <  1:163282/149‑1 (MQ=255)
                                                                                                                             cAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCgt                      <  3:284772/148‑1 (MQ=255)
                                                                                                                                    gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtc              <  3:316964/149‑1 (MQ=255)
                                                                                                                                    gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtc              >  5:315885/1‑149 (MQ=255)
                                                                                                                                       gCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCt           <  6:83323/149‑1 (MQ=255)
                                                                                                                                       gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCtt                                                                                      <  1:147647/74‑1 (MQ=38)
                                                                                                                                       gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCtt                                                                                      >  2:147647/1‑74 (MQ=38)
                                                                                                                                       gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCgtgt                    <  2:117241/140‑1 (MQ=255)
                                                                                                                                       gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCgtgt                    >  1:117241/1‑140 (MQ=255)
                                                                                                                                            gcacTTCAAGCGCCTGGCGAGTTCTGCTGTGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAg     <  7:74982/146‑1 (MQ=255)
                                                                                                                                             cgcTTGAAGTGCCGGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAg     <  6:213306/149‑1 (MQ=255)
                                                                                                                                                ttGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGAt  >  2:245792/1‑149 (MQ=255)
                                                                                                                                                ttGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGACGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtcg             <  1:118668/138‑1 (MQ=255)
                                                                                                                                                ttGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGACGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtcg             >  2:118668/1‑138 (MQ=255)
                                                                                                                                                 tGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtc              >  6:62293/1‑136 (MQ=255)
                                                                                                                                                 tGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtc              <  5:62293/136‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGAT  >  NC_002947/5253741‑5254033

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: