Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,948,117 T→A 6.0% H513L (CAC→CTC)  PP_1749 ← acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,948,1170TA 6.0% 189.1 / 2.9 67H513L (CAC→CTC) PP_1749acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base A (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

GTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGG  >  NC_002947/1947971‑1948262
                                                                                                                                                   |                                                                                                                                                 
gTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGTACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTg                                                                                                                                                  >  5:200333/1‑148 (MQ=255)
   cATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                                               <  7:21398/148‑1 (MQ=255)
   cATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTACCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                                               <  5:364359/148‑1 (MQ=255)
    aTGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCTGCACCCGGTGCTGGc                                                                                                                                             <  7:351856/149‑1 (MQ=255)
     tGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                                               <  4:337977/146‑1 (MQ=255)
     tGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                                               >  3:337977/1‑146 (MQ=255)
       gTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGtt                                                                                                                                          <  8:66936/149‑1 (MQ=255)
        ttGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCTGCACCCGGTGCTGGccgat                                                                                                                                         <  1:204720/149‑5 (MQ=255)
         tGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCg                                                                                                                                        >  5:428480/1‑149 (MQ=255)
            ccaggccagCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTg                                                                                                                                     <  5:290651/149‑1 (MQ=255)
              aggccagCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGAc                                                                                                                                   <  1:114514/149‑1 (MQ=255)
              aggccagCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGAc                                                                                                                                   <  3:50266/149‑1 (MQ=255)
                 ccagCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGtc                                                                                                                                <  2:387636/149‑1 (MQ=255)
                                gCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTgccgc                                                                                                                  <  6:50666/148‑1 (MQ=255)
                                 ccTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGccgccg                                                                                                                >  8:189723/1‑149 (MQ=255)
                                  cTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCgg                                                                                                               >  4:190200/1‑149 (MQ=255)
                                     gatgatCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGgtgt                                                                                                            <  5:106299/149‑1 (MQ=255)
                                     gatgatCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGgt                                                                                                              <  4:164752/147‑1 (MQ=255)
                                         atCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCa                                                                                                         <  2:147671/148‑1 (MQ=255)
                                          tCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAgg                                                                                                                         <  2:3781/131‑1 (MQ=255)
                                          tCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAgg                                                                                                                         >  1:3781/1‑131 (MQ=255)
                                            aCGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTTCCGCGCGCACGGCAGCATCGGCCAGCACCGGGAGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGtt                                                                                                     >  6:209093/1‑149 (MQ=255)
                                            aCGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGt                                                                                                      <  2:250397/148‑1 (MQ=255)
                                            aCGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGAGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGtt                                                                                                     >  6:202758/1‑149 (MQ=255)
                                             cGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGGTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTg                                                                                                    >  2:309569/1‑149 (MQ=255)
                                                  tCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGt                                                                                               >  1:55783/1‑149 (MQ=255)
                                                    cGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACGGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgc                                                                                             <  3:120800/149‑1 (MQ=255)
                                                    cGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTg                                                                                              <  4:164087/148‑1 (MQ=255)
                                                       cTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgcc                                                                                          <  2:194213/149‑1 (MQ=255)
                                                           ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCa                                                                                       <  8:274417/148‑1 (MQ=255)
                                                           ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCa                                                                                       >  7:274417/1‑148 (MQ=255)
                                                           ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCAc                                                                                      >  1:350021/1‑149 (MQ=255)
                                                            cGGCATCACGCACCATGCAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGGGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                     <  3:39234/149‑1 (MQ=255)
                                                            cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                     <  7:193786/149‑1 (MQ=255)
                                                               cATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGcc                                                                                  <  4:123012/149‑1 (MQ=255)
                                                                  cACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGCCaac                                                                               <  3:190200/149‑3 (MQ=255)
                                                                  cACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAgg                                                                               >  5:116682/1‑149 (MQ=255)
                                                                   aCGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCCGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGgc                                                                              >  5:422771/1‑149 (MQ=255)
                                                                       aCCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCGGTGCGCCGCACGGCCAGGCGCTg                                                                          <  2:430051/149‑1 (MQ=255)
                                                                        ccATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTgg                                                                         <  6:428480/149‑1 (MQ=255)
                                                                        ccATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTgg                                                                         <  1:309569/149‑1 (MQ=255)
                                                                               gTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGgccggccg                                                                  >  3:153363/1‑149 (MQ=255)
                                                                                  cAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCcggccggc                                                                <  6:422771/148‑1 (MQ=255)
                                                                                   aGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGAGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                              <  5:209093/149‑1 (MQ=255)
                                                                                   aGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGAGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                              <  5:202758/149‑1 (MQ=255)
                                                                                    gCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                              >  4:207130/1‑148 (MQ=255)
                                                                                    gCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAac                                                             <  6:200333/149‑1 (MQ=255)
                                                                                          ccaccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtcg                                                       <  3:287247/149‑1 (MQ=255)
                                                                                          ccaccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtcg                                                       >  8:262533/1‑149 (MQ=255)
                                                                                            accGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtcgt                                                      >  6:219048/1‑148 (MQ=255)
                                                                                              cGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCa                                                    >  5:420470/1‑148 (MQ=255)
                                                                                               ggAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATa                                                  <  3:207130/149‑1 (MQ=255)
                                                                                                gAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAgg                                                                               <  6:17243/119‑1 (MQ=255)
                                                                                                      ccAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACg                                           >  4:313516/1‑149 (MQ=255)
                                                                                                        aGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCATGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGcc                                         >  6:229154/1‑149 (MQ=255)
                                                                                                           tgccGTGCGCGCGGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgc                                                                                             >  3:360342/3‑94 (MQ=255)
                                                                                                           tgccGTGCGCGCGGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgc                                                                                             <  4:360342/92‑1 (MQ=255)
                                                                                                               gTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgcg                                  >  8:334806/1‑149 (MQ=255)
                                                                                                                tGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgcgc                                 <  5:333902/149‑1 (MQ=255)
                                                                                                                               cagcaTCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCg                   <  3:313516/149‑1 (MQ=255)
                                                                                                                                    tCGGCCAGCACCGGGGGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGGCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtc              <  5:36374/149‑1 (MQ=255)
                                                                                                                                       gCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtc           >  6:371892/1‑149 (MQ=255)
                                                                                                                                       gCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtc           >  6:331545/1‑149 (MQ=255)
                                                                                                                                            cACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCgtcgtc        >  5:185799/1‑147 (MQ=255)
                                                                                                                                            cACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCCTGCGCCGAACGGACAGGCGCAGGCCGGCCGGCAACACGTAGTCATAGCCAACGCCCGCTGCCCGCAACGTGCCCTCGGTTACgtcgtcaa      >  7:384895/1‑147 (MQ=255)
                                                                                                                                              ccGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGc    >  1:291681/1‑149 (MQ=255)
                                                                                                                                                 ggTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCgg  >  2:149662/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTACTCCGGCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGG  >  NC_002947/1947971‑1948262

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: