Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,068,691 T→G 8.2% S125R (AGC→CGC)  mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,068,6910TG 8.2% 168.4 / 5.5 58S125R (AGC→CGC) mdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base G (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGC  >  NC_002947/4068543‑4068837
                                                                                                                                                    |                                                                                                                                                  
tCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGCCGGGTTGATGTTCTTGTAGGGGGGCAtgct                                                                                                                                                    >  8:231409/1‑149 (MQ=255)
   ccTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTgcg                                                                                                                                                 >  1:347880/1‑149 (MQ=255)
    ctagcgccGACACCTGTGCCAGGCTCTGTGAAAGGAGGGTGTCGGCGAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGTCAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGCGCATGCTgcgc                                                                                                                                                <  1:67979/143‑1 (MQ=255)
       cccccGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCAt                                                                                                                                              <  2:354904/148‑1 (MQ=255)
        ccccGTCACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCCCagtg                                                                                                                                            <  3:28827/149‑4 (MQ=255)
          ccGACACCCGAGCCCGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGCTCTTCTTGTAGGTGGGCATGCGCAGCATGccc                                                                                                                                          <  4:68863/149‑2 (MQ=255)
            gACACCTGCGCCAGGCTCTGTGAAAGGATGGGGTCGGCCATGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACTGGTTGATCTTCTTGTAGGTGGGCATGCGCAGCATGcccac                                                                                                                                        <  2:60315/149‑4 (MQ=255)
                 cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCAtgctgc                                                                                                                                                  <  1:14717/134‑1 (MQ=255)
                 cTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCAtgctgc                                                                                                                                                  >  2:14717/1‑134 (MQ=255)
                    aGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGgcagca                                                                                                                                <  6:275798/149‑1 (MQ=255)
                     gCCAGGCTCGGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGTCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGTCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGcagcag                                                                                                                               <  6:352349/149‑1 (MQ=255)
                          gCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGc                                                                                                                          <  8:183024/149‑1 (MQ=255)
                               gtgAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCggg                                                                                                                       >  3:130213/1‑147 (MQ=255)
                                  aaaGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGc                                                                                                                   <  4:146108/148‑1 (MQ=255)
                                        tGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATc                                                                                                            >  1:79018/1‑149 (MQ=255)
                                               ggCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCtt                                                                                                     <  5:155150/149‑1 (MQ=255)
                                                gCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCtt                                                                                                     >  3:26306/1‑148 (MQ=255)
                                                gCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCAACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCtt                                                                                                     >  7:75361/1‑148 (MQ=255)
                                                  cAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCtt                                                                                                     <  5:255383/146‑1 (MQ=255)
                                                  cAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGAGGGTAGCGTTGATCGCCGCtt                                                                                                     >  6:255383/1‑146 (MQ=255)
                                                  cAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCGCAGCATGcccac                                                                                                                                        <  6:195714/111‑4 (MQ=255)
                                                  cAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCGCAGCATGcccac                                                                                                                                        >  5:195714/1‑108 (MQ=255)
                                                           cAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACggg                                                                                         <  3:261773/149‑1 (MQ=255)
                                                               tGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTgc                                                                                     >  2:17762/1‑149 (MQ=255)
                                                                 gCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCATGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTgcgc                                                                                   >  4:58155/1‑149 (MQ=255)
                                                                      tCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCAt                                                                                <  4:130213/147‑1 (MQ=255)
                                                                        tGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGAt                                                                            >  2:264449/1‑149 (MQ=255)
                                                                          cAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                       >  6:198211/1‑136 (MQ=255)
                                                                          cAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                       <  5:198211/136‑1 (MQ=255)
                                                                          cAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGt                                                                          >  6:339872/1‑149 (MQ=255)
                                                                           aaCACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTc                                                                         <  7:231409/149‑1 (MQ=255)
                                                                             cacaTCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCAc                                                                       >  3:5247/1‑149 (MQ=255)
                                                                             cacaTCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCAc                                                                       >  2:86343/1‑149 (MQ=255)
                                                                                   cGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGccc                                                                  >  3:7701/1‑148 (MQ=255)
                                                                                     aaGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAAc                                                               <  1:264449/149‑1 (MQ=255)
                                                                                      aGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACt                                                              <  1:37805/149‑1 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACg                                                                      >  2:350368/1‑137 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACg                                                                      <  1:350368/137‑1 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCg                                                            >  4:270399/1‑147 (MQ=255)
                                                                                          aGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGaa                                                          <  1:143288/149‑1 (MQ=255)
                                                                                            cGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGc                                                        <  4:7701/149‑1 (MQ=255)
                                                                                            cGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGc                                                        <  4:26306/149‑1 (MQ=255)
                                                                                             gACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGcc                                                       >  8:197343/1‑149 (MQ=255)
                                                                                                            ggCGCCTGCGACGGGTTGATCTTCATGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCgtgt                                        <  3:58155/149‑1 (MQ=255)
                                                                                                             gcgcCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTc                                       <  5:96856/149‑1 (MQ=255)
                                                                                                                ccTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGaa                                    >  8:143629/1‑149 (MQ=255)
                                                                                                                 cTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAAc                                   <  4:5247/149‑1 (MQ=255)
                                                                                                                 cTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAAc                                   >  4:261858/1‑149 (MQ=255)
                                                                                                                         ggTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTg                           <  1:86343/149‑1 (MQ=255)
                                                                                                                             gATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTgc                                                                                     <  5:229753/87‑1 (MQ=255)
                                                                                                                             gATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTgc                                                                                     >  6:229753/1‑87 (MQ=255)
                                                                                                                             gATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAgc                       >  7:282619/1‑149 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGCCTACCCGTGTCGACCCCGGGTTGGAGCTGCTGGTGCGCgt          >  4:205371/1‑149 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTgg                          >  2:182453/1‑133 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTgg                          <  1:182453/133‑1 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt          <  2:79018/149‑1 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt          <  3:261858/149‑1 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACAACCCGTGTCGAACCCTGGTTGGATCTGCTGGTGAGc            >  8:224396/1‑147 (MQ=255)
                                                                                                                                          gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTAAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACAACCCGTGTCGAACCCTGGTTGGAGCTCCTGGGGAGCgt          >  8:226475/1‑149 (MQ=255)
                                                                                                                                          agtggCATGCGCAGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt          >  2:97642/4‑149 (MQ=255)
                                                                                                                                                aTGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGAc    <  5:339872/149‑1 (MQ=255)
                                                                                                                                                  gctgcGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCAGGTTGGAGCTGCTGGTGAGCGTGGTGACGc  >  6:337766/1‑149 (MQ=255)
                                                                                                                                                    tgcGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAg             >  3:63723/1‑136 (MQ=255)
                                                                                                                                                    tgcGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAg             <  4:63723/136‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCCCCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGC  >  NC_002947/4068543‑4068837

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