Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,070,583 T→C 5.2% K526E (AAG→GAG)  mdtB ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,070,5830TC 5.2% 230.1 / 2.1 76K526E (AAG→GAG) mdtBmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/32);  new base C (2/2);  total (42/34)
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.87e-01

GCGAACGAGGTGGACTGCGGCGCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCGTCAGGGTCAGGGACAC  >  NC_002947/4070435‑4070730
                                                                                                                                                     |                                                                                                                                                   
gCGAACGAGGTGGACTGCGGCGCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCtt                                                                                                                                                     <  4:311767/149‑1 (MQ=255)
  gAAGGAGGTGGACTGCGGAGCTTCGGAAATAGCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCGCTTGGGCACGACCATGTACAGGAACACCGTGTGCTGTAGGCTGGCCACCGCCACCAGCAAGGTCAGGTTCTGATGCttta                                                                                                                                                   <  4:314088/149‑3 (MQ=255)
  gAACGAGGTGGACTGCGGCGCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCtt                                                                                                                                                     >  1:41969/1‑147 (MQ=255)
   aaCGAGGTGGACTGCGGCGCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACTCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGa                                                                                                                                                   <  1:69112/149‑1 (MQ=255)
        ggTGGACTGCGGCGCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCGTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCAcc                                                                                                                                              >  3:64255/1‑148 (MQ=255)
                 cggcgCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAggg                                                                                                                                    <  1:164540/149‑1 (MQ=255)
                 cggcgCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAgg                                                                                                                                     <  2:110424/148‑1 (MQ=255)
                            aaaTACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAg                                                                                                                         <  4:64255/149‑1 (MQ=255)
                              aTACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTCGAGCACCCACTGCAGGGCACTGCCGTAGTg                                                                                                                       >  5:21242/1‑149 (MQ=255)
                              aTACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTCGAGCACCCACTGCAGGGCACTGCCGTAGTg                                                                                                                       >  5:25980/1‑149 (MQ=255)
                                 cccTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTg                                                                                                                    >  2:43997/1‑149 (MQ=255)
                                  ccTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTg                                                                                                                    >  8:132569/1‑148 (MQ=255)
                                   cTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGa                                                                                                                  <  1:193804/149‑1 (MQ=255)
                                      gATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGa                                                                                                                  <  1:188727/146‑1 (MQ=255)
                                      gATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGa                                                                                                                  >  2:188727/1‑146 (MQ=255)
                                      gATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGa                                                                                                                  >  2:193804/1‑146 (MQ=255)
                                            cccGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGt                                                                                                         >  3:222786/1‑149 (MQ=255)
                                                     cTGCACAGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAgg                                                                                                <  7:96201/149‑1 (MQ=255)
                                                           cGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCAc                                                                                          >  4:260668/1‑149 (MQ=255)
                                                           cGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCAc                                                                                          >  8:259526/1‑149 (MQ=255)
                                                             ggAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTg                                                                                        >  8:287348/1‑149 (MQ=255)
                                                                 gaagCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGt                                                                                                         >  4:66449/1‑128 (MQ=255)
                                                                 gaagCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGt                                                                                                         <  3:66449/128‑1 (MQ=255)
                                                                       cTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCTCCACCAGCAAGGTCAGCGGCTGATGCTTGAGAACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGTCGCGGTAg                                                                               >  8:162091/1‑148 (MQ=255)
                                                                           ggCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCg                                                                          <  3:260668/149‑1 (MQ=255)
                                                                           ggCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGAGGTAGAAgag                                                                          >  4:86718/1‑147 (MQ=255)
                                                                           ggCACGACCATGTACAGGAACACAGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCg                                                                          <  2:41969/149‑1 (MQ=255)
                                                                             cACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCa                                                                                           >  5:291912/1‑130 (MQ=255)
                                                                             cACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCa                                                                                           <  6:291912/130‑1 (MQ=255)
                                                                             cACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGc                                                                        <  7:132569/149‑1 (MQ=255)
                                                                              aCGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGcc                                                                       <  7:259526/149‑1 (MQ=255)
                                                                              aCGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGcc                                                                       >  1:245487/1‑149 (MQ=255)
                                                                                  ccATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGt                                                                   >  5:327941/1‑149 (MQ=255)
                                                                                    aTGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttc                                                                 >  7:43259/1‑149 (MQ=255)
                                                                                     tGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttct                                                                <  1:43997/149‑1 (MQ=255)
                                                                                     tGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCAGTAGTGCTTGATCAGCCAGTCGATCCAGGCGCCACTTGCGCGGTAGAAGCGGCCTTGttct                                                                <  1:204399/149‑1 (MQ=255)
                                                                                      gTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttctt                                                               >  5:266799/1‑149 (MQ=255)
                                                                                       tACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGcgcg                                                                                   <  1:127847/128‑1 (MQ=255)
                                                                                       tACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGcgcg                                                                                   >  2:127847/1‑128 (MQ=255)
                                                                                         cAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttcttctt                                                            <  7:287348/149‑1 (MQ=255)
                                                                                          aGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttcttcttc                                                           >  7:226831/1‑149 (MQ=255)
                                                                                          aGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttcttcttc                                                           >  6:11717/1‑149 (MQ=255)
                                                                                          aGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttcttcttc                                                           >  5:34323/1‑149 (MQ=255)
                                                                                            gAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCt                                                         >  4:8052/1‑149 (MQ=255)
                                                                                            gAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCt                                                         >  3:125915/1‑149 (MQ=255)
                                                                                              acacCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTg                                                       <  4:222786/149‑1 (MQ=255)
                                                                                               cacCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGt                                                                                                         <  5:233779/98‑1 (MQ=255)
                                                                                               cacCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGt                                                                                                         >  6:233779/1‑98 (MQ=255)
                                                                                                    tGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTcg                                                 >  3:198152/1‑149 (MQ=255)
                                                                                                    tGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTcg                                                 <  4:125915/149‑1 (MQ=255)
                                                                                                       gCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCgcgc                                              <  1:48635/149‑1 (MQ=255)
                                                                                                            cAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCtt                                            >  3:313332/1‑147 (MQ=255)
                                                                                                            cAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGa                                          <  3:8052/149‑1 (MQ=255)
                                                                                                            cAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGa                                          <  6:266799/149‑1 (MQ=255)
                                                                                                            cAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTCGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGa                                          <  6:21242/149‑1 (MQ=255)
                                                                                                            cAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTCGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGa                                          <  6:25980/149‑1 (MQ=255)
                                                                                                               gcttggcaccgccaccAGCAAGGTCAGCGGCTGATGCTTGAGCACCGACTGCAGGGCACTGCCGTAGTGCTGGAGCAGCCAGTCGATCCAGGCGCGACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCa                                       <  3:136667/143‑1 (MQ=255)
                                                                                                               gCTCGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGCTCAGCCAGTCGATCCAGGCGCCACAGGCCCGGTAGAAGCGGCCTTCTTCTTCTTCCTTGGGCTCGCGCTTGAg                                         >  4:136667/1‑147 (MQ=255)
                                                                                                                 tGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGa                                     >  2:137167/1‑149 (MQ=255)
                                                                                                                   gccaccgccaccAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGAc                                    <  1:137167/148‑1 (MQ=255)
                                                                                                                    ccaccgccaccAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGAcgc                                  >  4:276390/1‑149 (MQ=255)
                                                                                                                         gccaccAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGca                             <  6:327941/149‑1 (MQ=255)
                                                                                                                          ccaccaGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGcac                            <  8:226831/149‑1 (MQ=255)
                                                                                                                           caccaGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGcaca                           >  6:8461/1‑149 (MQ=255)
                                                                                                                                gCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACatcatc                      <  8:43259/149‑1 (MQ=255)
                                                                                                                                 cAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGCAGGGGCCTTGTTCTTCTTCCTTGGCCTAGCGCTTGAGCAGAAGCGCGCACATCATCg                     >  7:236101/1‑149 (MQ=255)
                                                                                                                                        aGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTg                                                                    >  2:217657/1‑95 (MQ=255)
                                                                                                                                        aGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCgtcag              <  5:11717/149‑1 (MQ=255)
                                                                                                                                         gCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttcttctt                                                            <  8:69203/102‑1 (MQ=255)
                                                                                                                                         gCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGttcttctt                                                            >  7:69203/1‑102 (MQ=255)
                                                                                                                                         gCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCgtcagg             >  7:195962/1‑149 (MQ=255)
                                                                                                                                         gCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCgtcagg             >  8:260295/1‑149 (MQ=255)
                                                                                                                                               gATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCgtcagggtcagg       >  5:81165/1‑149 (MQ=255)
                                                                                                                                                aTGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCgtcagggtcagg       >  2:49113/1‑148 (MQ=255)
                                                                                                                                                  gCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCGTCAGGGTCAGGGa     <  6:34323/148‑1 (MQ=255)
                                                                                                                                                  gCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCGTCAGGGTCAGGGa     >  4:263513/1‑148 (MQ=255)
                                                                                                                                                   cTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACgcgcgc                              <  5:315242/122‑1 (MQ=255)
                                                                                                                                                   cTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACgcgcgc                              >  6:315242/1‑122 (MQ=255)
                                                                                                                                                   cTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCCGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCGTCAGGGTCAGGGaca   >  7:7455/1‑149 (MQ=255)
                                                                                                                                                    ttGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCGTCAGGGTCAGGGacac  >  7:176785/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGAACGAGGTGGACTGCGGCGCTTCGGAAATACCCTGGATCACCCCGGTATCCTGCACCGGGAAGAAGCCCTTGGGCACGACCATGTACAGGAACACCGTGAGCACCAGGCTGGCCACCGCCACCAGCAAGGTCAGCGGCTGATGCTTGAGCACCCACTGCAGGGCACTGCCGTAGTGCTGGATCAGCCAGTCGATCCAGGCGCCACTGGCGCGGTAGAAGCGGCCTTGTTCTTCTTCCTTGGGCTCGCGCTTGAGCAGACGCGCGCACATCATCGGCGTCAGGGTCAGGGACAC  >  NC_002947/4070435‑4070730

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: