Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,068,693 C→A 8.7% R124L (CGC→CTC) ‡ mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,068,6930CA 8.7% 170.5 / 4.9 56R124L (CGC→CTC) ‡mdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base A (2/3);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGC  >  NC_002947/4068546‑4068837
                                                                                                                                                   |                                                                                                                                                
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                                                                                                                                        >  5:195714/1‑108 (MQ=255)
                                               cAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCGCAGCATGcccac                                                                                                                                        <  6:195714/111‑4 (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                                                                                       <  5:198211/136‑1 (MQ=255)
                                                                       cAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCt                                                                                       >  6:198211/1‑136 (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:26306/149‑1 (MQ=255)
                                                                                         cGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGc                                                        <  4:7701/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          <  3:261858/149‑1 (MQ=255)
                                                                                                                                       gTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCgt          <  2:79018/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             <  4:63723/136‑1 (MQ=255)
                                                                                                                                                 tgcGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAg             >  3:63723/1‑136 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCTACCCCCGACACCTGAGCCAGGCTCTGTGAAAGGATGGTGTCGGCCAGGTCATACAGCTGGCCTTTCTGCAACACATCCGAAGTCAGCGACAGCACCATGATCGGCGCCTGCGACGGGTTGATCTTCTTGTAGGTGGGCATGCTGCGCATGCCACTGGGCAGCAGGTTGCGGGTAGCGTTGATCGCCGCTTGAACCTCACGGGCTGCGCCATCGATGTCACGGCCCAACTCGAAGCCGATGACTACCCGTGTCGAACCCTGGTTGGAGCTGCTGGTGAGCGTGGTGACGC  >  NC_002947/4068546‑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: ‑