Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,010,482 G→A 5.8% intergenic (‑1388/‑280) PP_2628 ← / → PP_2631 ABC transporter ATP‑binding protein/cellulose biosynthesis protein BcsF/YhjT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,010,4820GA 5.8% 207.6 / 2.0 69intergenic (‑1388/‑280)PP_2628/PP_2631ABC transporter ATP‑binding protein/cellulose biosynthesis protein BcsF/YhjT
Reads supporting (aligned to +/- strand):  ref base G (26/39);  new base A (2/2);  total (28/41)
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.82e-01

AGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGA  >  NC_002947/3010336‑3010626
                                                                                                                                                  |                                                                                                                                                
aGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGc                                                                                                                                                 >  8:259758/1‑148 (MQ=255)
 gcgcCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCCGCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGa                                                                                                                                               >  1:278282/1‑149 (MQ=255)
    ccACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGATCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACAt                                                                                                                                            <  2:31821/149‑1 (MQ=255)
    ccACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACAt                                                                                                                                            >  7:399642/1‑149 (MQ=255)
     cACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACAtt                                                                                                                                           <  5:188716/149‑1 (MQ=255)
       agggCGTTAGCCGTTGCCGTGGAAAGCATGTGGCACGGGGGGCGCAACGGCGAAATCGTCCACCAGTTTCTGGTGTTGCTTCCGGCGCCTGGCCTGACCCCGCTTCAGGCTGGCTCGCGTACGTTATTCCGCAGTGATGGCGACAtttc                                                                                                                                         <  8:265493/148‑3 (MQ=255)
            gTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCg                                                                                                                                    <  6:395678/149‑1 (MQ=255)
              tCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTg                                                                                                                                  >  6:222834/1‑149 (MQ=255)
                  cGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTg                                                                                                                              >  4:412565/1‑149 (MQ=255)
                       ccGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGaa                                                                                                                         <  6:84703/149‑1 (MQ=255)
                       ccGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGaa                                                                                                                         <  3:412565/149‑1 (MQ=255)
                       ccGGGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCGGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGa                                                                                                                          <  4:23148/148‑1 (MQ=255)
                        cGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTTGTTGGCGAAg                                                                                                                        >  3:337703/1‑149 (MQ=255)
                        cGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAg                                                                                                                        >  2:381331/1‑149 (MQ=255)
                         gTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGt                                                                                                                       <  7:42955/149‑1 (MQ=255)
                         gTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGt                                                                                                                       <  3:48787/149‑1 (MQ=255)
                          tGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTg                                                                                                                      <  5:13504/149‑1 (MQ=255)
                              aaGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgt                                                                                                                  >  8:393761/1‑149 (MQ=255)
                              aaGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgt                                                                                                                  >  8:92345/1‑149 (MQ=255)
                               aGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgtt                                                                                                                 >  5:114905/1‑149 (MQ=255)
                                gCATGTGGCACGGGGTGCGCAACGGCGAAATCGGCCACCAGTTGCTGGTGTTGCGGCCGGCGCCTGGCCTGACCCCGCTGCAGGCGTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttc                                                                                                                <  5:357653/149‑1 (MQ=255)
                                   tgtgGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctg                                                                                                             <  6:332981/149‑1 (MQ=255)
                                   tgtgGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctg                                                                                                             <  5:222834/149‑1 (MQ=255)
                                     tgGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctgt                                                                                                            <  5:279879/148‑1 (MQ=255)
                                       gCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctgttcc                                                                                                         >  2:351001/1‑149 (MQ=255)
                                       gCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctgttc                                                                                                          <  1:38842/148‑1 (MQ=255)
                                        cACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctgttcct                                                                                                        <  8:399642/149‑1 (MQ=255)
                                                         cGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTcc                                                                                        >  7:233534/1‑148 (MQ=255)
                                                           aaaTCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGaa                                                                                     <  1:351001/149‑1 (MQ=255)
                                                             aTCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAgg                                                                                   <  4:149886/149‑1 (MQ=255)
                                                                gTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAgggg                                                                                 <  2:81581/148‑1 (MQ=255)
                                                                gTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGt                                                                                <  4:274381/149‑1 (MQ=255)
                                                                 tCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTc                                                                               <  7:40929/149‑1 (MQ=255)
                                                                     ccaGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGc                                                                           <  3:155218/149‑1 (MQ=255)
                                                                      caGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAATTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGc                                                                           >  7:179831/1‑148 (MQ=255)
                                                                        gTTGCTGGTGTTGCGTCCGGCGCGTGGCCTTACCCGGCTGCAGGATTGGTCGCGTACGGTATTCCGCCGTGATGGCGACTTTGCCTGCGTGGTTGGCGAAGTTCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAg                                                                        <  1:313887/149‑1 (MQ=255)
                                                                            cTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAgcgc                                                                     <  8:52981/148‑1 (MQ=255)
                                                                               gtgtTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGa                                                                 >  7:339868/1‑149 (MQ=255)
                                                                                tgtTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGAc                                                                <  7:92345/149‑1 (MQ=255)
                                                                                tgtTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCAAGCAAGCGCTTGAc                                                                <  7:393761/149‑1 (MQ=255)
                                                                                 gtTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGAc                                                                >  5:360262/1‑148 (MQ=255)
                                                                                  ttGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGAcc                                                               <  1:381331/148‑1 (MQ=255)
                                                                                  ttGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACca                                                              >  1:168576/1‑149 (MQ=255)
                                                                                         cGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCtt                                                        <  6:114905/148‑1 (MQ=255)
                                                                                          ggCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTc                                                       >  5:101221/1‑148 (MQ=255)
                                                                                               ctggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTg                                                 >  7:197103/1‑149 (MQ=255)
                                                                                                tggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGt                                                >  2:254097/1‑149 (MQ=255)
                                                                                                tggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGACTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGt                                                >  2:260614/1‑149 (MQ=255)
                                                                                                   cctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCAt                                             <  5:211749/149‑1 (MQ=255)
                                                                                                    ctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGc                                                                           <  4:128920/118‑1 (MQ=255)
                                                                                                    ctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGc                                                                           >  3:128920/1‑118 (MQ=255)
                                                                                                        ccccGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTTGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGa                                          >  7:186054/1‑147 (MQ=255)
                                                                                                        ccccGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGGTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGaaa                                        <  2:278282/149‑1 (MQ=255)
                                                                                                             cTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCt                                   <  7:334473/149‑1 (MQ=255)
                                                                                                               gCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGt                                 <  6:360262/149‑1 (MQ=255)
                                                                                                                cAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGtt                                <  3:199093/149‑1 (MQ=255)
                                                                                                                    cTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTtcatc                            <  6:194634/149‑1 (MQ=255)
                                                                                                                           gcgTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATg                                            >  5:292055/1‑126 (MQ=255)
                                                                                                                           gcgTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATg                                            <  6:292055/126‑1 (MQ=255)
                                                                                                                            cgTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCagga                    <  8:179831/149‑1 (MQ=255)
                                                                                                                                   tATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACtt             >  5:236180/1‑149 (MQ=255)
                                                                                                                                   tATTCCGCCGTGATGACGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCa                                                     <  5:241704/109‑1 (MQ=255)
                                                                                                                                   tATTCCGCCGTGATGACGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCa                                                     <  5:235692/109‑1 (MQ=255)
                                                                                                                                   tATTCCGCCGTGATGACGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCa                                                     >  6:235692/1‑109 (MQ=255)
                                                                                                                                   tATTCCGCCGTGATGACGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCa                                                     >  6:241704/1‑109 (MQ=255)
                                                                                                                                     ttCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGc           <  8:339868/149‑1 (MQ=255)
                                                                                                                                       ccgccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCg         >  8:282905/1‑149 (MQ=255)
                                                                                                                                        cgccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCgg        <  7:259758/149‑1 (MQ=255)
                                                                                                                                          ccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCg      >  8:402194/1‑149 (MQ=255)
                                                                                                                                              gacggCGCCACTGCCTGCCTGGCTGGCGCAGCGCGCTGCCCCTTCCAGCTCGCCTGCCGCCCCGCCGGGCTTGAGCAAGCGCTTGACCACAACTCCCAGCTGTCATGGACAGAGCTGTTCATCAGCCAGGAGGTCCTTGCCGGCGCCGa  <  8:186054/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGA  >  NC_002947/3010336‑3010626

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: