Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,452,781 C→T 7.7% Q315* (CAG→TAG)  PP_1271 → putative Multidrug efflux MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,452,7810CT 7.7% 187.3 / 4.5 65Q315* (CAG→TAG) PP_1271putative Multidrug efflux MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base T (3/2);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CAACGAGTGGTCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGC  >  NC_002947/1452642‑1452928
                                                                                                                                           |                                                                                                                                                   
cAACGAGTGGTCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGTAGATCCAgg                                                                                                                                            >  5:292152/1‑149 (MQ=255)
cAACGAGTGGTCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAgg                                                                                                                                            >  1:308042/1‑149 (MQ=255)
 aaCGAGTGGTCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAggg                                                                                                                                           >  1:215389/1‑149 (MQ=255)
     aGTGGTCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGTAGATCCAGGGCTAc                                                                                                                                       <  6:292152/149‑1 (MQ=255)
          tCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGccc                                                                                                                                  >  7:45841/1‑149 (MQ=255)
            cGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCgg                                                                                                                                >  7:222285/1‑149 (MQ=255)
             gAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGgc                                                                                                                               <  7:269353/149‑1 (MQ=255)
              aGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGgcg                                                                                                                              >  1:281836/1‑149 (MQ=255)
              aGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGACGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGTAGATCCAGGGCTACCGCCCGgcg                                                                                                                              >  5:11380/1‑149 (MQ=255)
                 cgttgccgttTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAg                                                                                                                           <  5:220237/149‑1 (MQ=255)
                  gttgccgttTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGa                                                                                                                          >  8:44483/1‑149 (MQ=255)
                       cgttTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGc                                                                                                                     >  8:290466/1‑149 (MQ=255)
                       cgttTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGc                                                                                                                     >  6:327058/1‑149 (MQ=255)
                               cAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTtga                                                                                                              <  7:42891/148‑1 (MQ=255)
                                      gcTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCtg                                                                                                      >  6:281483/1‑149 (MQ=255)
                                          aTTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGACGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGTAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCtggt                                                                                                    <  6:11380/147‑1 (MQ=255)
                                           ttGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCCCAGCCCAGCCCGTTGATGATGCTGTTGGc                                                                                                 >  4:252331/1‑149 (MQ=255)
                                             gCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACGAGCCCGTTGATGATGCTGGTGGccc                                                                                               >  3:32388/1‑148 (MQ=255)
                                               aGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATg                                                                                             <  8:45841/149‑1 (MQ=255)
                                                  ggcgCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGCCCAGCCCGTTGATGATGCTGGTGGCCATGCCg                                                                                          >  5:183621/1‑149 (MQ=255)
                                                     gcAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGGGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCaa                                                                                       <  1:96086/149‑1 (MQ=255)
                                                     gcAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCaa                                                                                       <  7:44483/149‑1 (MQ=255)
                                                        aCTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACt                                                                                     >  3:129204/1‑148 (MQ=255)
                                                         cTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGa                                                                                   >  7:211545/1‑149 (MQ=255)
                                                          ttGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGAt                                                                                  >  4:91392/1‑149 (MQ=255)
                                                            gAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCg                                                                                >  5:129957/1‑149 (MQ=255)
                                                               tttCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACAGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGccc                                                                             >  2:225972/1‑149 (MQ=255)
                                                                 tCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTg                                                                           >  3:124049/1‑149 (MQ=255)
                                                                   gCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGc                                                                          <  2:195877/148‑1 (MQ=255)
                                                                     ccTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGc                                                                       >  3:49537/1‑149 (MQ=255)
                                                                         gTGCCCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGTAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAc                                                                   >  4:35962/1‑149 (MQ=255)
                                                                           gACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAcgg                                                                 <  2:215389/149‑1 (MQ=255)
                                                                             cccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCAt                                                                                              >  6:179439/1‑118 (MQ=255)
                                                                             cccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCAt                                                                                              <  5:179439/118‑1 (MQ=255)
                                                                             cccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAcggcg                                                               <  5:156610/149‑1 (MQ=255)
                                                                              ccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCGTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAcggcgg                                                              >  1:228934/1‑149 (MQ=255)
                                                                              ccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAcggcgg                                                              >  5:17474/1‑149 (MQ=255)
                                                                                tGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAcg                                                                  <  4:51260/143‑1 (MQ=255)
                                                                                tGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAcg                                                                  >  3:51260/1‑143 (MQ=255)
                                                                                  gCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCa                                                            <  7:290466/147‑1 (MQ=255)
                                                                                        gTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCaa                                                    <  2:281836/149‑1 (MQ=255)
                                                                                         tGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCaa                                                    >  6:249345/1‑148 (MQ=255)
                                                                                            tGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATc                                                >  4:43553/1‑149 (MQ=255)
                                                                                              gTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCcg                                              >  6:108144/1‑149 (MQ=255)
                                                                                                   gCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATcc                                               >  8:43620/1‑143 (MQ=255)
                                                                                                   gCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATcc                                               <  7:43620/143‑1 (MQ=255)
                                                                                                    cTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGCTCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGt                                        <  1:225972/149‑1 (MQ=255)
                                                                                                      ttCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTg                                       >  6:147446/1‑148 (MQ=255)
                                                                                                        cGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGa                                     <  6:17474/148‑1 (MQ=255)
                                                                                                        cGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGa                                     <  2:228934/148‑1 (MQ=255)
                                                                                                         gggTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGCTCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACt                                   >  2:57240/1‑149 (MQ=255)
                                                                                                              tgGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGc                              <  3:91392/149‑1 (MQ=255)
                                                                                                                 gcagcaTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGcc                                                                              <  6:202611/98‑1 (MQ=255)
                                                                                                                 gcagcaTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGcc                                                                              >  5:202611/1‑98 (MQ=255)
                                                                                                                  cagcaTCACCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACTGCGGCACTGTGCAACATCCGCGCAGTGGCCTGCCGCTggg                          >  6:173762/1‑149 (MQ=255)
                                                                                                                     caTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGc                       <  5:194264/149‑1 (MQ=255)
                                                                                                                       tCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTg                     <  5:249345/149‑1 (MQ=255)
                                                                                                                          ccTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACAACCCCGTTGATGATTCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCCGTGGACTGCCGCTGGGTGCTggg                   >  4:3263/1‑148 (MQ=255)
                                                                                                                           cTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCGACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCa                 >  1:72709/1‑149 (MQ=255)
                                                                                                                            tCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCAt                >  4:294112/1‑149 (MQ=255)
                                                                                                                              ggCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTg              <  4:129204/149‑1 (MQ=255)
                                                                                                                                  tACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCt          >  4:303906/1‑149 (MQ=255)
                                                                                                                                   aCCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTg         <  3:252331/149‑1 (MQ=255)
                                                                                                                                    ccTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCTACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTgg        <  2:308042/149‑1 (MQ=255)
                                                                                                                                     cTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGc       <  4:49537/149‑1 (MQ=255)
                                                                                                                                          gcAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGc  >  3:251038/1‑149 (MQ=255)
                                                                                                                                          gcAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCAGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGc  >  8:106914/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CAACGAGTGGTCCGAGCCGTTGCCGTTTTTCCAGTTGCGCTTATTGCAGCGGCGCAACTTGAGTTTCGCCCTGGTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGC  >  NC_002947/1452642‑1452928

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: