Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,839,324 T→A 7.2% intergenic (‑60/‑200) 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,3240TA 7.2% 361.4 / 14.5 125intergenic (‑60/‑200)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 T (54/62);  new base A (6/3);  total (60/65)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

Only 100 of 126 total aligned reads displayed.

TGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCA  >  NC_000913/2839176‑2839472
                                                                                                                                                    |                                                                                                                                                    
  gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                   <  1:1006965/150‑1 (MQ=255)
  gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                   <  1:793096/150‑1 (MQ=255)
  gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                   <  1:786277/150‑1 (MQ=255)
   cTGACTTAGCCATTACCTGAAAGGACGCGGGAAACTGTCTCTTTTGAAACCCCGGCGCGCTTCGACACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAAAGCATTGTACTGCACCGTTTGCCGATTGTCCTTGcac                                                                                                                                                  <  2:1661775/150‑1 (MQ=255)
       cGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACCCGCAATCAGCGCAGTGTACTGCACGGTTTGCCGATTGTCGTTGCACAATc                                                                                                                                              >  1:183589/1‑150 (MQ=255)
          aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACTCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggc                                                                                                                                           <  2:878099/150‑1 (MQ=255)
          aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                           <  1:890228/150‑1 (MQ=255)
          aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                           <  1:1428243/150‑1 (MQ=255)
           gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACCCAGTGTACTGCCCCGTTTGCCGATTGTCCTTGCACAATcggcg                                                                                                                                          >  2:1366322/1‑150 (MQ=255)
             cATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCGTTACACGCAATCAACGCAGTGTACTGCACCGTTTTCCGATTGTCCTTGCACAATCGGCggg                                                                                                                                        >  2:1452001/1‑150 (MQ=255)
                 tccTGAAAGCCCGCGGGAAACGGGCGCGGTTGAAACCGCGGCGCGCTTCGCCACCTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGAGTACTCCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaaa                                                                                                                                    <  2:879557/149‑1 (MQ=255)
                  ccTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCCTCATTTTGATCCCGTTACACGCAATCAACGCAGTGGACGGCACCGTTTGCCGTTTGTCCTTGCACCAGCGGCGGGtaaaa                                                                                                                                   >  1:421979/1‑150 (MQ=255)
                   cTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCCGCATCGTCGTCCTCATTTTCATCCCTTTACACGCAATCCACCCAGTGTTCTGGACCGTTTGCCGCTTGTCCGTTCACAATCGGCGGGGAAAat                                                                                                                                  >  1:1384927/1‑150 (MQ=255)
                          cACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAgg                                                                                                                           >  1:1194415/1‑150 (MQ=255)
                            cgcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTg                                                                                                                         >  2:1761031/1‑150 (MQ=255)
                            cgcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTg                                                                                                                         <  2:1007264/150‑1 (MQ=255)
                               gggAAACGGTCGCTTTTGAAACACCGGCGCGCTTCGCCACTTCCAGAATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                                                      <  1:862063/150‑1 (MQ=255)
                                     cGGTAGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCTATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTc                                                                                                                <  1:1228332/150‑1 (MQ=255)
                                       gTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcac                                                                                                              >  1:1296469/1‑150 (MQ=255)
                                        tCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcaca                                                                                                             <  1:562494/150‑1 (MQ=255)
                                        tCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCAGCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGGTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcaca                                                                                                             <  1:109636/150‑1 (MQ=255)
                                         cGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACaa                                                                                                            <  1:1464978/150‑1 (MQ=255)
                                          gCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACaaa                                                                                                           <  2:514540/150‑1 (MQ=255)
                                            ttttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAatt                                                                                                         >  1:458098/1‑149 (MQ=255)
                                            ttttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAata                                                                                                         <  1:508762/150‑1 (MQ=255)
                                              ttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAAtata                                                                                                       <  2:680341/150‑1 (MQ=255)
                                               tGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCTCCGTTTGCCGATTGTCCTTGCACAGTCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaa                                                                                                      >  1:327741/1‑150 (MQ=255)
                                                gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTTACCGGTTTCACAAATATaaa                                                                                                     >  2:423478/1‑150 (MQ=255)
                                                 aaaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaa                                                                                                    <  2:1319576/150‑1 (MQ=255)
                                                 aaaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaa                                                                                                    <  1:1674735/150‑1 (MQ=255)
                                                  aaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaaa                                                                                                   <  1:1490050/150‑1 (MQ=255)
                                                  aaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACCCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaaa                                                                                                   <  1:1375859/150‑1 (MQ=255)
                                                    ccccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAAt                                                                                                 >  2:485341/1‑150 (MQ=255)
                                                     cccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATg                                                                                                <  1:26645/150‑1 (MQ=255)
                                                     cccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATg                                                                                                <  2:1776499/150‑1 (MQ=255)
                                                     cccGGCGCGCTTCGCCACTTCCAGCATAGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTTCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATg                                                                                                >  1:174680/1‑150 (MQ=255)
                                                      ccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGa                                                                                               <  1:467519/150‑1 (MQ=255)
                                                      ccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGa                                                                                               >  1:815635/1‑150 (MQ=255)
                                                       cGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGaa                                                                                              <  1:1765300/150‑1 (MQ=255)
                                                        ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                                             <  1:255551/150‑1 (MQ=255)
                                                        ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                                             <  2:1332031/150‑1 (MQ=255)
                                                            cgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAtt                                                                                         >  2:669568/1‑150 (MQ=255)
                                                             gcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTc                                                                                        <  1:1629335/150‑1 (MQ=255)
                                                              cTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCa                                                                                       >  1:208405/1‑150 (MQ=255)
                                                                  gCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActct                                                                                   >  2:592192/1‑150 (MQ=255)
                                                                   ccACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctc                                                                                  >  2:1095420/1‑150 (MQ=255)
                                                                    cACTTCCATCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctct                                                                                 <  2:796201/150‑1 (MQ=255)
                                                                    cACTTCCAGCATCGTCGTCATCATTTTAATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTACACAATCGGCGGGAAAAAAATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctct                                                                                 >  2:1798780/1‑150 (MQ=255)
                                                                     aCTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCtt                                                                                <  1:1818901/150‑1 (MQ=255)
                                                                          cAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCCCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGCACCATTCCCTCGCTTGCTTa                                                                           >  2:340690/1‑150 (MQ=255)
                                                                                cgtcgtCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGt                                                                     <  2:1129428/150‑1 (MQ=255)
                                                                                   cgtcATCATTTTCATCCCTTTACACGCAATCCACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCAGAAATATAAAAAATGAACAATTCACTCTCTTGATTATTTAGTGAc                                                                  >  1:147453/1‑150 (MQ=255)
                                                                                   cgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGAc                                                                  >  1:1034192/1‑150 (MQ=255)
                                                                                      catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                                               <  2:900977/150‑1 (MQ=255)
                                                                                      catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTAATGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAAAATTCACTCTCTTGCTTATTAAGTGaaaaa                                                               >  1:1249686/1‑149 (MQ=255)
                                                                                          atattCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAACTGAACAATTCACTCTCTTGCTTATTTAGTGACAACTAtt                                                           <  1:791006/147‑1 (MQ=255)
                                                                                                 tCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtt                                                    >  2:1417204/1‑150 (MQ=255)
                                                                                                 tCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGAACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtt                                                    >  2:615134/1‑150 (MQ=255)
                                                                                                  cccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAttt                                                   >  1:288549/1‑150 (MQ=255)
                                                                                                  cccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAttt                                                   >  2:389151/1‑150 (MQ=255)
                                                                                                   ccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAATAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAattt                                                  >  2:263269/1‑150 (MQ=255)
                                                                                                   ccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtttt                                                  >  2:859308/1‑150 (MQ=255)
                                                                                                    cTTTACCCGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAATATTCATGATTTtg                                                 >  1:1377874/1‑150 (MQ=255)
                                                                                                      ttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                               <  2:1650167/150‑1 (MQ=255)
                                                                                                         cacGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaaa                                            >  2:759098/1‑150 (MQ=255)
                                                                                                           cGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAcc                                          <  2:1339994/150‑1 (MQ=255)
                                                                                                             cAATCAACGCAGAGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                        <  1:1111196/150‑1 (MQ=255)
                                                                                                              aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGt                                       >  1:731270/1‑150 (MQ=255)
                                                                                                                tCAACGCTGTGTACTGCACCGTTTGCCGATTGTCCTTGTACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGttt                                     <  2:1738709/150‑1 (MQ=255)
                                                                                                                tCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGttt                                     >  2:1400699/1‑150 (MQ=255)
                                                                                                                tCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGttt                                     >  2:508506/1‑150 (MQ=255)
                                                                                                                 cAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                    <  1:1095420/150‑1 (MQ=255)
                                                                                                                 cAACGCAGTGTACTGCACAGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                    <  1:368758/150‑1 (MQ=255)
                                                                                                                 aaaCGCACTTTAATGCACCGTTTGCCGATTGTCCTTGAACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                    <  2:1108539/149‑1 (MQ=255)
                                                                                                                   aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                                  <  2:183589/150‑1 (MQ=255)
                                                                                                                   aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                                  <  1:1048629/150‑1 (MQ=255)
                                                                                                                     gCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTaa                                >  1:1018640/1‑150 (MQ=255)
                                                                                                                         tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTcc                            <  1:1475575/150‑1 (MQ=255)
                                                                                                                           tACTGCACCGTTTGCCGATAGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGt                          <  1:95414/150‑1 (MQ=255)
                                                                                                                            aCTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGtt                         >  1:929951/1‑150 (MQ=255)
                                                                                                                            aCTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGtt                         <  1:717062/150‑1 (MQ=255)
                                                                                                                             cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGttt                        >  2:1639848/1‑150 (MQ=255)
                                                                                                                              tGCACCGTTTGCCTATTGTCCTTGCACAATCGGCGGTAAAAATATTCAGTTGACCGGTTTCACAAATATAAAAACTGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTc                       <  1:515833/150‑1 (MQ=255)
                                                                                                                              tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTc                       <  1:592192/150‑1 (MQ=255)
                                                                                                                              tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTATCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAATATTCATGATTTTGTGAAACATTTTTCTTAATTCCGttta                       >  1:331653/1‑149 (MQ=255)
                                                                                                                                cACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAg                     <  2:1724615/150‑1 (MQ=255)
                                                                                                                                 aCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGc                    >  1:505488/1‑150 (MQ=255)
                                                                                                                                   cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCAt                  >  2:1591732/1‑150 (MQ=255)
                                                                                                                                      ttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCgg               >  1:545432/1‑150 (MQ=255)
                                                                                                                                      ccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATATAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCAGTTTCAGCATCgg               >  1:796627/3‑150 (MQ=255)
                                                                                                                                      ccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCgg               >  1:1754172/3‑150 (MQ=255)
                                                                                                                                       cgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGc              >  1:1544957/2‑150 (MQ=255)
                                                                                                                                       cgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGc              >  2:345671/2‑150 (MQ=255)
                                                                                                                                        gCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCa             >  2:1379242/1‑150 (MQ=255)
                                                                                                                                         ccGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTACGTTTCAGCATCGGCAt            >  1:1715547/1‑150 (MQ=255)
                                                                                                                                          cGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCAtt           >  2:456179/1‑150 (MQ=255)
                                                                                                                                          cGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAACGGTTTCTTAATTCAGTTTCATCATCGGCAtt           >  1:226427/1‑150 (MQ=255)
                                                                                                                                            aTTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCATGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCAtttt         <  1:1115263/150‑1 (MQ=255)
                                                                                                                                                tCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCg     <  1:1551760/150‑1 (MQ=255)
                                                                                                                                                   ttGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCa  >  2:599859/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCA  >  NC_000913/2839176‑2839472

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: