Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,984,695 C→A 7.5% G142V (GGC→GTC)  ycaC‑II ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,984,6950CA 7.5% 250.1 / 8.2 80G142V (GGC→GTC) ycaC‑IIputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (45/29);  new base A (4/2);  total (49/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 = 6.08e-01

CGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTC  >  NC_002947/5984548‑5984843
                                                                                                                                                   |                                                                                                                                                    
cGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGc                                                                                                                                                      <  2:55339/148‑1 (MQ=255)
 gCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGccc                                                                                                                                                    >  1:119596/1‑149 (MQ=255)
   tGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGACCGa                                                                                                                                                  >  6:379507/1‑149 (MQ=255)
        cagcagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGACCGATGcct                                                                                                                                             <  5:379507/149‑3 (MQ=255)
           cagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGTCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCgg                                                                                                                                          >  1:97775/1‑149 (MQ=255)
               ttACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcac                                                                                                                                      <  1:202743/149‑1 (MQ=255)
                 aCCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcacca                                                                                                                                    >  5:168577/1‑149 (MQ=255)
                  cccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcacca                                                                                                                                    <  1:12908/148‑1 (MQ=255)
                   ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                   >  8:77567/1‑148 (MQ=255)
                   ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                   >  5:122515/1‑148 (MQ=255)
                   ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                   >  7:65365/1‑148 (MQ=255)
                     aGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGCACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGaa                                                                                                                                >  5:269132/1‑149 (MQ=255)
                     aGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGaa                                                                                                                                >  4:269372/1‑149 (MQ=255)
                     aGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGACCGATGCCTCGGGcaccttcaa                                                                                                                                <  8:73422/149‑6 (MQ=37)
                               tGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCaaaa                                                                                                                      <  6:168577/149‑1 (MQ=255)
                                 tCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAAcc                                                                                                                    >  3:129457/1‑149 (MQ=255)
                                   gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCg                                                                                                                                           >  7:276282/1‑124 (MQ=255)
                                   gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCg                                                                                                                                           <  8:276282/124‑1 (MQ=255)
                                   gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCt                                                                                                                   <  5:322356/148‑1 (MQ=255)
                                        gcCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCGTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCAGCGGTCACCACGAACACTTCAACACCTTCCGc                                                                                                              >  2:98363/1‑148 (MQ=255)
                                           aGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAgc                                                                                                          >  3:109548/1‑149 (MQ=255)
                                             tCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAgcgc                                                                                                        >  2:249484/1‑149 (MQ=255)
                                             tCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAgcgc                                                                                                        >  2:206086/1‑149 (MQ=255)
                                                cgGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCg                                                                                                                                                   >  6:330322/1‑103 (MQ=255)
                                                cgGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCg                                                                                                                                                   <  5:330322/103‑1 (MQ=255)
                                                cgGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGa                                                                                                     <  2:97775/149‑1 (MQ=255)
                                                     cAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCACAACCTTCCGCCAGCGCCGACAggg                                                                                                >  7:205094/1‑149 (MQ=255)
                                                     cAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAgg                                                                                                 <  3:269372/148‑1 (MQ=255)
                                                       gTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTc                                                                                              >  5:64801/1‑149 (MQ=255)
                                                          cGCAGGTCACCGAGAACCATTTCATCATCTGCGCACCGGCGTGGGTCATGCGTACCCAAGCCGGTTGCTGCACGGTTTCGTTGAAGGTGCCCGTGGCATCGGTAAGCAGGAACTCTTCAAAACGTTCCGCCAGCGCCGACAGGGTCgg                                                                                            <  1:69412/148‑1 (MQ=255)
                                                           gCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCggg                                                                                           <  8:158965/148‑1 (MQ=255)
                                                             aGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAg                                                                                        <  6:64801/149‑1 (MQ=255)
                                                               gCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGCGAAGGc                                                                                      >  6:244292/1‑149 (MQ=255)
                                                                ccACCGAGAACCAGTTCATCATCTGCGCACCCGGCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGcc                                                                                     <  5:199470/149‑1 (MQ=255)
                                                                 cACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCa                                                                                    >  5:347161/1‑149 (MQ=255)
                                                                  aCCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCggg                                                                                           >  4:193724/1‑141 (MQ=255)
                                                                  aCCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCggg                                                                                           <  3:193724/141‑1 (MQ=255)
                                                                       gaACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGATTGCTGAACGGTTTCGTTGAAGTTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGcaca                                                                              >  5:381028/1‑149 (MQ=255)
                                                                          ccAGTTCATCATCTACGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCACACCCTTCCGCCAGCGCCGACATGGTCGGGAAGGCCACGccaccat                                                                           >  4:15069/1‑143 (MQ=38)
                                                                           cAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATc                                                                          >  6:273767/1‑149 (MQ=255)
                                                                           cAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATc                                                                          >  2:357050/1‑149 (MQ=255)
                                                                                 atcatcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAACACCTTCCGCCAGCGCCGACAGGGTCGGGAACGCCACGCACACATCGGTcac                                                                    >  5:365542/1‑149 (MQ=255)
                                                                                 atcatcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGGTTCGTTGAAGGTGCGCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTcac                                                                    <  6:240557/149‑1 (MQ=255)
                                                                                    atcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTcacca                                                                  >  5:271227/1‑148 (MQ=255)
                                                                                        gcgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAgg                                                                                       <  1:163939/123‑1 (MQ=255)
                                                                                        gcgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAgg                                                                                       >  2:163939/1‑123 (MQ=255)
                                                                                        gcgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCg                                                             >  5:204949/1‑149 (MQ=255)
                                                                                         cgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACAcc                                                              <  6:347161/147‑1 (MQ=255)
                                                                                              ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCgat                                                        >  3:137850/1‑148 (MQ=255)
                                                                                              ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCgat                                                        >  1:318159/1‑148 (MQ=255)
                                                                                                 ccTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATc                                                    <  4:109548/149‑1 (MQ=255)
                                                                                                     ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCCGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGtt                                                <  6:269132/149‑1 (MQ=255)
                                                                                                     ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAg                                                  <  7:259639/147‑1 (MQ=255)
                                                                                                     ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGtt                                                <  6:122515/149‑1 (MQ=255)
                                                                                                            cGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCg                                         <  7:77567/149‑1 (MQ=255)
                                                                                                             gTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCgg                                        >  1:192531/1‑149 (MQ=255)
                                                                                                                cccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGTGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGccc                                     <  5:244292/149‑1 (MQ=255)
                                                                                                                 ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCa                                                   <  8:103128/134‑1 (MQ=255)
                                                                                                                 ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCa                                                   >  7:103128/1‑134 (MQ=255)
                                                                                                                 ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCg                                    >  3:14412/1‑149 (MQ=255)
                                                                                                                       cGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCt                              >  5:121174/1‑149 (MQ=255)
                                                                                                                          ttGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGa                           >  7:111333/1‑149 (MQ=255)
                                                                                                                            gctgcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGa                           <  6:271227/147‑1 (MQ=255)
                                                                                                                             ctgcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGAtt                         <  2:119596/148‑1 (MQ=255)
                                                                                                                             ctgcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGAtt                         <  6:204949/148‑1 (MQ=255)
                                                                                                                               gcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGcc                      >  8:96079/1‑149 (MQ=255)
                                                                                                                                cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCCCGCACACATCGGTCACCACCCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                     >  6:42923/1‑149 (MQ=255)
                                                                                                                                cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                     <  5:273767/149‑1 (MQ=255)
                                                                                                                                cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                     >  7:32838/1‑149 (MQ=255)
                                                                                                                                       tCGTTGAAGGTGCCAGAGGAATCTGTCACCACGAAAACTTCAAAACCTTCCGCCAGCGCCTACAGGGTCGGGAAGGCCACTCAGACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGa               <  6:381028/148‑1 (MQ=255)
                                                                                                                                        cGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAg             >  8:101734/1‑149 (MQ=255)
                                                                                                                                         gTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGt            >  3:310275/1‑149 (MQ=255)
                                                                                                                                          ttGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTc           >  6:278779/1‑149 (MQ=255)
                                                                                                                                            gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTg                   <  1:79351/139‑1 (MQ=255)
                                                                                                                                            gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTg                   >  2:79351/1‑139 (MQ=255)
                                                                                                                                              tggTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCCCCACACGGGCGATGCTCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg       >  2:296962/2‑149 (MQ=255)
                                                                                                                                               ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCGTTGACGAAGTCTTCGt      >  6:86419/1‑149 (MQ=255)
                                                                                                                                               ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGt      >  7:355137/1‑149 (MQ=255)
                                                                                                                                               ggTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCAACACACCGGCGATGATCAGTTGCTTCCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGt      >  7:200541/1‑149 (MQ=25)
                                                                                                                                                gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt     >  1:148453/1‑149 (MQ=255)
                                                                                                                                                gTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt     >  6:137461/1‑149 (MQ=255)
                                                                                                                                                  gCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGt   >  7:213292/1‑149 (MQ=255)
                                                                                                                                                   cccGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTc  <  8:355137/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTC  >  NC_002947/5984548‑5984843

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