Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,046,227 C→G 7.9% M717I (ATG→ATC)  PP_0906 ← RND family multidrug transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,046,2270CG 7.9% 181.8 / 7.1 63M717I (ATG→ATC) PP_0906RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (2/3);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

CACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGT  >  NC_002947/1046083‑1046371
                                                                                                                                                |                                                                                                                                                
cACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATc                                                                                                                                               <  4:322270/148‑1 (MQ=255)
 aCGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGtcgatcga                                                                                                                                             >  2:232505/1‑149 (MQ=255)
  cGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCCACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGGGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGac                                                                                                                                            <  3:26734/149‑1 (MQ=255)
   gtagtaGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATc                                                                                                                                               >  1:158999/1‑145 (MQ=255)
   gtagtaGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATc                                                                                                                                               <  2:158999/145‑1 (MQ=255)
      gtagtTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACAccc                                                                                                                                        >  2:217671/1‑149 (MQ=255)
      gtagtTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACAcc                                                                                                                                         <  1:120340/148‑1 (MQ=255)
      gtagtTGTCCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTCTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACAccc                                                                                                                                        <  4:313896/149‑1 (MQ=255)
       tagtTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccct                                                                                                                                       >  4:1401/1‑148 (MQ=255)
         gttgttCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccctcg                                                                                                                                     <  3:1401/149‑4 (MQ=255)
         gttgttCAGCCAGCCGGGGTTGTAGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGACGCCGAGGGTGTCCATCGACACCcaca                                                                                                                                     >  5:205678/1‑149 (MQ=255)
            gttCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGGGGGTGTCGATGGACAccctcggc                                                                                                                                   <  5:251508/148‑6 (MQ=255)
              tCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAGGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTgc                                                                                                                                >  6:299425/1‑149 (MQ=255)
                  ccagccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACAccc                                                                                                                                        >  2:13427/1‑137 (MQ=255)
                  ccagccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGcgccg                                                                                                                            >  7:292401/1‑149 (MQ=255)
                    agccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGAGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGcgccgcg                                                                                                                          >  8:195811/1‑148 (MQ=255)
                      ccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCtt                                                                                                                        <  1:161053/149‑1 (MQ=255)
                       cGGGGTTGTCGCTGGAAGGCCGTTCAACCTGGGCGATCACCGTGTCACTGCCGCCCTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCttt                                                                                                                       <  5:318746/149‑1 (MQ=255)
                          ggTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGcc                                                                                                                    <  4:118258/149‑1 (MQ=255)
                                  cgGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAGGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGAt                                                                                                            <  5:299425/149‑1 (MQ=255)
                                  cgGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGAt                                                                                                            <  8:292401/149‑1 (MQ=255)
                                       aGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGGGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTcca                                                                                                       <  7:315932/149‑1 (MQ=255)
                                         tccGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccac                                                                                                      >  6:4601/2‑148 (MQ=255)
                                         gCCGATCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccac                                                                                                      >  3:147266/1‑148 (MQ=255)
                                           cGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTGGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGCGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccaccac                                                                                                   <  5:278660/149‑1 (MQ=255)
                                            gTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccctc                                                                                                                                      >  5:79155/1‑111 (MQ=255)
                                            gTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccctc                                                                                                                                      <  6:79155/113‑3 (MQ=255)
                                            gTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCTGTCGATGTccaccacc                                                                                                  <  6:39511/149‑1 (MQ=255)
                                                aaCCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACtt                                                                                              >  5:143892/1‑149 (MQ=255)
                                                    tGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCggg                                                                                          >  6:134796/1‑149 (MQ=255)
                                                    tGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCggg                                                                                          >  3:39155/1‑149 (MQ=255)
                                                        cGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTg                                                                                      <  4:147266/149‑1 (MQ=255)
                                                           tCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCAAGCAGGGTCGACAAGGTGCCGCCGAGGGTTTCCATCGACACCCCCATCTGCGCCGCCTTTCCCCGGTCGATCACCACCACCACTTCGGCTTTgac                                                                                    >  2:365309/1‑146 (MQ=255)
                                                             aCCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCg                                                                                   >  4:63396/1‑147 (MQ=255)
                                                             aCCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGaa                                                                                 <  5:4601/149‑1 (MQ=255)
                                                               cTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGCCACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAgg                                                                               >  4:274208/1‑149 (MQ=255)
                                                               cTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAgg                                                                               <  3:274208/149‑1 (MQ=255)
                                                               cTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAgg                                                                               >  5:56678/1‑149 (MQ=255)
                                                               cTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACAACTTCGGGTTTGTCGAAgg                                                                               >  6:94931/1‑149 (MQ=255)
                                                                ttGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAgg                                                                               >  4:237055/1‑148 (MQ=255)
                                                                   tAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCaa                                                                            <  4:202168/148‑1 (MQ=255)
                                                                       tGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTc                                                                                                               >  7:39479/1‑109 (MQ=255)
                                                                       tGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTc                                                                                                               <  8:39479/109‑1 (MQ=255)
                                                                       tGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccaccac                                                                                                   >  4:327153/1‑121 (MQ=255)
                                                                       tGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccaccac                                                                                                   <  3:327153/121‑1 (MQ=255)
                                                                       tGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATc                                                                        >  4:95302/1‑148 (MQ=255)
                                                                            ccTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGACAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGCGTTTGTCGAAGGCAAGATCGATATc                                                                  >  1:283359/1‑149 (MQ=255)
                                                                                cacgggTCAACCGGTGGATTCCCGCCGCGCCCCGCCGGGCCGCCACGGTGCCGCCTAGGGTGTCCATCCACCCCCCCATCTGCGCCGCCTTTGCCGGGTCTATGCCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCgcg                                                               <  4:39155/145‑3 (MQ=255)
                                                                                cAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGa                                                                <  7:69324/147‑1 (MQ=255)
                                                                                   gggTGAAACGGTTGATTTCCGCCTCGCACAGCAGGGTCGACAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGTCAAGATCGATATCGAGGAAg                                                           >  3:204226/1‑149 (MQ=255)
                                                                                         aaCGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAAc                                                     <  3:63396/149‑1 (MQ=255)
                                                                                           cGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACtt                                                   <  5:222988/149‑1 (MQ=255)
                                                                                              ttGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGcc                                                <  6:143892/149‑1 (MQ=255)
                                                                                              ttGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGcc                                                >  3:322868/1‑149 (MQ=255)
                                                                                                gATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGccc                                               >  4:252160/1‑148 (MQ=255)
                                                                                                  tttCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGAt                                            >  3:244791/1‑149 (MQ=255)
                                                                                                        ccTCGCCCAGCCGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACAGCCCCATCTGCGCCGCCTTGGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGa                                      <  6:205678/149‑1 (MQ=255)
                                                                                                        ccTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGa                                      <  3:237055/149‑1 (MQ=255)
                                                                                                                      gTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGAt                        <  5:134796/149‑1 (MQ=255)
                                                                                                                          ccAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGa                                                             >  7:7022/1‑108 (MQ=255)
                                                                                                                          ccAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGa                                                             <  8:7022/108‑1 (MQ=255)
                                                                                                                               gTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCAGCACTTCGGGTTTTTCGAAGGCAAGATCGATATCGTGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGc               <  1:365309/149‑1 (MQ=255)
                                                                                                                                tGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATTTCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGcc              <  2:283359/149‑1 (MQ=255)
                                                                                                                                            tgtCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGATGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGt  <  6:56678/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGT  >  NC_002947/1046083‑1046371

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