BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,625,067	A→C	5.3%	T99P (ACC→CCC)	ydeJ →	conserved protein YdeJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,625,067	0	A	C	5.3%	350.3 / 8.5	112	T99P (ACC→CCC)	ydeJ	conserved protein YdeJ
Reads supporting (aligned to +/- strand): ref base A (50/56); new base C (3/3); total (53/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 114 total aligned reads displayed.

CCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATT  >  NC_000913/1624918‑1625214

cccTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCaa                                                                                                                                                     >  1:78475/1‑150 (MQ=255)

 ccTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAAc                                                                                                                                                    >  2:179273/1‑150 (MQ=255)

  cTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAAcc                                                                                                                                                   <  1:241384/150‑1 (MQ=255)

   tgtgtgCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCCGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCg                                                                                                                                                  >  1:525747/1‑150 (MQ=255)

      gtgCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTg                                                                                                                                               <  1:924617/150‑1 (MQ=255)

      gtgCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTg                                                                                                                                               <  1:270241/150‑1 (MQ=255)

      gtgCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTg                                                                                                                                               >  2:1553057/1‑150 (MQ=255)

       tgCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGc                                                                                                                                              <  2:347619/150‑1 (MQ=255)

        gCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGcc                                                                                                                                             >  1:1054869/1‑150 (MQ=255)

           gCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATa                                                                                                                                          <  2:1500087/150‑1 (MQ=255)

              gAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTTAGTGAGAAAGTGGCAGCAGAAATGGCACCGGTTGCCAtttct                                                                                                                                       <  2:688116/150‑5 (MQ=255)

              gAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAAACTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATagag                                                                                                                                       >  1:844310/1‑150 (MQ=255)

               aaGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCACCGGTTGCCAtttctg                                                                                                                                      <  1:778579/150‑6 (MQ=255)

                aGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCg                                                                                                                                     <  2:425486/150‑1 (MQ=255)

                aGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCg                                                                                                                                     <  2:46672/150‑1 (MQ=255)

                   tACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGc                                                                                                                                  >  1:358692/1‑150 (MQ=255)

                     cacCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAGGCCAGCAATGGCTGGAACGATTTTGGGGGGTGAGTGAGAACTTGGCGGCGGAAATGGAACCCGGTGCCATGCAGCGTgggg                                                                                                                                >  2:210964/1‑150 (MQ=255)

                      acCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGa                                                                                                                               <  1:780317/150‑1 (MQ=255)

                       cccAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGAt                                                                                                                              >  2:1098908/1‑150 (MQ=255)

                       cccAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGAt                                                                                                                              >  2:15584/1‑150 (MQ=255)

                       cccAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGAt                                                                                                                              >  1:521408/1‑150 (MQ=255)

                       cccAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGAt                                                                                                                              >  2:205673/1‑150 (MQ=255)

                        ccAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATg                                                                                                                             >  2:1143465/1‑150 (MQ=255)

                        ccAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATg                                                                                                                             >  1:1175774/1‑150 (MQ=255)

                          aaaTTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCt                                                                                                                           <  2:1668298/150‑1 (MQ=255)

                             ttttACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGAt                                                                                                                        <  2:1656069/150‑1 (MQ=255)

                                  cGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAg                                                                                                                   <  2:1328226/150‑1 (MQ=255)

                                  cGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAg                                                                                                                   <  2:1501931/150‑1 (MQ=255)

                                      gCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTAtt                                                                                                               >  2:1046013/1‑150 (MQ=255)

                                           cTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGTCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCAt                                                                                                          >  1:303108/1‑150 (MQ=255)

                                             ttGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTa                                                                                                        <  1:1487549/150‑1 (MQ=255)

                                              tGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTAc                                                                                                       <  1:179273/150‑1 (MQ=255)

                                                 tACTTTCAGCGATCATGTATAGTTGAAAATCCTCTGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAATGTGGGAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTGAGTATTGCCATTACCgg                                                                                                    <  1:437657/150‑1 (MQ=255)

                                                 tACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCgg                                                                                                    <  1:1273511/150‑1 (MQ=255)

                                                  aCTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGc                                                                                                   <  1:1163561/150‑1 (MQ=255)

                                                   cTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCt                                                                                                  >  2:864941/1‑150 (MQ=255)

                                                   cTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCt                                                                                                  >  1:911036/1‑150 (MQ=255)

                                                    tttCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTa                                                                                                 >  1:123186/1‑150 (MQ=255)

                                                    tttCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTa                                                                                                 <  1:1143465/150‑1 (MQ=255)

                                                    tttCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTa                                                                                                 >  1:959291/1‑150 (MQ=255)

                                                     ttCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTAc                                                                                                >  1:35535/1‑150 (MQ=255)

                                                     ttCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTAc                                                                                                >  2:1243959/1‑150 (MQ=255)

                                                        aCCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTAcggc                                                                                             >  1:1391692/1‑150 (MQ=255)

                                                         ccGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTAcggcg                                                                                            <  1:1651901/150‑1 (MQ=255)

                                                               aGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCgg                                                                                      <  2:911036/150‑1 (MQ=255)

                                                                ggCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGa                                                                                     <  1:1098908/150‑1 (MQ=255)

                                                                ggCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGa                                                                                     <  1:1046013/150‑1 (MQ=255)

                                                                   aaaGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAggg                                                                                  <  1:563474/150‑1 (MQ=255)

                                                                    aaGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGc                                                                                 >  1:342669/1‑150 (MQ=255)

                                                                    aaGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGc                                                                                 >  1:425958/1‑150 (MQ=255)

                                                                      gATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCgg                                                                               <  1:1553057/150‑1 (MQ=255)

                                                                       aTGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGt                                                                              >  2:596508/1‑150 (MQ=255)

                                                                        tGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCgg                                                                                                                                                 <  2:965539/82‑1 (MQ=255)

                                                                         gAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCNACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGa                                                                            <  2:597021/150‑1 (MQ=255)

                                                                         gAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTAGTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGa                                                                            <  1:213360/150‑1 (MQ=255)

                                                                          aaaaTCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGaa                                                                           <  2:521408/150‑1 (MQ=255)

                                                                          aaaaTCCTCAGCGTAAGCCAGCAACCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGaa                                                                           <  2:844310/150‑1 (MQ=255)

                                                                              tCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATg                                                                       <  2:358692/150‑1 (MQ=255)

                                                                                cTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGt                                                                     >  2:1360546/1‑150 (MQ=255)

                                                                                cTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGt                                                                     <  2:499736/150‑1 (MQ=255)

                                                                                   aGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACg                                                                  <  2:1606180/150‑1 (MQ=255)

                                                                                    gCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGc                                                                 >  2:428057/1‑150 (MQ=255)

                                                                                    gCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGGGAAGATGGTACGc                                                                 >  1:693446/1‑150 (MQ=255)

                                                                                     cGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGcc                                                                <  2:105486/150‑1 (MQ=255)

                                                                                        aaGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGc                                                             >  1:597010/1‑150 (MQ=255)

                                                                                         aGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCg                                                            >  2:715153/1‑150 (MQ=255)

                                                                                         aGCCAGCAATCTCAGGAACGATATTCTGCGGTGAGTCAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTCCGCCcgcg                                                            >  2:247281/1‑150 (MQ=255)

                                                                                            cAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGt                                                         >  2:1347385/1‑150 (MQ=255)

                                                                                             aGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTa                                                        <  1:715153/150‑1 (MQ=255)

                                                                                              gCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTAc                                                       >  2:215666/1‑150 (MQ=255)

                                                                                               cAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGGGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTAcc                                                      <  1:250965/150‑1 (MQ=255)

                                                                                                  tctcTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTc                                                   <  2:967180/150‑1 (MQ=255)

                                                                                                     ctTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTgg                                                <  2:548746/150‑1 (MQ=255)

                                                                                                      ttGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGt                                               >  1:1254080/1‑150 (MQ=255)

                                                                                                          aCGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGc                                           >  1:494856/1‑150 (MQ=255)

                                                                                                          aCGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGc                                           >  1:850206/1‑150 (MQ=255)

                                                                                                            gATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCTGAGTGCGGTTAAGATTGTACGACAGAGGGTACCGTCTGGTTTGCGt                                         <  2:959291/150‑1 (MQ=255)

                                                                                                            gATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGAGGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGt                                         <  1:1126700/150‑1 (MQ=255)

                                                                                                             atatTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTg                                        <  2:461894/150‑1 (MQ=255)

                                                                                                               atTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGc                                      <  2:1620199/150‑1 (MQ=255)

                                                                                                               atTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGc                                      <  1:205673/150‑1 (MQ=255)

                                                                                                                 tCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCat                                    >  2:1126111/1‑150 (MQ=255)

                                                                                                                    gCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATAtt                                 >  2:1579848/1‑150 (MQ=255)

                                                                                                                     cGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTa                                >  2:1582129/1‑150 (MQ=255)

                                                                                                                          agtgagAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGc                           <  1:15584/150‑1 (MQ=255)

                                                                                                                           gtgagAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGcc                          >  1:878251/1‑150 (MQ=255)

                                                                                                                              agaAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTATGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGa                       >  1:1354236/1‑150 (MQ=255)

                                                                                                                              agaAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGa                       <  1:503790/150‑1 (MQ=255)

                                                                                                                                aaaGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAAc                     <  1:864941/150‑1 (MQ=255)

                                                                                                                                aaaGGGGCAGCAGAAATGGCAACCGGCGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAAc                     <  2:895221/150‑1 (MQ=255)

                                                                                                                                 aaGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACt                    <  2:1206234/150‑1 (MQ=255)

                                                                                                                                          cagAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCgg           >  2:558550/1‑150 (MQ=255)

                                                                                                                                           ctctaTGGCACCGGTTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAAc                     <  2:204353/135‑1 (MQ=255)

                                                                                                                                           ctctaTGGCACCGGTTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAAc                     >  1:204353/5‑139 (MQ=255)

                                                                                                                                           ctctaTGGCACCGGTTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGt          >  2:1016325/5‑150 (MQ=255)

                                                                                                                                            gAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGtt         <  1:498683/150‑1 (MQ=255)

                                                                                                                                             aaaTGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTa        <  1:215666/150‑1 (MQ=255)

                                                                                                                                                tGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGt          <  1:1650320/145‑1 (MQ=255)

                                                                                                                                                tGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGt          >  2:1650320/1‑145 (MQ=255)

                                                                                                                                                   cAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCAtt  >  2:388407/1‑150 (MQ=255)

CCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGAGAAAGTGGCAGCAGAAATGGCAACCGGTGCCATAGAGCGTGCGGATGCTGATGTCAGTATTGCCATTACCGGCTACGGCGGACCGGAGGGCGGTGAAGATGGTACGCCAGCGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATT  >  NC_000913/1624918‑1625214

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: `‑`