Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,046,230 C→G 7.9% S716S (TCG→TCC)  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,2300CG 7.9% 185.1 / 4.8 61S716S (TCG→TCC) PP_0906RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

CACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCGG  >  NC_002947/1046083‑1046376
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                                           <  6:79155/113‑3 (MQ=255)
                                            gTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccctc                                                                                                                                           >  5:79155/1‑111 (MQ=255)
                                            gTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCTGTCGATGTccaccacc                                                                                                       <  6:39511/149‑1 (MQ=255)
                                                aaCCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACtt                                                                                                   >  5:143892/1‑149 (MQ=255)
                                                    tGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCggg                                                                                               >  3:39155/1‑149 (MQ=255)
                                                    tGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCggg                                                                                               >  6:134796/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                                                                                    >  5:56678/1‑149 (MQ=255)
                                                               cTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAgg                                                                                    <  3:274208/149‑1 (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)
                                                                                                                                                 ctCGACACCCCCAGCGGAGCCGCCTTTGCCCGGTCGATGTCCAGCACCACTTCGGGTTTGGCGAAGGAAATATAGATATCTAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCgg  <  4:204226/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAGGGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCGG  >  NC_002947/1046083‑1046376

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