Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,065,815 T→A 8.2% S49C (AGC→TGC)  PP_3582 ← RND transporter outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,065,8150TA 8.2% 174.2 / 4.8 61S49C (AGC→TGC) PP_3582RND transporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base A (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 = 5.18e-01

GGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATA  >  NC_002947/4065672‑4065960
                                                                                                                                                |                                                                                                                                                 
ggTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCa                                                                                                                                              <  6:161446/149‑1 (MQ=255)
      gTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagca                                                                                                                                          <  1:210740/147‑5 (MQ=255)
          ctgcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCg                                                                                                                                    >  7:70516/1‑149 (MQ=255)
          ctgcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCg                                                                                                                                    <  3:107120/149‑1 (MQ=255)
            gcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggt                                                                                                                                  <  5:3558/149‑1 (MQ=255)
             cGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                                 <  3:202666/149‑1 (MQ=255)
             cGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                                 <  2:260550/149‑1 (MQ=255)
             cGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggt                                                                                                                                  <  2:138185/148‑1 (MQ=255)
                      tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                                 <  1:151163/140‑1 (MQ=255)
                      tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                                 >  2:151163/1‑140 (MQ=255)
                      tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCggg                                                                                                                        >  6:118303/1‑149 (MQ=255)
                        ttCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGc                                                                                                                      >  3:35498/1‑149 (MQ=255)
                          cGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGc                                                                                                                      <  2:242200/147‑1 (MQ=255)
                           gACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCa                                                                                                                   <  8:70516/149‑1 (MQ=255)
                             ccAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAgg                                                                                                                 >  1:172998/1‑149 (MQ=255)
                             ccAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAgg                                                                                                                 <  2:33365/149‑1 (MQ=255)
                              cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                                >  1:253068/1‑149 (MQ=255)
                              cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                                <  3:84115/149‑1 (MQ=255)
                              cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                                >  4:265640/1‑149 (MQ=255)
                               aGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGt                                                                                                               <  3:265640/149‑1 (MQ=255)
                                  gcATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACa                                                                                                            <  8:308197/149‑1 (MQ=255)
                                               cATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAg                                                                                               <  4:35498/149‑1 (MQ=255)
                                                aTCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAgg                                                                                              >  1:8669/1‑149 (MQ=255)
                                                aTCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAgg                                                                                              >  4:239986/1‑149 (MQ=255)
                                                    ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                                          >  6:70770/1‑149 (MQ=255)
                                                    ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                                          >  5:34316/1‑149 (MQ=255)
                                                    ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagc                                                                                                                                           >  2:72777/1‑97 (MQ=38)
                                                    ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagc                                                                                                                                           <  1:72777/100‑4 (MQ=38)
                                                     cGTAGATTTCCCACCAGGCGCCGCGGGCGAGTGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                                          >  1:202916/1‑149 (MQ=255)
                                                     cGTAGATTTACCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACg                                                                                         >  5:291212/1‑149 (MQ=255)
                                                          aTTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACaa                                                                                    >  3:313645/1‑149 (MQ=255)
                                                                caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaa                                                                                                                                             >  2:223162/1‑85 (MQ=37)
                                                                caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaa                                                                                                                                             <  1:223162/86‑2 (MQ=37)
                                                                caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCT‑GGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCa                                                                             >  6:78919/1‑149 (MQ=255)
                                                                     ggCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTTCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGa                                                                         <  2:172998/149‑1 (MQ=255)
                                                                          cgcgGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAgg                                                                    >  7:243689/1‑149 (MQ=255)
                                                                               gCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCTGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGa                                                               <  5:193199/149‑1 (MQ=255)
                                                                                 gA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgc                                                             <  2:253068/149‑1 (MQ=255)
                                                                                      gCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgccacg                                                         >  4:296604/1‑149 (MQ=255)
                                                                                                 cGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCa                                              <  3:239986/149‑1 (MQ=255)
                                                                                                 cGAGGCCTGGATCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGt                                                >  1:34058/1‑147 (MQ=255)
                                                                                                  gAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAg                                             >  8:153327/1‑149 (MQ=255)
                                                                                                   aGGCCT‑GGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTg                                           <  5:78919/149‑1 (MQ=255)
                                                                                                     gCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGc                                          >  8:144551/1‑149 (MQ=255)
                                                                                                     gCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGc                                          >  7:62655/1‑149 (MQ=255)
                                                                                                       cTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgccacggc                                                       >  3:277837/1‑134 (MQ=255)
                                                                                                       cTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgccacggc                                                       <  4:277837/134‑1 (MQ=255)
                                                                                                         gggTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGc                                      <  8:243689/149‑1 (MQ=255)
                                                                                                          ggTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCa                                     <  6:291212/149‑1 (MQ=255)
                                                                                                              cAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAg                                    <  5:337720/146‑1 (MQ=255)
                                                                                                              cAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAg                                    >  6:337720/1‑146 (MQ=255)
                                                                                                                  cTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACgg                             >  6:184843/1‑149 (MQ=255)
                                                                                                                   ttCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGt                            >  4:182979/1‑149 (MQ=255)
                                                                                                                           gCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTggg                    >  6:20865/1‑149 (MQ=255)
                                                                                                                               gAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGt                  <  5:118303/147‑1 (MQ=255)
                                                                                                                                aaCTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTg               <  2:8669/149‑1 (MQ=255)
                                                                                                                                 aCTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTgg              <  8:62655/149‑1 (MQ=255)
                                                                                                                                 aCTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGATAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCtt                                        <  6:122926/123‑1 (MQ=255)
                                                                                                                                       ccgccgTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCaaaa        <  5:184843/149‑1 (MQ=255)
                                                                                                                                         gccgTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGt      <  8:290875/149‑1 (MQ=255)
                                                                                                                                             tGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATa  <  4:161069/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATA  >  NC_002947/4065672‑4065960

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