Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 538,397 A→G 10.5% I259V (ATC→GTC)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947538,3970AG10.5% 173.9 / 8.9 66I259V (ATC→GTC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (3/4);  total (31/35)
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.99e-01

GCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGTGCA  >  NC_002947/538250‑538545
                                                                                                                                                   |                                                                                                                                                    
gCCTCGGTGCTGCTGCGCGCCCTGGGTTAGAGCACTGAAGAAGTGCTGAACACCTTCTACACCAACAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACAt                                                                                                                                                     <  4:103684/149‑1 (MQ=255)
gCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACAt                                                                                                                                                     <  8:232172/149‑1 (MQ=255)
 ccTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATc                                                                                                                                                    <  3:106245/149‑1 (MQ=255)
  cTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATcc                                                                                                                                                   <  5:191869/149‑1 (MQ=255)
  cTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATc                                                                                                                                                    >  6:97349/1‑148 (MQ=255)
        gctgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGAc                                                                                                                                              >  1:336090/1‑148 (MQ=255)
         ctgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActg                                                                                                                                            <  1:339837/149‑3 (MQ=255)
          tgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGa                                                                                                                                            >  1:51427/1‑148 (MQ=255)
          tgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGa                                                                                                                                            >  6:80429/1‑148 (MQ=255)
            ctgcGCGCCCTGGGTTACAGCACTTAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGaaa                                                                                                                                          <  5:23822/148‑1 (MQ=255)
             tgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAAcc                                                                                                                                        >  3:75048/1‑149 (MQ=255)
              gcgcgcCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgca                                                                                                                                          <  2:143685/146‑5 (MQ=255)
              gcgcgcCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgca                                                                                                                                          >  1:143685/1‑142 (MQ=255)
               cgcgcCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCGACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCgg                                                                                                                                      >  1:241311/1‑149 (MQ=255)
                  gcCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActg                                                                                                                                            <  3:51680/140‑3 (MQ=255)
                  gcCCTCGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActg                                                                                                                                            >  4:51680/1‑138 (MQ=255)
                   cccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCaaa                                                                                                                                  <  5:80429/149‑1 (MQ=255)
                    ccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGAct                                                                                                                                             >  7:302533/1‑136 (MQ=255)
                    ccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGAct                                                                                                                                             <  8:302533/137‑2 (MQ=255)
                         gTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtc                                                                                                                            <  5:97349/149‑1 (MQ=255)
                           tACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtcgt                                                                                                                          >  4:267985/1‑149 (MQ=255)
                               gCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAg                                                                                                                      <  3:267985/149‑1 (MQ=255)
                                  cTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCaa                                                                                                                   <  7:48168/149‑1 (MQ=255)
                                                 acacCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTAcc                                                                                                    >  3:164408/1‑149 (MQ=255)
                                                 acacCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTAcc                                                                                                    <  7:239200/149‑1 (MQ=255)
                                                    ccTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAgg                                                                                                                 <  4:157735/133‑1 (MQ=255)
                                                    ccTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAgg                                                                                                                 >  3:157735/1‑133 (MQ=255)
                                                     cTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc                                                                                                >  3:150248/1‑149 (MQ=255)
                                                     cTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc                                                                                                <  2:336090/149‑1 (MQ=255)
                                                      ttCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgcg                                                                                               >  7:208221/1‑149 (MQ=255)
                                                           caccacCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCcac                                                                                           <  2:241311/148‑1 (MQ=255)
                                                              caccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATc                                                                                        >  5:33623/1‑148 (MQ=255)
                                                               accgACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCaa                                                                                      >  2:12642/1‑149 (MQ=255)
                                                                ccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAAc                                                                                     <  2:51427/149‑1 (MQ=255)
                                                                 caACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAt                                                                                                                             <  7:318001/108‑1 (MQ=255)
                                                                 caACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAt                                                                                                                             >  8:318001/1‑108 (MQ=255)
                                                                 caACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAAc                                                                                     >  8:53802/1‑148 (MQ=255)
                                                                   aCGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAg                                                                                  >  2:312236/1‑149 (MQ=255)
                                                                            acaTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAg                                                                         <  4:164408/149‑1 (MQ=255)
                                                                            acaTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAg                                                                         >  1:93934/1‑149 (MQ=255)
                                                                                   cGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAg                                                                         >  8:96600/1‑142 (MQ=255)
                                                                                   cGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAg                                                                         <  7:96600/142‑1 (MQ=255)
                                                                                   cGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGc                                                                   <  7:174056/148‑1 (MQ=255)
                                                                                     gCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGt                                                                 >  2:78341/1‑148 (MQ=255)
                                                                                     gCGAAAAACTCAGCCTGGAACTGGTACCGCAGAGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTc                                                                <  1:53702/149‑1 (MQ=255)
                                                                                       gAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCaa                                                              >  5:292732/1‑149 (MQ=255)
                                                                                        aaaaaCTCAGCCTGGGACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAg                                                             <  1:40679/149‑1 (MQ=255)
                                                                                        aaaaaCTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCTCATCAACCAGCTCGAGAAGGCCGGCGTCAAg                                                             <  1:312236/149‑1 (MQ=255)
                                                                                         aaaaCTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGc                                                            <  6:292732/149‑1 (MQ=255)
                                                                                            aCTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAAc                                                         <  7:53802/149‑1 (MQ=255)
                                                                                                   cTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGt                                                  <  4:150248/149‑1 (MQ=255)
                                                                                                   cTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGt                                                  >  6:13451/1‑149 (MQ=255)
                                                                                                     ggAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTc                                                >  3:61879/1‑149 (MQ=255)
                                                                                                         cTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAAt                                            <  1:12642/149‑1 (MQ=255)
                                                                                                            gTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCTGCGTCAAGCAACTGGACGTTCCAATg                                           <  1:78341/147‑1 (MQ=255)
                                                                                                            gTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATgg                                          >  7:93965/1‑148 (MQ=255)
                                                                                                            gTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATg                                           <  8:208221/147‑1 (MQ=255)
                                                                                                                                     ttGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGc                 >  4:58410/1‑148 (MQ=255)
                                                                                                                                        cAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAg             >  4:11289/1‑149 (MQ=255)
                                                                                                                                        cAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCCCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCaa              <  6:33546/148‑1 (MQ=255)
                                                                                                                                         aGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAgg            >  6:202677/1‑149 (MQ=255)
                                                                                                                                         aGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAgg            <  4:75048/149‑1 (MQ=255)
                                                                                                                                           tCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGc           >  3:4242/1‑148 (MQ=255)
                                                                                                                                               ggACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCg      >  4:129913/1‑149 (MQ=255)
                                                                                                                                                gACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGt     <  2:93934/149‑1 (MQ=255)
                                                                                                                                                  cATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGTGc   <  5:202677/149‑1 (MQ=255)
                                                                                                                                                   aTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGTGCa  >  6:18796/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATCGTGCA  >  NC_002947/538250‑538545

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: