Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,839,326 C→G 5.0% intergenic (‑62/‑198) ascG ← / → ascF asc operon transcriptional repressor; prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,839,3260CG 5.0% 409.4 / 10.1 134intergenic (‑62/‑198)ascG/ascFasc operon transcriptional repressor; prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base C (57/70);  new base G (4/3);  total (61/73)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 140 total aligned reads displayed.

GCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGT  >  NC_000913/2839178‑2839475
                                                                                                                                                    |                                                                                                                                                     
gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                      <  1:670275/150‑1 (MQ=255)
gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGGCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTATCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                      <  2:192558/150‑1 (MQ=255)
 cTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCCCCGTTTGCCGATTGTCCTTGcac                                                                                                                                                     >  2:175737/1‑150 (MQ=255)
 cTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGcac                                                                                                                                                     >  2:1858843/1‑150 (MQ=255)
   gACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACaa                                                                                                                                                   <  1:510671/150‑1 (MQ=255)
        aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                              <  2:552687/150‑1 (MQ=255)
        aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                              <  1:1502619/150‑1 (MQ=255)
        aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggc                                                                                                                                              <  1:1476668/150‑1 (MQ=255)
        aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggc                                                                                                                                              <  1:1536591/150‑1 (MQ=255)
        aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGACACTACCAGCATCGTCGTCATCATTTTCCTCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTGTGCCGATTGTGCAAGGACAATcggc                                                                                                                                              <  2:265790/150‑1 (MQ=255)
          ccATTACCTGAAAGCACGCGGGAACCGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGGCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggcgg                                                                                                                                            <  1:1756993/150‑1 (MQ=255)
          ccATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggcgg                                                                                                                                            >  2:95412/1‑150 (MQ=255)
                ccTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaaaa                                                                                                                                      >  2:78331/1‑150 (MQ=255)
                  tGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACCCGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAatt                                                                                                                                    >  1:1745118/1‑149 (MQ=255)
                   gAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGCTTGTCCTTGCACAATCGGCGGGAAAAatat                                                                                                                                   >  2:1858823/1‑150 (MQ=255)
                   gAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCAGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAatat                                                                                                                                   <  2:1519468/150‑1 (MQ=255)
                          cgcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTAGCACAATCGGCCGGAAAAATATTCAGGTg                                                                                                                            >  1:285990/1‑150 (MQ=255)
                            cgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACCCGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAc                                                                                                                          >  1:1817526/1‑150 (MQ=255)
                             gggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                                                         <  2:66172/150‑1 (MQ=255)
                             gggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                                                         <  2:1377873/150‑1 (MQ=255)
                             gggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTAATCCCTTTACACGCAATAAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                                                         >  1:661750/1‑150 (MQ=255)
                                   cGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTc                                                                                                                   <  1:1733986/150‑1 (MQ=255)
                                    ggTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTca                                                                                                                  >  1:985464/1‑150 (MQ=255)
                                     gTCGCTTTTTAAACACCGGCGCGCTTCTACACTTACAACATAGTCGTCATCATTATCATCCCTTTACACGCAATCAAACCACTGTACTGCAACCTTTGCCGATAGTCCATGCACAATCGGCAGTAAAAAAATTCACCTGAACGGTTAcaa                                                                                                                 >  2:558725/1‑149 (MQ=255)
                                       cGCATTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTAATCCCTTTACACGAAATCAACGCAGTGTACTGCACCGTTTTCCGATTGTCCTTGCACACTCGGCGGGAAAAATATTCATGTGACCGGTTTCACaa                                                                                                               >  2:744654/1‑150 (MQ=255)
                                          ttttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCCACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAata                                                                                                            <  2:1132827/150‑1 (MQ=255)
                                           tttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTAACAAatat                                                                                                           >  2:1185221/1‑150 (MQ=255)
                                                aaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaaa                                                                                                      >  1:575191/1‑150 (MQ=255)
                                                  ccccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAAt                                                                                                    <  1:1578052/150‑1 (MQ=255)
                                                  ccccGGCGCGCTTCGCCACTTACAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACAGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAAt                                                                                                    >  2:603283/1‑150 (MQ=255)
                                                   cccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATg                                                                                                   <  2:1653570/150‑1 (MQ=255)
                                                    ccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGa                                                                                                  <  1:1588168/150‑1 (MQ=255)
                                                           gcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTc                                                                                           <  1:1325390/150‑1 (MQ=255)
                                                            cTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCa                                                                                          >  1:651228/1‑150 (MQ=255)
                                                             ttCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATAATCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                                                         <  2:1006436/150‑1 (MQ=255)
                                                              tCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAct                                                                                        <  2:739970/150‑1 (MQ=255)
                                                              tCGCCACTTCCAGCATCGTCGTCATCATTTTCATACCTATACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCCGCAAAAATATTCAGGTGACCGGTTCCACAAATATAAAAAATGAACTATTCaat                                                                                        >  1:1166183/1‑148 (MQ=255)
                                                               cGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActc                                                                                       >  1:1912239/1‑150 (MQ=255)
                                                                gCCACTTCCAGCATCGTCGTCATCATTTTCTTCCCTTTACACGCCATCAACGCAGTGTACAGCACCGTTTGCCGATTGTACTTGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActct                                                                                      >  2:368148/1‑150 (MQ=255)
                                                                gCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActct                                                                                      >  2:1576641/1‑150 (MQ=255)
                                                                  cACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctct                                                                                    <  1:1668482/150‑1 (MQ=255)
                                                                     ttCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATGAACGCAGTGTACTGGACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATTAACAATTCACTCTCTTGc                                                                                 <  1:837317/150‑1 (MQ=255)
                                                                       ccAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCtt                                                                               >  2:1087483/1‑150 (MQ=255)
                                                                        catcATCGTCGGCATGATTTTAATCCCTTTACACGCAATCAACCCTTTGTACTTCTCCGTTTTGGGATTGGCGTTGCTCATTCGGCGGGAAAAATATTCAGTGGACCGGTTTCACTAATATTAAAAATCAACATTTCACTCTCTTGCttc                                                                              <  2:630768/147‑2 (MQ=255)
                                                                        cAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTa                                                                              <  2:959861/150‑1 (MQ=255)
                                                                        cAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTa                                                                              <  1:465357/150‑1 (MQ=255)
                                                                           cATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTAttt                                                                           <  1:1344568/150‑1 (MQ=255)
                                                                                  gtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAATATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACa                                                                    <  2:439985/150‑1 (MQ=255)
                                                                                  gtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTAACAAATATAAAAAATGCACAATTCACTCTCTTGCTTATTTAGTGACa                                                                    <  1:755302/150‑1 (MQ=255)
                                                                                    catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                                                  <  1:1851214/150‑1 (MQ=255)
                                                                                    catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                                                  >  2:216977/1‑150 (MQ=255)
                                                                                                   tttACACTAAATCATCTCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTTACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgt                                                   <  1:1113871/150‑1 (MQ=255)
                                                                                                   tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgt                                                   <  1:593835/150‑1 (MQ=255)
                                                                                                    tttcacGCATTCAAGGCTGTGTACTGCAACGTTTTCCGTTTGTCCTTGCACAATCGGCGGGAAAAATATTACTGTGAACTGTTTCACAATTATAAAAAATGAACAATTCACTCTCTTGCTTATTAAGTGACAACTATTCATGATTTtgtg                                                  <  1:905622/147‑1 (MQ=255)
                                                                                                    ttccacGCAATCAACGCAGTGTACTTCAGCGGTTGCCGTTTGTCCATGCACAAGCGGCGGGAAAAATAATCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                                  <  2:365157/147‑1 (MQ=255)
                                                                                                    ttACACGCAATCAACGCAGTGTTCAGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGTAAAATAATCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                                  <  2:1481810/150‑1 (MQ=255)
                                                                                                    ttACAAGCAATCAACGCAGTGTACTGCACCGTTTGCAGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTGCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                                  <  2:1163188/150‑1 (MQ=255)
                                                                                                      acacGCAATCAACGCAGTGTACTGCACCGATTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaa                                                <  1:375329/150‑1 (MQ=255)
                                                                                                       cacGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACATCTATTCATGATTTTGTGaaa                                               >  1:1446743/1‑150 (MQ=255)
                                                                                                       cacGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaaa                                               <  1:1549060/150‑1 (MQ=255)
                                                                                                       cacGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaaa                                               >  2:615292/1‑150 (MQ=255)
                                                                                                       cacGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaaa                                               >  1:734324/1‑150 (MQ=255)
                                                                                                         cGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAATcc                                             >  2:1776213/1‑150 (MQ=255)
                                                                                                             aTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACTATTCACTCTCTTGCTTATTTCGTGACAACTATTCATGATTTTGTGAAACCGGtt                                         <  2:1343999/150‑1 (MQ=255)
                                                                                                              tCAACGCAGTGTACTGCACAGTTTGCCGATAGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGATTCACAAATATAAAAAATGAAAAATTCACTCTCTAGCTAATTTAGTGAAAACTAATCTAGATTTTGTGAAACCGGttt                                        >  2:1581175/1‑150 (MQ=255)
                                                                                                               cAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                       <  1:1705535/150‑1 (MQ=255)
                                                                                                               cAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                       <  1:1799854/150‑1 (MQ=255)
                                                                                                                aaCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCt                                      >  2:921114/1‑150 (MQ=255)
                                                                                                                aaCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCt                                      >  2:1241112/1‑150 (MQ=255)
                                                                                                                 ccGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                                     <  2:844760/149‑1 (MQ=255)
                                                                                                                  cGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTa                                    <  2:1072600/150‑1 (MQ=255)
                                                                                                                        gtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                              >  2:143242/1‑150 (MQ=255)
                                                                                                                        gtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAAAATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                              <  1:117774/150‑1 (MQ=255)
                                                                                                                         tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGt                             >  2:415810/1‑150 (MQ=255)
                                                                                                                          aCTTCACCGTTTGCCGATTGTCCTTGCACAATCGGCGTGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCATATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGtt                            <  1:266949/150‑1 (MQ=255)
                                                                                                                            tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATCTTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTc                          <  2:493327/150‑1 (MQ=255)
                                                                                                                            tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTc                          >  2:456279/1‑150 (MQ=255)
                                                                                                                              cACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAg                        <  2:842017/150‑1 (MQ=255)
                                                                                                                                ccgattgccgattgTCCTTGCACAATAGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCa                      <  1:340924/150‑1 (MQ=255)
                                                                                                                                ccGTTTGCCTATTGTCATTTCAAAATCGGCGGGAAAAATATTCAGGTGACCTGTTTCACAAATATAAATAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCa                      <  1:225336/150‑1 (MQ=255)
                                                                                                                                ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCa                      <  1:824529/150‑1 (MQ=255)
                                                                                                                                  gTTTGCCGATTGTCCTTTCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATc                    <  1:907973/150‑1 (MQ=255)
                                                                                                                                   tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCg                   >  2:2008/1‑150 (MQ=255)
                                                                                                                                    ccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCgg                  >  1:43086/3‑150 (MQ=255)
                                                                                                                                     tGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGc                 <  1:216977/150‑1 (MQ=255)
                                                                                                                                      gCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCa                >  2:1488543/1‑150 (MQ=255)
                                                                                                                                        cGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCAtt              >  2:193472/1‑150 (MQ=255)
                                                                                                                                            tGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGTTTTCTTAATTCCGTTTCAGCATCGGCATTTTTc          >  1:193922/1‑150 (MQ=255)
                                                                                                                                            tGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTc          <  1:175737/150‑1 (MQ=255)
                                                                                                                                            tGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTc          <  2:661750/150‑1 (MQ=255)
                                                                                                                                             gTCCTTGCACAATCGGCGGGAAAAATATTCAGGTTACCGGTTCCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTAAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc         <  2:1166183/150‑1 (MQ=255)
                                                                                                                                             gTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc         >  2:440616/1‑150 (MQ=255)
                                                                                                                                              tCCTTGCACAATCGGCGGGAAAAATATTCCTGGGACCGGTTTCACAAATAAAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGTTTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCg        <  2:1227573/150‑1 (MQ=255)
                                                                                                                                               ccTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTTTGAATCCGGTTTCTTAATTACTTTTCAGCATCGGCATTTTTCCGt       >  2:935381/1‑150 (MQ=255)
                                                                                                                                               ccTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGt       >  1:401418/1‑150 (MQ=255)
                                                                                                                                                cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGta      >  2:1700881/1‑149 (MQ=255)
                                                                                                                                                cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCAAAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTTACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTc      <  2:79208/150‑1 (MQ=255)
                                                                                                                                                  tGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCAc    >  2:1622831/1‑150 (MQ=255)
                                                                                                                                                    cacaATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGt  >  2:461077/1‑150 (MQ=255)
                                                                                                                                                    cacaATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGt  >  1:872486/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGT  >  NC_000913/2839178‑2839475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: