Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,127 T→G 5.1% intergenic (+26/+31) bglA → / ← ygfF 6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,1270TG 5.1% 271.6 / 2.6 92intergenic (+26/+31)bglA/ygfF6‑phospho‑beta‑glucosidase A/putative NAD(P)‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (39/48);  new base G (2/3);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGT  >  NC_000913/3044979‑3045276
                                                                                                                                                    |                                                                                                                                                     
aaaCGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACtt                                                                                                                                                      <  2:1425097/150‑1 (MQ=255)
 aaCGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttg                                                                                                                                                     >  1:1378893/1‑150 (MQ=255)
   cGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgtt                                                                                                                                                   <  2:144781/150‑1 (MQ=255)
    gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATAGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                                  <  1:217950/150‑1 (MQ=255)
          gCTTTATCTATGTGAATAAACTTGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCg                                                                                                                                            <  2:415328/150‑1 (MQ=255)
          gCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCg                                                                                                                                            <  2:470905/150‑1 (MQ=255)
           cTTTATCTATGTTCATAAACATGACGACGGTACTTGCGATATGTCGCGTTCACGTAAGAAGAGCTATAACTGGTACAAAGAGGTGATCGCCAGCAACGGCGAGAAGCTTTAAGTCGATGTAGTACCGGATGCAATACTTGTTGCATCCgg                                                                                                                                           <  1:842200/150‑1 (MQ=255)
           cTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCgg                                                                                                                                           <  1:1332842/150‑1 (MQ=255)
                 cTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGAc                                                                                                                                     >  2:1485931/1‑150 (MQ=255)
                  tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACa                                                                                                                                    <  1:1610152/150‑1 (MQ=255)
                  tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACa                                                                                                                                    >  1:1322914/1‑150 (MQ=255)
                    tgtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAAc                                                                                                                                  >  2:1107267/1‑150 (MQ=255)
                      tgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACtt                                                                                                                                >  1:610470/1‑150 (MQ=255)
                      tgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCATTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACtt                                                                                                                                >  2:130563/1‑150 (MQ=255)
                       gAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTc                                                                                                                               <  2:608601/150‑1 (MQ=255)
                       gAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTc                                                                                                                               <  2:1101393/150‑1 (MQ=255)
                        aaTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTAGAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTTAGTCCATGAAGTACCGGATGCAATACTGGTTGCATCCGGGAAGACAACTTcc                                                                                                                              >  2:1084709/1‑150 (MQ=255)
                        aaTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                              <  2:1035737/150‑1 (MQ=255)
                         aTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTccc                                                                                                                             <  1:1219320/150‑1 (MQ=255)
                         aTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTccc                                                                                                                             <  1:314296/150‑1 (MQ=255)
                          tAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCt                                                                                                                            <  2:499482/150‑1 (MQ=255)
                           aaaCATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTAAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTg                                                                                                                           <  2:473950/150‑1 (MQ=255)
                            aaCATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGt                                                                                                                          <  2:381755/150‑1 (MQ=255)
                            aaCATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGt                                                                                                                          >  1:819293/1‑150 (MQ=255)
                             aCATGACGACGGTACTGGCGATATGTCGCGCTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGtt                                                                                                                         <  1:1555394/150‑1 (MQ=255)
                               aTGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGtttt                                                                                                                       <  2:804874/150‑1 (MQ=255)
                                tGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTa                                                                                                                      >  2:328502/1‑150 (MQ=255)
                                 gacgacGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGGGATTGCCAGCAACGTCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACATCCCTGTTTTAt                                                                                                                     <  1:1319382/150‑1 (MQ=255)
                                   cgacgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAttt                                                                                                                   <  1:1703122/150‑1 (MQ=255)
                                   cgacgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAttt                                                                                                                   <  2:918361/150‑1 (MQ=255)
                                     acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                                 <  2:842450/150‑1 (MQ=255)
                                          aCTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcccgccc                                                                                                            >  2:970418/1‑150 (MQ=255)
                                             ggCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCcgcccgcc                                                                                                         >  1:1064508/1‑150 (MQ=255)
                                             ggCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGACGATGAAGTACCGGATGCAATACTTGTTGCCTCCGGTCAGACAACTTGACTGTTTTATTTCCcgcccgcc                                                                                                         >  2:484677/1‑150 (MQ=255)
                                              gCGATATGTCGCGTTCACGTAAGAAGAGCTTAAACAGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGAAGAAGTACCGGATGCAAAACTTGTTGCATCCGGTCAGACAACTTCCCTTATTTATATGCCGCACGCCa                                                                                                        >  1:1072570/1‑150 (MQ=255)
                                                      tCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATa                                                                                                >  1:1412068/1‑150 (MQ=255)
                                                          gTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAAc                                                                                            <  2:1398428/150‑1 (MQ=255)
                                                               cGTAAAAAGCGCTTTACCTGGTACAAAGAGGTGATTGCCATCAACGGCGATAAGCTTTAAGTCGATGAAGTACCGGATGCACTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTccc                                                                                       <  2:230707/150‑1 (MQ=255)
                                                                   agaagaGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCACATCGATAAAACTTCCCGTGa                                                                                   >  1:1544478/1‑150 (MQ=255)
                                                                   agaagaGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGa                                                                                   >  2:665958/1‑150 (MQ=255)
                                                                      agagCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGt                                                                                >  2:416581/1‑150 (MQ=255)
                                                                      agagCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGt                                                                                >  1:424903/1‑150 (MQ=255)
                                                                         gCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAag                                                                             <  1:665958/150‑1 (MQ=255)
                                                                           tttAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAagag                                                                           <  2:455778/150‑1 (MQ=255)
                                                                             ttaCTGGTACAAAGAGGTGATTGCCAGCAAAGGCGAGTAGCTTTAAGTCGATGAAGTACCGGAGTCAATACTTGTTGCATCCGGTCAGACAAGTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGc                                                                         <  1:931206/148‑1 (MQ=255)
                                                                              aaCTGGTACAAAGAGGGGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCt                                                                        <  1:1260613/150‑1 (MQ=255)
                                                                               aCTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTAGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTACCTGTTTTATTTACCGCCCGCCAAATCGATAAAAATTCCCGTTACGTAAGAGGCtt                                                                       >  2:529892/1‑150 (MQ=255)
                                                                                cTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTTCATCAGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCttt                                                                      >  2:924175/1‑150 (MQ=255)
                                                                                  ggTACAAAGAGGTGATTTCCAGCAACGGCGAGAAGCTTTAAGTGGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTAt                                                                    <  2:1347668/150‑1 (MQ=255)
                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCa                                                                  >  1:791867/1‑150 (MQ=255)
                                                                                    tACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACGTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCCTGACGTAAGAGGCTTTATCa                                                                  >  1:1273887/1‑150 (MQ=255)
                                                                                      cAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACt                                                                >  1:1312019/1‑150 (MQ=255)
                                                                                       aaaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                               <  2:1064508/150‑1 (MQ=255)
                                                                                         agagGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtag                                                             <  2:610470/150‑1 (MQ=255)
                                                                                         agagGTGATTGACAGCAACGGCGAGAAGCTTTAAGTCGATGTAGTATCGGATGCAATACTTGTTTCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCTCCCGCCAAATCGATAAAACTTCCCGTGACGTTAGTGGCTTTATCACTtag                                                             <  1:104661/150‑1 (MQ=255)
                                                                                            ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCTGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGTCTTTATCACTtagtag                                                          <  1:886112/150‑1 (MQ=255)
                                                                                            ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCTGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtagtag                                                          >  2:314657/1‑150 (MQ=255)
                                                                                            ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtagtag                                                          <  2:497562/150‑1 (MQ=255)
                                                                                            ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtagtag                                                          >  1:293168/1‑150 (MQ=255)
                                                                                              cgATTGCCAGCAACGGCGAGAAGCTTTAACGCGATGAAGTACCGGATGCAAAACTTCTAGCATCCGGTCAGAAAACTTCCCTGTTTTATATCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTACGAGGCTTAATCACTTAGTAGcc                                                        <  2:1589609/149‑1 (MQ=255)
                                                                                                aTTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGCGGGGTTATCACTTAGTAGCCAg                                                      >  2:299082/1‑150 (MQ=255)
                                                                                                aTTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAg                                                      <  2:538038/150‑1 (MQ=255)
                                                                                                aTTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAg                                                      <  2:835275/150‑1 (MQ=255)
                                                                                                aTTGCAAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGAATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAg                                                      <  2:544846/150‑1 (MQ=255)
                                                                                                 ttGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGa                                                     >  2:1492575/1‑150 (MQ=255)
                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGCGGCTTTATCACTTAGTAGCCAGACTAta                                                >  1:121491/1‑149 (MQ=255)
                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                               >  2:238841/1‑135 (MQ=255)
                                                                                                      aGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                               <  1:238841/135‑1 (MQ=255)
                                                                                                         aaCGGCGAGAAGCTTTAAGTCGATGAAGTAGCGGTTGCACTACTTGTTGCATCCGGTCAGACATCTTCCCTGTTTTATATGCCTCCCGCCAAATCGATAAAACTTCCCGTGACGTTAGAGGCTTTATCACTTAGTAGCCAGACAATGGcc                                             <  2:1065025/150‑1 (MQ=255)
                                                                                                         aaCGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGcc                                             <  2:720809/150‑1 (MQ=255)
                                                                                                         aaCGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGcc                                             <  2:1104156/150‑1 (MQ=255)
                                                                                                           cGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTg                                           >  2:654500/1‑150 (MQ=255)
                                                                                                            ggCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTAATCACTTAGTAGCCAGACAATGGCCTgc                                          <  2:819293/150‑1 (MQ=255)
                                                                                                                   aGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACct                                   <  2:424903/150‑1 (MQ=255)
                                                                                                                   aGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCGTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACct                                   <  1:1050587/150‑1 (MQ=255)
                                                                                                                     cTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACctct                                 <  2:583752/150‑1 (MQ=255)
                                                                                                                        tAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGCCGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCtcttct                              >  2:1273094/1‑150 (MQ=255)
                                                                                                                         aaGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAGCTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGATGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGCCCtcttcgg                             >  1:1054213/1‑148 (MQ=255)
                                                                                                                         aaGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTg                             <  2:982811/150‑1 (MQ=255)
                                                                                                                          aGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGc                            >  2:773883/1‑150 (MQ=255)
                                                                                                                          aGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCACGTGACGTAAGAGGCTTTAACACTTAGTAGCCAGACAATGGCCTGCGCGACATCTTCTGc                            >  2:1638299/1‑150 (MQ=255)
                                                                                                                           ggcgcagAATTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                           <  1:885929/144‑1 (MQ=255)
                                                                                                                           gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                           <  1:745638/150‑1 (MQ=255)
                                                                                                                           gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                           <  1:598776/150‑1 (MQ=255)
                                                                                                                                tGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTcc                      <  2:1179401/150‑1 (MQ=255)
                                                                                                                                   aGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACATCCCTGTTTTATTTCCCGCCCGCCAAATCGATAACACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCGTCTGCCTGTccacc                   >  1:668100/1‑150 (MQ=255)
                                                                                                                                       ccGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCAGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc               <  1:1221764/150‑1 (MQ=255)
                                                                                                                                        cGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCt              >  2:813497/1‑150 (MQ=255)
                                                                                                                                          gATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGc            >  1:1038138/1‑150 (MQ=255)
                                                                                                                                           aTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCGAGACTATGGCGTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCa           >  1:1382959/1‑150 (MQ=255)
                                                                                                                                             gCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATg         >  1:600025/1‑150 (MQ=255)
                                                                                                                                              cAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                           >  1:1278263/1‑131 (MQ=255)
                                                                                                                                              cAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                           <  2:1278263/131‑1 (MQ=255)
                                                                                                                                               aaTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCAgggg       >  2:138732/1‑150 (MQ=255)
                                                                                                                                                  aCTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGAt    >  1:1051091/1‑150 (MQ=255)
                                                                                                                                                   cTTGTTGCATCCGGTCAGACAACTTCGCTGTTTTATATCCCGCCCGCCAAATCGTTAAAACTTCCCGTGACGTAAGTGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATg   <  1:317436/150‑1 (MQ=255)
                                                                                                                                                    ttgttgCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGATGCTTTATCACTTAGTAGCCAGTCATTGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGt  <  1:507507/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGATGT  >  NC_000913/3044979‑3045276

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: