Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,690,623 G→A 5.7% A662V (GCG→GTG)  bcsB ← regulator of cellulose synthase, cyclic di‑GMP binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,690,6230GA 5.7% 225.8 / 2.7 70A662V (GCG→GTG) bcsBregulator of cellulose synthase, cyclic di‑GMP binding
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

ACCGCGACCGAACCGAACATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATC  >  NC_000913/3690479‑3690768
                                                                                                                                                |                                                                                                                                                 
aCCGCGACCGAACCGAACATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCAt                                                                                                                                               <  5:86221/149‑1 (MQ=255)
  cgcgACCGAACCGAACATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCg                                                                                                                                             <  4:341154/149‑1 (MQ=255)
         gaaccgaacATGGTGGCTCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCgg                                                                                                                                       >  8:70082/1‑148 (MQ=255)
         gaaccgaacATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCACCATCGCACCgg                                                                                                                                       >  7:402071/1‑148 (MQ=255)
           accgaacATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAag                                                                                                                                    <  4:314975/149‑1 (MQ=255)
              gaacATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAgg                                                                                                                                 >  4:161361/1‑149 (MQ=255)
                 cATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGt                                                                                                                                >  8:266642/1‑147 (MQ=255)
                 cATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGt                                                                                                                                >  2:359493/1‑147 (MQ=255)
                  aTGGTGGCGCGTTTGCCGCTATCGTTCACCGCCTCGTGAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGACGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAg                                                                                                                             <  3:344600/149‑1 (MQ=255)
                    ggtggCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCg                                                                                                                           >  7:62466/1‑149 (MQ=255)
                     gtggCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCg                                                                                                                           <  7:7417/148‑1 (MQ=255)
                        gcgcGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGa                                                                                                                       >  2:337601/1‑149 (MQ=255)
                          gcgTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCACTCACCGCCGCCATCGCACCGGAAGCGGTCAGCGTTGAcc                                                                                                                     >  6:362306/1‑149 (MQ=255)
                           cgTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCg                                                                                                                    <  3:161361/149‑1 (MQ=255)
                              ttGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGt                                                                                                                 <  2:322502/149‑1 (MQ=255)
                                  cGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGCCTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTGGACCGGGTTTCt                                                                                                             >  6:403116/1‑149 (MQ=255)
                                   gCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCAGCCAACAGCGCAATCACGCTGCGCTGCTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTg                                                                                                            >  1:194063/1‑149 (MQ=255)
                                       tCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGcc                                                                                                          <  1:368016/147‑1 (MQ=255)
                                         gTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCACCATCGCACCGGAagag                                                                                                                                  >  4:250046/1‑121 (MQ=255)
                                         gTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCACCATCGCACCGGAagag                                                                                                                                  <  3:250046/121‑1 (MQ=255)
                                             acggtaTCGTTAAGCATTTCTTAACCGCGTGGGCTATCTGCCTACAGCGCAATCAGGCTGCTCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGGGATCGCACCGGAAGAGGGCAGCGTGGTCCGGGTTTCTGCCGCGCGATc                                                                                                  <  5:16271/144‑1 (MQ=255)
                                             aCCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATc                                                                                                  <  3:33079/149‑1 (MQ=255)
                                               cGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGa                                                                                                    <  6:171786/145‑1 (MQ=255)
                                               cGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGa                                                                                                    >  5:171786/1‑145 (MQ=255)
                                               cGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGc                                                                                                >  5:282822/1‑149 (MQ=255)
                                                 cATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGctc                                                                                              <  8:62466/149‑1 (MQ=255)
                                                 cATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGctc                                                                                              >  4:14153/1‑149 (MQ=255)
                                                 cATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGctc                                                                                              <  1:148229/149‑1 (MQ=255)
                                                 cATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGctc                                                                                              >  5:207893/1‑149 (MQ=255)
                                                      ttAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTc                                                                                         >  2:244474/1‑149 (MQ=255)
                                                       tAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGc                                                                                                <  7:21421/141‑1 (MQ=255)
                                                       tAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGc                                                                                                >  8:21421/1‑141 (MQ=255)
                                                         aGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCg                                                                                       <  7:266642/148‑1 (MQ=255)
                                                          gCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCCCCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCgg                                                                                      >  6:10080/1‑148 (MQ=255)
                                                           cATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGca                                                                                    <  6:207893/149‑1 (MQ=255)
                                                            aTTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCTATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGca                                                                                    >  4:81113/1‑148 (MQ=255)
                                                               tCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGCGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAAt                                                                                >  4:112880/1‑149 (MQ=255)
                                                                 aTAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGc                                                                              >  6:86911/1‑149 (MQ=255)
                                                                 aTAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCAAAATGc                                                                              >  6:27270/1‑149 (MQ=255)
                                                                  tAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGcc                                                                             >  5:227970/1‑149 (MQ=255)
                                                                    cccGCGTGGGCTAGCTGCCCCCAGCGCAATCACGCTGCGCGGGGCGTGATACGGCGACTGGAAGCCAATCCCCGCCGCCATCGCACCGGAAGAGGGCAGCGTGGACCGGGTTTCTGCCGCGCGAGCGCTCTCGTCCGGCACAATGCCg                                                                            <  7:290806/147‑1 (MQ=255)
                                                                     ccGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCggg                                                                          >  1:352990/1‑149 (MQ=255)
                                                                       gcgTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGa                                                                                                                       >  4:381766/1‑102 (MQ=255)
                                                                       gcgTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGa                                                                                                                       <  3:381766/102‑1 (MQ=255)
                                                                       gcgTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCgggg                                                                         >  8:266897/1‑148 (MQ=255)
                                                                          tGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACg                                                                     <  7:70082/149‑1 (MQ=255)
                                                                                 tCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCt                                                               <  1:359493/148‑1 (MQ=255)
                                                                                     ccAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCTCCGGAAGAGGTCAGCGTTGACCGGGGTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGTGTCTg                                                              <  5:10080/145‑1 (MQ=255)
                                                                                          aGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCa                                                         >  1:44542/1‑145 (MQ=255)
                                                                                          aGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCa                                                         <  2:44542/145‑1 (MQ=255)
                                                                                            cgcAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAGCGGGGTCTGGCGCATCGgtg                                                   >  7:36373/1‑149 (MQ=255)
                                                                                             gcAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCACCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGgtgt                                                  <  8:402071/149‑1 (MQ=255)
                                                                                              cAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGtt                                                 <  7:266897/149‑1 (MQ=255)
                                                                                                 tCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGtt                                                 <  7:159541/146‑1 (MQ=255)
                                                                                                 tCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGGGtt                                                 >  8:159541/1‑146 (MQ=255)
                                                                                                   aCGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCAc                                            <  5:69586/149‑1 (MQ=255)
                                                                                                          gctgGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCtt                                     >  1:18777/1‑149 (MQ=255)
                                                                                                           ctgGTCGTTATACGGCGACTGGACGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCCCCCAGCttt                                    >  3:56402/1‑149 (MQ=255)
                                                                                                              gTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCgg                                 >  3:148750/1‑149 (MQ=255)
                                                                                                              gTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCgg                                 >  5:239652/1‑149 (MQ=255)
                                                                                                              gTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGCGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCgg                                 >  3:229598/1‑149 (MQ=255)
                                                                                                                 gTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCg                              >  7:401022/1‑149 (MQ=255)
                                                                                                                  ttATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGc                             >  7:59168/1‑149 (MQ=255)
                                                                                                                      aCGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGCCCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTg                          >  5:10472/1‑148 (MQ=255)
                                                                                                                           gACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCaa                    <  5:86911/149‑1 (MQ=255)
                                                                                                                           gACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCaa                    <  8:59168/149‑1 (MQ=255)
                                                                                                                             cTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATa                  <  6:227970/149‑1 (MQ=255)
                                                                                                                                gAAGCCAATCACCGCCGCCATCGCACCGTAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGt               <  2:194063/149‑1 (MQ=255)
                                                                                                                                        tCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGc       <  2:18777/149‑1 (MQ=255)
                                                                                                                                          aCCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCt      >  5:282818/1‑148 (MQ=255)
                                                                                                                                            cgccgcCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGAt           >  1:248440/1‑141 (MQ=255)
                                                                                                                                            cgccgcCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGAt           <  2:248440/141‑1 (MQ=255)
                                                                                                                                             gccgccATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGATCGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATc  <  8:36373/149‑1 (MQ=255)
                                                                                                                                             gccgccATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATc  <  6:282822/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACCGCGACCGAACCGAACATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCACGCTGCGCTGGTCGTTATACGGCGACTGGAAGCCAATCACCGCCGCCATCGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATC  >  NC_000913/3690479‑3690768

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: