Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,926,513 T→G 7.1% D118A (GAC→GCC)  PP_1726 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,926,5130TG 7.1% 203.7 / 2.6 66D118A (GAC→GCC) PP_1726ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base G (3/2);  total (35/31)
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

TTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGGTGT  >  NC_002947/1926375‑1926660
                                                                                                                                          |                                                                                                                                                   
ttGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgct                                                                                                                                           >  5:63876/1‑149 (MQ=255)
ttGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgct                                                                                                                                           >  3:11658/1‑149 (MQ=255)
ttGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgct                                                                                                                                           >  3:115107/1‑149 (MQ=255)
    cTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggc                                                                                                                                       >  8:60421/1‑149 (MQ=255)
    cTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggc                                                                                                                                       >  5:217400/1‑149 (MQ=255)
     tttCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccg                                                                                                                                          >  4:277564/1‑143 (MQ=255)
     tttCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccg                                                                                                                                          <  3:277564/145‑3 (MQ=255)
         tcCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGt                                                                                                                                   >  1:58294/1‑148 (MQ=255)
           cAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGt                                                                                                                                   <  7:60421/146‑1 (MQ=255)
           cAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTaa                                                                                                                                 <  6:63876/148‑1 (MQ=255)
           cAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTaa                                                                                                                                 <  2:58294/148‑1 (MQ=255)
            aGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAgg                                                                                                                               >  1:65352/1‑149 (MQ=255)
              gTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAgg                                                                                                                                               <  4:312373/131‑1 (MQ=255)
              gTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAgg                                                                                                                                               >  3:312373/1‑131 (MQ=255)
                    ccAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGa                                                                                                                       <  1:136612/149‑1 (MQ=255)
                     cAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctg                                                                                                                            >  2:283696/1‑143 (MQ=255)
                     cAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctg                                                                                                                            <  1:283696/143‑1 (MQ=255)
                              ggCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTg                                                                                                                                             <  6:101656/117‑1 (MQ=255)
                              ggCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTg                                                                                                                                             >  5:101656/1‑117 (MQ=255)
                               gCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGt                                                                                                             >  6:73523/1‑148 (MQ=255)
                                      ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGa                                                                                                     <  4:115107/149‑1 (MQ=255)
                                      ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGa                                                                                                     <  4:11658/149‑1 (MQ=255)
                                      ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGAGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTCGCCCATTCCGGCACCTGGTCCCAGGTGCTTCGCTTGTAAGCCTGGCTGACCCCCTTGGTCACAGCGa                                                                                                     >  6:106113/1‑149 (MQ=255)
                                       tctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccg                                                                                                                                          <  6:194471/111‑3 (MQ=255)
                                       tctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccg                                                                                                                                          >  5:194471/1‑109 (MQ=255)
                                         ttCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               <  8:181227/136‑1 (MQ=255)
                                         ttCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               >  7:181227/1‑136 (MQ=255)
                                            cTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGc                                                                                               >  4:36763/1‑149 (MQ=255)
                                               gACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGcc                                                                                              <  6:217400/147‑1 (MQ=255)
                                                 ccAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAg                                                                                          >  1:13701/1‑149 (MQ=255)
                                                     gTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggc                                                                                                                              <  2:227577/109‑1 (MQ=255)
                                                     gTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggc                                                                                                                              >  1:227577/1‑109 (MQ=255)
                                                           gttgATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCgg                                                                                                 <  2:286715/132‑1 (MQ=255)
                                                           gttgATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCgg                                                                                                 >  1:286715/1‑132 (MQ=255)
                                                              gatgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGt                                                                             >  3:23333/1‑149 (MQ=255)
                                                               atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGt                                                                             <  3:36763/148‑1 (MQ=255)
                                                                        gCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTc                                                                   <  5:73523/149‑1 (MQ=255)
                                                                               tGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTg                                                            <  1:140349/149‑1 (MQ=255)
                                                                                gCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTgg                                                           >  4:197114/1‑149 (MQ=255)
                                                                                          gCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCtt                                                 >  5:98102/1‑149 (MQ=255)
                                                                                            cAGCTCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTtct                                               <  5:106113/149‑1 (MQ=255)
                                                                                             aGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCGTGTAAGGCTGGCTGACGCCCTTGGGCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTtctc                                              <  5:257063/149‑1 (MQ=255)
                                                                                              gcgcCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCg                                             >  4:56049/1‑149 (MQ=255)
                                                                                                   cAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTc                                        >  3:91916/1‑149 (MQ=255)
                                                                                                    aGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTc                                        >  3:312498/1‑148 (MQ=255)
                                                                                                     gTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGtt                                                       >  5:255956/1‑132 (MQ=255)
                                                                                                     gTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGtt                                                       <  6:255956/132‑1 (MQ=255)
                                                                                                         ccGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAAc                                    >  8:309939/1‑147 (MQ=255)
                                                                                                         ccGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAAc                                    <  7:309939/147‑1 (MQ=255)
                                                                                                          cGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTg                                 <  2:13701/149‑1 (MQ=255)
                                                                                                            tccttgtccttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGc                               <  3:197114/149‑1 (MQ=255)
                                                                                                                  tccttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTc                         <  4:23333/149‑1 (MQ=255)
                                                                                                                    cttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCt                       >  7:172607/1‑149 (MQ=255)
                                                                                                                         cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGcc                  >  6:74781/1‑149 (MQ=255)
                                                                                                                          ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTgg                     >  1:16109/1‑145 (MQ=255)
                                                                                                                          ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTgg                     <  2:16109/145‑1 (MQ=255)
                                                                                                                            aTTCCGGCACCTGGTCCCAGGTGCTTGGCTTCTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAg               >  5:338314/1‑149 (MQ=255)
                                                                                                                               caagCACCTGGGACCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCt             >  7:308256/4‑148 (MQ=255)
                                                                                                                                cGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCTGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTc            <  1:217191/148‑1 (MQ=255)
                                                                                                                                 ggCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCAt          >  4:293318/1‑149 (MQ=255)
                                                                                                                                   cACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGt        <  3:56049/149‑1 (MQ=255)
                                                                                                                                     ccTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTc                                        >  6:196016/1‑115 (MQ=255)
                                                                                                                                     ccTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTTATGTCg      <  8:172607/149‑1 (MQ=255)
                                                                                                                                         gTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCAt          <  7:314871/141‑1 (MQ=255)
                                                                                                                                         gTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCAt          >  8:314871/1‑141 (MQ=255)
                                                                                                                                         gTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGgtgt  <  2:65352/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGGTGT  >  NC_002947/1926375‑1926660

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