Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,690,773 T→C 6.0% K612R (AAA→AGA)  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,7730TC 6.0% 197.2 / 3.4 67K612R (AAA→AGA) bcsBregulator of cellulose synthase, cyclic di‑GMP binding
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base C (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTT  >  NC_000913/3690628‑3690914
                                                                                                                                                 |                                                                                                                                             
cGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCttt                                                                                                                                              >  4:247296/1‑147 (MQ=255)
  cACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCTCTATCGCTCTCGTCCGGAACAATGCCGGGGAACGTGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAg                                                                                                                                           <  5:370454/148‑1 (MQ=255)
            ggTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCTGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCggg                                                                                                                                 <  3:247296/148‑1 (MQ=255)
            ggTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCggg                                                                                                                                 <  6:105367/148‑1 (MQ=255)
            ggTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCggg                                                                                                                                 <  1:76094/148‑1 (MQ=255)
            ggTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCggg                                                                                                                                 <  3:226766/148‑1 (MQ=255)
                     tGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGt                                                                                                                                      <  1:319675/134‑1 (MQ=255)
                     tGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGt                                                                                                                                      >  2:319675/1‑134 (MQ=255)
                     tGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATaccacc                                                                                                                         >  6:31759/1‑147 (MQ=255)
                      gACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCg                                                                                                                        <  7:156362/147‑1 (MQ=255)
                               tCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCatg                                                                                                               <  8:74020/147‑1 (MQ=255)
                               tCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCatg                                                                                                               <  2:128266/147‑1 (MQ=255)
                                 tGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCatgatg                                                                                                            >  4:86290/1‑148 (MQ=255)
                                 tGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCatgatg                                                                                                            >  4:307922/1‑148 (MQ=255)
                                  gCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTGTCAGTTTGTCCGGGATACCACCGATGATCATGATGt                                                                                                           >  1:187119/1‑148 (MQ=255)
                                   ccGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGt                                                                                                           <  5:31759/147‑1 (MQ=255)
                                   ccGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTc                                                                                                          <  6:303343/148‑1 (MQ=255)
                                   ccGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTGCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTc                                                                                                          <  2:130838/148‑1 (MQ=255)
                                            cGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCttt                                                                                                 >  2:347649/1‑148 (MQ=255)
                                             gCTCTCGTCCGGCACAATGCCGGGGAACGGGGGCTGGCGCAGCGGTGGTTTCACCCAGCTTTCGGGCGCCTGCCCCAATAGGGCGAGCTGCTTCTCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTg                                                                                                <  2:354878/148‑1 (MQ=255)
                                                 tcGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGcc                                                                                              <  1:347649/146‑1 (MQ=255)
                                                   gTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGccc                                                                                             <  3:307922/145‑1 (MQ=255)
                                                    gccGGCACAATGCCGGGGAACGGGTTCTGGCGCATCGGTGTTTGCACCCAGCTTTCGGTCGCCTGCACCAATAGCTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGaa                                                                                         <  6:257370/147‑1 (MQ=255)
                                                       ggCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCgg                                                                                                        >  7:297643/1‑130 (MQ=255)
                                                       ggCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCgg                                                                                                        <  8:297643/130‑1 (MQ=255)
                                                        gCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGt                                                                                     <  3:86290/148‑1 (MQ=255)
                                                               gCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGt                                                                                                                                          >  7:336142/1‑88 (MQ=255)
                                                               gCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGt                                                                                                                                          <  8:336142/88‑1 (MQ=255)
                                                                     gAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCg                                                                        >  7:182025/1‑148 (MQ=255)
                                                                            gTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCg                                                                 >  4:68788/1‑148 (MQ=255)
                                                                              cTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTc                                                               <  5:332412/148‑1 (MQ=255)
                                                                                 gcgcATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAg                                                             <  7:312199/147‑1 (MQ=255)
                                                                                       cGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCatg                                                                                                               >  3:20/1‑91 (MQ=255)
                                                                                       cGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCatg                                                                                                               <  4:20/91‑1 (MQ=255)
                                                                                       cGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGa                                                                                          <  8:131414/112‑1 (MQ=255)
                                                                                       cGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGa                                                                                          >  7:131414/1‑112 (MQ=255)
                                                                                       cGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCAACGATGATCATGATGTCGGCATATTTTCACTGAATGGTGCTGCCTTCATCGGTCCCCGACAGATTAATc                                                      >  6:297282/1‑148 (MQ=255)
                                                                                          tgtTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGcc                                                   <  3:68788/148‑1 (MQ=255)
                                                                                           gtTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCg                                                  >  2:274723/1‑148 (MQ=255)
                                                                                             tttCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCgg                                                 >  1:2661/1‑147 (MQ=255)
                                                                                               tCACCCAGCTTTCGGTCTCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGa                                               <  6:86902/147‑1 (MQ=255)
                                                                                                  cccAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGcc                                           >  8:86060/1‑148 (MQ=255)
                                                                                                    cAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCg                                         >  8:116400/1‑148 (MQ=255)
                                                                                                       cTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCTTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCTCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCt                                      >  8:50307/1‑148 (MQ=255)
                                                                                                        tttCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCtg                                     >  2:97584/1‑148 (MQ=255)
                                                                                                             gTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTgcgc                                <  5:297282/148‑1 (MQ=255)
                                                                                                                   tGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCTATGATCATGATGTCTGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTACTCGCCGGGACGCCCGTCTGTGCGCCGATa                           >  6:97948/1‑147 (MQ=255)
                                                                                                                     cACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCTCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATaaa                         <  7:50307/147‑1 (MQ=255)
                                                                                                                     cACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATaaa                         <  1:274723/147‑1 (MQ=255)
                                                                                                                        cAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCa                     <  7:116400/148‑1 (MQ=255)
                                                                                                                        cAATAGGTCGATCTGCTTATCGTCTGTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGATTGGTGCTGCCCTCAGCGTTCCCCGGCAGCTTAGTCGCCGGGAAGCCCGGCTGTGGGCCGAGAAACCCa                     >  2:292102/1‑148 (MQ=255)
                                                                                                                         aaTAGGTCGATCTGCTTATCGTCTCTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCg                                                                        <  3:18638/96‑1 (MQ=255)
                                                                                                                         aaTAGGTCGATCTGCTTATCGTCTCTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCg                                                                        >  4:13404/1‑96 (MQ=255)
                                                                                                                         aaTAGGTCGATCTGCTTATCGTCTCTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCg                                                                        >  4:18638/1‑96 (MQ=255)
                                                                                                                         aaTAGGTCGATCTGCTTATCGTCTCTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCg                                                                        <  3:13404/96‑1 (MQ=255)
                                                                                                                          aTAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAAc                   <  7:86060/148‑1 (MQ=255)
                                                                                                                              gTCGATCTGCTTATCGTCTTTCAGTTTGTCGGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCt                                      <  2:100533/125‑1 (MQ=255)
                                                                                                                              gTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCt                                      >  1:100533/1‑125 (MQ=255)
                                                                                                                               tCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTAt              >  3:153359/1‑148 (MQ=255)
                                                                                                                               tCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGACATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGAtt                                                          >  5:357298/1‑104 (MQ=255)
                                                                                                                               tCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGACATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGAtt                                                          <  6:357298/104‑1 (MQ=255)
                                                                                                                                cGAGCTGCTTACCGTCTTTCCGCTCGTCCGGCATACCACCGATCCCCATGATGTCGGCATCTTGGCCCTGAATCGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGACCCCCGGCTGCGCGCCGATAAAACCAACAGTAtt             <  8:182025/148‑1 (MQ=255)
                                                                                                                                          aTCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGTAGCCCGTc                                       >  1:196903/1‑112 (MQ=255)
                                                                                                                                          aTCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGTAGCCCGTc                                       <  2:196903/112‑1 (MQ=255)
                                                                                                                                          aTCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGt    <  4:153359/147‑1 (MQ=255)
                                                                                                                                           tCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGttt  >  7:230893/1‑148 (MQ=255)
                                                                                                                                              tCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAAcc                      >  7:357039/1‑125 (MQ=255)
                                                                                                                                              tCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAAcc                      <  8:357039/125‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CGCACCGGAAGAGGTCAGCGTTGACCGGGTTTCTGCCGCGCGATCGCTCTCGTCCGGCACAATGCCGGGGAACGGGGTCTGGCGCATCGGTGTTTTCACCCAGCTTTCGGTCGCCTGCACCAATAGGTCGATCTGCTTATCGTCTTTCAGTTTGTCCGGGATACCACCGATGATCATGATGTCGGCATCTTTGCCCTGAATGGTGCTGCCATCATCGGTCACCGTCAGATTAATCGCCGGGAAGCCCGTCTGTGCGCCGATAAAACCAACAGTATTCAGCAACGTTT  >  NC_000913/3690628‑3690914

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: ‑