Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 687,473 C→T 5.5% S113L (TCG→TTG)  ydgK → Bcr/CflA family multidrug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947687,4730CT 5.5% 222.7 / 2.5 72S113L (TCG→TTG) ydgKBcr/CflA family multidrug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base T (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GACCGCTTTGGCCGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCC  >  NC_002947/687325‑687619
                                                                                                                                                     |                                                                                                                                                  
gACCGCTTTGGCCGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTc                                                                                                                                                    >  5:180131/1‑149 (MQ=255)
    gCTTTGGCCGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCgg                                                                                                                                                <  4:336335/149‑1 (MQ=255)
     cTTTGGCCGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCggg                                                                                                                                               <  5:14090/149‑1 (MQ=255)
      tttGGCCGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCATTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGc                                                                                                                                              >  5:7395/1‑149 (MQ=255)
           ccGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCg                                                                                                                                         >  8:384336/1‑149 (MQ=255)
            cgtcGCACACCTATGATGTTCTGCGTTACGCTGTTCACCCTTGCATCACTGGCCTGTGCCTACTCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGt                                                                                                                                        <  8:327814/146‑1 (MQ=255)
            cGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTTGCGGGCCATCGt                                                                                                                                        >  5:28380/1‑149 (MQ=255)
            cGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTTGCGGGCCATCGt                                                                                                                                        >  5:34412/1‑149 (MQ=255)
               cGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGt                                                                                                                                        >  8:156708/1‑146 (MQ=255)
               cGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCa                                                                                                                                      <  7:156708/148‑1 (MQ=255)
                gCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGc                                                                                                                                    >  3:46036/1‑149 (MQ=255)
                  aaaCCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGa                                                                                                                                  >  7:26706/1‑149 (MQ=255)
                   aaCCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGa                                                                                                                                  >  2:149960/1‑148 (MQ=255)
                       gCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTTGCGGGCCATCGTCAGCGACAAGt                                                                                                                             <  6:28380/149‑1 (MQ=255)
                       gCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTTGCGGGCCATCGTCAGCGACAAGt                                                                                                                             <  6:34412/149‑1 (MQ=255)
                                  gCGTTACGCTGTTCACCCTGGCATCATTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTa                                                                                                                  <  6:7395/149‑1 (MQ=255)
                                     ttACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCt                                                                                                               <  1:149960/149‑1 (MQ=255)
                                     ttACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTAATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGAGACCCGGTAGCt                                                                                                               >  5:349859/1‑149 (MQ=255)
                                     ttACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGAAACCATGGAGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTc                                                                                                                                                    >  7:354706/1‑112 (MQ=255)
                                       aCGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCtt                                                                                                              <  4:184036/148‑1 (MQ=255)
                                        cGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCa                                                                                                            <  4:19581/149‑1 (MQ=255)
                                            gTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGcc                                                                                                         >  7:72007/1‑148 (MQ=255)
                                               cACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCaaa                                                                                                      <  8:72007/148‑1 (MQ=255)
                                                aCCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCCAAgt                                                                                                    >  2:409140/1‑149 (MQ=255)
                                                aCCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAgt                                                                                                    <  8:26706/149‑1 (MQ=255)
                                                  ccTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCgg                                                                                                                    >  4:284793/1‑131 (MQ=255)
                                                  ccTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCgg                                                                                                                    <  3:284793/131‑1 (MQ=255)
                                                    tGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGAc                                                                                                                        >  3:5670/1‑125 (MQ=255)
                                                    tGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGAc                                                                                                                        <  4:5670/125‑1 (MQ=255)
                                                        aTCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGCTTTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCg                                                                                             >  1:106672/1‑149 (MQ=255)
                                                        aTCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGc                                                                                            <  5:380224/149‑1 (MQ=255)
                                                             tGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTa                                                                                                                  <  6:241649/122‑1 (MQ=255)
                                                             tGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTa                                                                                                                  >  5:241649/1‑122 (MQ=255)
                                                               gCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGAt                                                                                     >  5:13292/1‑149 (MQ=255)
                                                                 cTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGc                                                                                   <  2:49472/149‑1 (MQ=255)
                                                                      ccTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCaaa                                                                                                      <  2:284873/125‑1 (MQ=255)
                                                                      ccTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCaaa                                                                                                      >  1:284873/1‑125 (MQ=255)
                                                                      ccTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGt                                                                               >  1:405320/1‑148 (MQ=255)
                                                                         aCGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATg                                                                           >  4:389700/1‑149 (MQ=255)
                                                                              cccAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCgg                                                                                                                    >  1:263238/1‑103 (MQ=255)
                                                                              cccAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCgg                                                                                                                    <  2:263238/103‑1 (MQ=255)
                                                                                 aaCCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTc                                                                                              >  6:63936/1‑122 (MQ=255)
                                                                                 aaCCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTc                                                                                              <  5:63936/122‑1 (MQ=255)
                                                                                   ccTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGgcgc                                                                 <  6:180131/149‑1 (MQ=255)
                                                                                            cTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTgg                                                        <  2:405320/149‑1 (MQ=255)
                                                                                            cTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTgg                                                        >  3:84743/1‑149 (MQ=255)
                                                                                                gtggtggCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCAcc                                                    <  4:84743/149‑1 (MQ=255)
                                                                                                    tggCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGACGATCCTGGCACCGATg                                                >  5:189392/1‑149 (MQ=255)
                                                                                                        gcgc‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTtgg                                            >  8:133999/1‑149 (MQ=255)
                                                                                                          gc‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTCGTCATGGGCTTGGCGCCGAACCTGGCACCGATGTtgggt                                          >  7:63874/1‑149 (MQ=255)
                                                                                                           c‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCTAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTtgggtg                                         <  1:409140/149‑1 (MQ=255)
                                                                                                           c‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTggggtg                                         >  6:322393/1‑149 (MQ=255)
                                                                                                             ttCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGTACCGATGtt                                              >  6:399596/1‑143 (MQ=255)
                                                                                                             ttCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGTACCGATGtt                                              <  5:399596/143‑1 (MQ=255)
                                                                                                             ttCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTtgggtgg                                        <  4:46036/149‑1 (MQ=255)
                                                                                                                aTCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATgg                                                                          <  7:153589/112‑1 (MQ=255)
                                                                                                                aTCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATgg                                                                          >  8:153589/1‑112 (MQ=255)
                                                                                                                 tCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGAGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTgggtggggtg                                    >  5:114311/1‑149 (MQ=255)
                                                                                                                  ccAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGc                                   >  4:350002/1‑149 (MQ=255)
                                                                                                                   cAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCt                                  >  4:149544/1‑149 (MQ=255)
                                                                                                                      gCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGt                               <  3:350002/149‑1 (MQ=255)
                                                                                                                             gtggtTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTg                        >  4:389325/1‑149 (MQ=255)
                                                                                                                              tggtTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTgg                       <  3:389325/149‑1 (MQ=255)
                                                                                                                                gtTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGcc                     >  3:85475/1‑149 (MQ=255)
                                                                                                                                    tgCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGGGCTGGTGAACCTGGCCGGc                  <  3:220525/148‑1 (MQ=255)
                                                                                                                                          cATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTtg                                             <  6:220631/115‑1 (MQ=255)
                                                                                                                                          cATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTtg                                             >  5:220631/1‑115 (MQ=255)
                                                                                                                                            tGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCg         <  3:149544/149‑1 (MQ=255)
                                                                                                                                               tGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATc      <  7:133999/149‑1 (MQ=255)
                                                                                                                                                  tttCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTc   >  8:153689/1‑149 (MQ=255)
                                                                                                                                                  tttCGCGGGACATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTc   >  6:264928/1‑149 (MQ=255)
                                                                                                                                                    tCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTcc  <  2:106672/148‑1 (MQ=255)
                                                                                                                                                    tCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTc   <  6:13292/147‑1 (MQ=255)
                                                                                                                                                    tCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTc   <  3:389700/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GACCGCTTTGGCCGACGCAAACCGCTGATGTTCGGCGTTACGCTGTTCACCCTGGCATCACTGGCCTGTGCCTACGCCCCCAACCTCGATACCTTGGTGGTGGCGCGC‑TTCATCCAGGCACTGGGTGGTTGTGCCGGCATGGTGCTTTCGCGGGCCATCGTCAGCGACAAGTGCGACCCGGTAGCTTCAGCCAAAGTGTTTTCGCAACTGATGCTGGTCATGGGCTTGGCGCCGATCCTGGCACCGATGTTGGGTGGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCC  >  NC_002947/687325‑687619

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