Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,839,323 T→A 7.2% intergenic (‑59/‑201) 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,3230TA 7.2% 335.0 / 14.0 125intergenic (‑59/‑201)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 (53/62);  new base A (6/3);  total (60/65)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

Only 100 of 128 total aligned reads displayed.

CCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCA  >  NC_000913/2839174‑2839472
                                                                                                                                                     |                                                                                                                                                     
ccTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCCTCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCt                                                                                                                                                       >  1:1650167/1‑150 (MQ=255)
ccTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGAGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCt                                                                                                                                                       <  2:942937/150‑1 (MQ=255)
  tGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTg                                                                                                                                                     <  2:572234/150‑1 (MQ=255)
  gggCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTg                                                                                                                                                     <  1:1338903/149‑1 (MQ=255)
    gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                   <  1:793096/150‑1 (MQ=255)
    gCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGca                                                                                                                                                   <  1:1006965/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)
            aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                           <  1:890228/150‑1 (MQ=255)
             gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTAATCCCTTTACACGCAATAAACCCAGTGTACTGCACCGTTTGACGATTGTCCTTGCACAATaggcg                                                                                                                                          >  2:1373966/1‑150 (MQ=255)
              ccATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCAGCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggcaa                                                                                                                                         <  1:323125/150‑3 (MQ=255)
               cATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCGTTACACGCAATCAACGCAGTGTACTGCACCGTTTTCCGATTGTCCTTGCACAATCGGCggg                                                                                                                                        >  2:1452001/1‑150 (MQ=255)
                   tccTGAAAGCCCGCGGGAAACGGGCGCGGTTGAAACCGCGGCGCGCTTCGCCACCTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGAGTACTCCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaaa                                                                                                                                    <  2:879557/149‑1 (MQ=255)
                    cgtGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCAATTTCATCCCATTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCTTGCACAATCGGCGGGaaaaa                                                                                                                                   >  2:1551760/3‑150 (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:1007264/150‑1 (MQ=255)
                              cgcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTg                                                                                                                         >  2:1761031/1‑150 (MQ=255)
                                 gggAAACGGTCGCTTTTGAAACACCGGCGCGCTTCGCCACTTCCAGAATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                                                      <  1:862063/150‑1 (MQ=255)
                                       cGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTc                                                                                                                <  2:1770391/150‑1 (MQ=255)
                                       cGGTAGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCTATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTc                                                                                                                <  1:1228332/150‑1 (MQ=255)
                                         gTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcac                                                                                                              >  2:1674735/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)
                                              ttttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAata                                                                                                         <  1:508762/150‑1 (MQ=255)
                                              ttttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAata                                                                                                         <  2:1196616/150‑1 (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)
                                                       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                                                                                             <  2:1332031/150‑1 (MQ=255)
                                                              cgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAtt                                                                                         >  2:669568/1‑150 (MQ=255)
                                                              cgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACGGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAtt                                                                                         >  2:669574/1‑150 (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)
                                                                      cACTTCCAGCATCGTAGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctct                                                                                 <  2:448957/150‑1 (MQ=255)
                                                                           ccAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCtt                                                                            <  1:1518811/150‑1 (MQ=255)
                                                                            cAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCCCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGCACCATTCCCTCGCTTGCTTa                                                                           >  2:340690/1‑150 (MQ=255)
                                                                                   gtcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTg                                                                    <  2:864450/150‑1 (MQ=255)
                                                                                     cgtcATCATTTTCATCCCTTTACACGCAATCCACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCAGAAATATAAAAAATGAACAATTCACTCTCTTGATTATTTAGTGAc                                                                  >  1:147453/1‑150 (MQ=255)
                                                                                     cgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGAc                                                                  >  1:1034192/1‑150 (MQ=255)
                                                                                        catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                                               <  2:318309/150‑1 (MQ=255)
                                                                                        catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                                               <  2:900977/150‑1 (MQ=255)
                                                                                        catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTAATGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAAAATTCACTCTCTTGCTTATTAAGTGaaaaa                                                               >  1:1249686/1‑149 (MQ=255)
                                                                                            atattCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAACTGAACAATTCACTCTCTTGCTTATTTAGTGACAACTAtt                                                           <  1:791006/147‑1 (MQ=255)
                                                                                                   tCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGAACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtt                                                    >  2:615134/1‑150 (MQ=255)
                                                                                                    cccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAttt                                                   >  1:288549/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)
                                                                                                       tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgt                                                <  1:1688503/150‑1 (MQ=255)
                                                                                                        ttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                               <  1:1452001/150‑1 (MQ=255)
                                                                                                          acacGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaa                                             <  1:1761031/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)
                                                                                                                  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                                  <  1:1048629/150‑1 (MQ=255)
                                                                                                                       gCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTaa                                >  1:1018640/1‑150 (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)
                                                                                                                                    ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTAAGATTAAATGTTTAAAAAATATAAAAAATGAACAATTAACTCTATAAATTAATTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAgaa                   >  1:220406/1‑148 (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)
                                                                                                                                        ccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCgg               >  1:434194/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)
                                                                                                                                            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)
                                                                                                                                                     |                                                                                                                                                     
CCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCA  >  NC_000913/2839174‑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: