Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,124 T→C 6.9% intergenic (+23/+34) bglA → / ← ygfF 6‑phospho‑beta‑glucosidase A/putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,1240TC 6.9% 255.7 / 7.6 87intergenic (+23/+34)bglA/ygfF6‑phospho‑beta‑glucosidase A/putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base C (3/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGA  >  NC_000913/3044975‑3045273
                                                                                                                                                     |                                                                                                                                                     
cAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAAt                                                                                                                                                       >  1:89331/1‑150 (MQ=255)
cAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAAt                                                                                                                                                       <  2:1514077/150‑1 (MQ=255)
  gCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATAc                                                                                                                                                     >  1:623185/1‑150 (MQ=255)
       cGCCACGGCTTTATCTATGTGAATAAACATGACTACGGTACTCGCGTTATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgtt                                                                                                                                                <  2:785334/150‑1 (MQ=255)
        gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTGCAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                               <  1:178409/150‑1 (MQ=255)
        gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                               <  2:1782883/150‑1 (MQ=255)
        gCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACttgttg                                                                                                                                               <  1:747502/150‑1 (MQ=255)
         cTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGc                                                                                                                                              >  2:357977/1‑150 (MQ=255)
             ggCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATcc                                                                                                                                          <  2:586191/150‑1 (MQ=255)
             ggCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATcc                                                                                                                                          <  2:203722/150‑1 (MQ=255)
              gCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCg                                                                                                                                         <  1:149537/150‑1 (MQ=255)
               cTTTATCTATGTGAATAAATATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCgg                                                                                                                                        <  2:210372/150‑1 (MQ=255)
               cTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAACAAGTATTGCATCCgg                                                                                                                                        <  2:400313/150‑1 (MQ=255)
                   aTCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAg                                                                                                                                    >  2:1009351/1‑150 (MQ=255)
                     cTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGAc                                                                                                                                  <  1:357977/150‑1 (MQ=255)
                     cTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGAc                                                                                                                                  >  1:485051/1‑150 (MQ=255)
                      tATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACa                                                                                                                                 >  1:68813/1‑150 (MQ=255)
                       aTGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACaa                                                                                                                                >  1:596853/1‑150 (MQ=255)
                        tgtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAAc                                                                                                                               >  2:1715703/1‑150 (MQ=255)
                         gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACt                                                                                                                              >  1:214607/1‑150 (MQ=255)
                            aaTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                           >  2:1046140/1‑150 (MQ=255)
                            aaTAAACATGACGACGGTACTGGCGAGATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTcc                                                                                                                           <  1:1843266/150‑1 (MQ=255)
                             aTAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTccc                                                                                                                          <  1:93632/150‑1 (MQ=255)
                                aaCATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGt                                                                                                                       >  1:1499600/1‑150 (MQ=255)
                                    tGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTa                                                                                                                   <  1:702539/150‑1 (MQ=255)
                                       cgacgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTAttt                                                                                                                <  2:618784/150‑1 (MQ=255)
                                         acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              <  1:1542623/150‑1 (MQ=255)
                                         acgGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcc                                                                                                              <  1:1656922/150‑1 (MQ=255)
                                            gTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcccgc                                                                                                           <  2:1403488/150‑1 (MQ=255)
                                               cTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcccgcccg                                                                                                        >  1:676640/1‑150 (MQ=255)
                                               cTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATCAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTcccgcccg                                                                                                        <  1:161698/150‑1 (MQ=255)
                                                tGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCccgcccgc                                                                                                       >  1:187798/1‑150 (MQ=255)
                                                 ggCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCcgcccgcc                                                                                                      <  1:1766283/150‑1 (MQ=255)
                                                        tGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCgt                                                                                               >  2:116447/1‑149 (MQ=255)
                                                         gTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGAt                                                                                              >  1:273229/1‑150 (MQ=255)
                                                           cgcgTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATaa                                                                                            >  1:190675/1‑150 (MQ=255)
                                                              gTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAAc                                                                                         <  1:1715703/150‑1 (MQ=255)
                                                              gTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAAc                                                                                         <  1:1009351/150‑1 (MQ=255)
                                                                  aCGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTcc                                                                                     <  2:214607/150‑1 (MQ=255)
                                                                    gTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCg                                                                                   <  2:485051/150‑1 (MQ=255)
                                                                    gTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCg                                                                                   <  2:273229/150‑1 (MQ=255)
                                                                      aagaagAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTg                                                                                 <  2:676640/150‑1 (MQ=255)
                                                                      aagaagAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTg                                                                                 <  2:68813/150‑1 (MQ=255)
                                                                       agaagaGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGAAATACATCTAGCATCCCCTCAGACAACTTCCATGTTTTATTTCCCGACCGCCAAATCCATAAAACTTCCCGTGa                                                                                >  2:8903/1‑150 (MQ=255)
                                                                          agagCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGt                                                                             >  1:70418/1‑150 (MQ=255)
                                                                          agagCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGt                                                                             >  2:24494/1‑150 (MQ=255)
                                                                                      ggTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTAt                                                                 <  1:1046140/150‑1 (MQ=255)
                                                                                      ggTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCACGCACGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTAt                                                                 >  2:617646/1‑150 (MQ=255)
                                                                                         aCAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCAc                                                              >  2:227072/1‑150 (MQ=255)
                                                                                           aaaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACtt                                                            <  2:56289/150‑1 (MQ=255)
                                                                                            aaGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTta                                                           >  2:1919841/1‑150 (MQ=255)
                                                                                              cagGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGAAGAAGTACCGGATGCAATACATGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtagt                                                         <  1:380008/149‑1 (MQ=255)
                                                                                                ggTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTCTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTtagtag                                                       >  2:464443/1‑150 (MQ=255)
                                                                                                  tGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGcc                                                     <  2:65647/150‑1 (MQ=255)
                                                                                                       gCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACATCTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTTGCCAGACa                                                <  2:623185/150‑1 (MQ=255)
                                                                                                        ccAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACaa                                               >  2:977451/1‑150 (MQ=255)
                                                                                                           gCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATgg                                            >  1:113682/1‑150 (MQ=255)
                                                                                                            cAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGc                                           <  1:24494/150‑1 (MQ=255)
                                                                                                             aaCGGCGAGAATCTTTAAGTGGATGAAGTACCGGATGGAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGTTAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGcc                                          <  2:89331/150‑1 (MQ=255)
                                                                                                                  cGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTgcgc                                     >  1:873217/1‑150 (MQ=255)
                                                                                                                  cGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTACCACTTAGTAGCCAGACAATGGCCTgcgc                                     >  2:1392422/1‑150 (MQ=255)
                                                                                                                   gagaAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGcgcg                                    <  2:187798/150‑1 (MQ=255)
                                                                                                                        gCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACctc                               >  1:439469/1‑150 (MQ=255)
                                                                                                                         cTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACctct                              >  1:516086/1‑150 (MQ=255)
                                                                                                                          tttAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGc                                                      <  2:815171/125‑1 (MQ=255)
                                                                                                                          tttAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGc                                                      >  1:815171/1‑125 (MQ=255)
                                                                                                                          tttAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCtctt                             >  1:1905817/1‑150 (MQ=255)
                                                                                                                           ttAAGTCTATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATTAAACTTGCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGAGAATTGCCTGATCGACCtcttc                            <  2:1221079/150‑1 (MQ=255)
                                                                                                                             aaGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTg                          >  2:497263/1‑150 (MQ=255)
                                                                                                                               gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGGCAAATCGTTATAAGTTCGCGTGACGTAAGAGGCTTTATCACTTAGTAGCGAGACATTGGCCTGCGCGACCTCTTCTGcc                        <  2:1499600/150‑1 (MQ=255)
                                                                                                                               gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                        >  2:907514/1‑150 (MQ=255)
                                                                                                                               gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                        >  1:514534/1‑150 (MQ=255)
                                                                                                                               gTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTCTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                        <  1:464443/150‑1 (MQ=255)
                                                                                                                                 cGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGcc                        <  1:1447482/148‑1 (MQ=255)
                                                                                                                                 cGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTg                      >  2:460868/1‑150 (MQ=255)
                                                                                                                                 cGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTg                      <  2:1658775/150‑1 (MQ=255)
                                                                                                                                 cGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCAAGACAATGGCCTGCGCGACCTCTTCTGcc                        >  2:1447482/1‑148 (MQ=255)
                                                                                                                                     gAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTcca                  >  2:1886749/1‑150 (MQ=255)
                                                                                                                                       aGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTccacc                >  2:776980/1‑150 (MQ=255)
                                                                                                                                         tACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCcaccac              <  1:1714800/150‑1 (MQ=255)
                                                                                                                                           ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc            >  1:441313/1‑150 (MQ=255)
                                                                                                                                           ccGGATGCAACAAGTATTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGc            >  1:609381/1‑150 (MQ=255)
                                                                                                                                              gATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGc         <  1:227072/150‑1 (MQ=255)
                                                                                                                                              gATGCAACAAGTATTGCATCTGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGc         >  2:787128/1‑150 (MQ=255)
                                                                                                                                               aTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTTTCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCa        <  2:1304077/150‑1 (MQ=255)
                                                                                                                                               aTGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCa        <  2:190675/150‑1 (MQ=255)
                                                                                                                                                 gCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATg      >  2:558486/1‑150 (MQ=255)
                                                                                                                                                    aTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCa        <  2:803802/145‑1 (MQ=255)
                                                                                                                                                    aTACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCa        >  1:803802/1‑145 (MQ=255)
                                                                                                                                                     tACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGa  <  2:147802/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCAAACGCTACGGCTTTATCTATGTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTATTTCCCGCCCGCCAAATCGATAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGCATGGGGA  >  NC_000913/3044975‑3045273

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