BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	404,117	G→T	7.5%	G139C (GGC→TGC)	dgcC →	diguanylate cyclase DgcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	404,117	0	G	T	7.5%	443.3 / 14.6	133	G139C (GGC→TGC)	dgcC	diguanylate cyclase DgcC
Reads supporting (aligned to +/- strand): ref base G (56/67); new base T (4/6); total (60/73)
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.42e-01

Only 100 of 138 total aligned reads displayed.

AGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTT  >  NC_000913/403970‑404266

 gggCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGcc                                                                                                                                                    >  1:1669241/1‑150 (MQ=255)

  ggCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccc                                                                                                                                                   >  2:2152640/1‑150 (MQ=255)

    ccGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccccc                                                                                                                                                 >  1:1805694/1‑150 (MQ=255)

     cGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCg                                                                                                                                                <  2:860184/150‑1 (MQ=255)

     cGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCg                                                                                                                                                <  1:2451606/150‑1 (MQ=255)

      gTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGt                                                                                                                                               >  1:1572614/1‑150 (MQ=255)

       tCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATc                                                                                                                                              <  2:2478831/150‑1 (MQ=255)

       tCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTc                                                                                                                                              >  1:775313/1‑150 (MQ=255)

       tCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTc                                                                                                                                              <  1:2447780/150‑1 (MQ=255)

        cGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATca                                                                                                                                             <  1:2257019/150‑2 (MQ=255)

        cGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCt                                                                                                                                             <  1:1074309/150‑1 (MQ=255)

        cGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCt                                                                                                                                             >  2:1333605/1‑150 (MQ=255)

        cGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCt                                                                                                                                             <  2:389500/150‑1 (MQ=255)

        cGATCCGCGTAGCCGGGAAATTTACACCTTCACAACCGATGCAGTATTAGCGGGTATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCt                                                                                                                                             <  2:1678805/150‑1 (MQ=255)

         gATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTAGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTg                                                                                                                                            <  2:736878/150‑1 (MQ=255)

           tCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATcaaat                                                                                                                                          <  2:1989043/150‑5 (MQ=255)

              gCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATcaaattca                                                                                                                                       <  1:734913/150‑8 (MQ=255)

                 tAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCACGTCTGTTTGTcgc                                                                                                                                    >  1:1639835/1‑150 (MQ=255)

                  aGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTcgcg                                                                                                                                   <  1:2313417/150‑1 (MQ=255)

                        gAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTg                                                                                                                             >  1:861553/1‑150 (MQ=255)

                         aaaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTgg                                                                                                                            <  2:2437402/150‑1 (MQ=255)

                             ttACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGtt                                                                                                                        <  1:2012230/150‑1 (MQ=255)

                                cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                     >  1:2119407/1‑150 (MQ=255)

                                cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                     >  2:1571817/1‑150 (MQ=255)

                                 aaCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtg                                                                                                                    <  1:904536/150‑1 (MQ=255)

                                   cTTAAAAAACGATGCAGTATTAGAGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGACTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtctt                                                                                                                  >  2:1185749/1‑147 (MQ=255)

                                     tAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtggtgg                                                                                                                <  1:1571817/150‑1 (MQ=255)

                                       aaaaCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGtt                                                                                                              >  2:570620/1‑150 (MQ=255)

                                         aaCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                                            >  1:1360752/1‑150 (MQ=255)

                                          aCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTcc                                                                                                           <  1:1374447/150‑1 (MQ=255)

                                            cGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTg                                                                                                         >  1:2460972/1‑150 (MQ=255)

                                             gATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGc                                                                                                        >  2:1571833/1‑150 (MQ=255)

                                               tGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTGTGAATTTGATGGGGGCAGGCGGCCCCCGTGTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCt                                                                                                      <  1:1555788/150‑1 (MQ=255)

                                               tGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCTTGTCTGGTGTTTATGGTGGTTTCCTGCCt                                                                                                      <  2:2278638/150‑1 (MQ=255)

                                               tGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCt                                                                                                      <  1:1553184/150‑1 (MQ=255)

                                                 cAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                                    >  1:1545968/1‑150 (MQ=255)

                                                 cAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                                    <  1:307992/150‑1 (MQ=255)

                                                  aGTATTAGCGGGAATGTGGGTCGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGAGTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGGCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGt                                                                                                   <  2:799216/150‑1 (MQ=255)

                                                  aGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGGCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGt                                                                                                   <  2:2006047/150‑1 (MQ=255)

                                                    tATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCa                                                                                                 <  2:1901479/150‑1 (MQ=255)

                                                     aTTAGCGGGCATGTGGGTAGGCGTAATGGGCGGAACCGTGCGGCCTTCCACCGGGATGTGTATGATTATGGGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAc                                                                                                <  1:944327/150‑1 (MQ=255)

                                                        aGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGCTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCt                                                                                             >  2:1588370/1‑150 (MQ=255)

                                                         gCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTc                                                                                            >  2:1644922/1‑150 (MQ=255)

                                                           gggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGa                                                                                          >  2:1169131/1‑150 (MQ=255)

                                                           gggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGa                                                                                          <  1:1920552/150‑1 (MQ=255)

                                                            ggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCCCCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCTGGTGGAGGATGTGTGTTTGTCGCGGGTCTGGTGTTGATGTTTGTTTTGTGCGTTGTGACGCTCGAg                                                                                         <  2:2298001/150‑1 (MQ=255)

                                                             gAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGc                                                                                        >  2:1218114/1‑150 (MQ=255)

                                                             gAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGGCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGc                                                                                        <  1:514714/150‑1 (MQ=255)

                                                                  tgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACg                                                                                   >  2:2254351/1‑150 (MQ=255)

                                                                  tgGGTAGGCGTAATGGGCGT‑AACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACgg                                                                                  <  2:1761753/150‑1 (MQ=255)

                                                                    ggTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                                 <  1:859891/150‑1 (MQ=255)

                                                                     gTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGc                                                                                >  1:1187853/1‑150 (MQ=255)

                                                                      tAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCa                                                                               <  2:258948/150‑1 (MQ=255)

                                                                        ggCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCAtt                                                                             >  2:2019045/1‑150 (MQ=255)

                                                                          cGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTAc                                                                           <  1:804832/150‑1 (MQ=255)

                                                                              aTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCgtg                                                                       >  1:536253/1‑150 (MQ=255)

                                                                                gggCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTc                                                                     >  1:1213549/1‑150 (MQ=255)

                                                                                 ggCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATTAGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGATTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCg                                                                    >  2:961253/1‑150 (MQ=255)

                                                                                    gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                                 <  2:1536157/150‑1 (MQ=255)

                                                                                    gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                                 <  2:1226426/150‑1 (MQ=255)

                                                                                    gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGGTGATGGTGGGTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                                 <  2:1130847/150‑1 (MQ=255)

                                                                                       aaCGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAAt                                                              >  1:2483259/1‑150 (MQ=255)

                                                                                          gTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGt                                                           >  2:138989/1‑150 (MQ=255)

                                                                                            gctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTgc                                                         >  2:2108138/1‑150 (MQ=255)

                                                                                              tgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTgcgc                                                       <  1:2118604/150‑1 (MQ=255)

                                                                                                 cTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATGGTGcgccgc                                                    >  2:1799837/1‑150 (MQ=255)

                                                                                                 cTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccgc                                                    >  1:1359147/1‑150 (MQ=255)

                                                                                                    ccACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTgg                                                 >  2:1244142/1‑150 (MQ=255)

                                                                                                    ccACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTgg                                                 <  2:1816589/150‑1 (MQ=255)

                                                                                                      aCCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCTAGCTGACGGGCATTACAGTGTCGTTCAATAGTGCGCCTCTGGaa                                               <  2:1805694/150‑1 (MQ=255)

                                                                                                      aCCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGaa                                               >  1:2272923/1‑150 (MQ=255)

                                                                                                       ccGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                                              <  1:36267/150‑1 (MQ=255)

                                                                                                           gATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtg                                          <  1:2344045/150‑1 (MQ=255)

                                                                                                           gATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtg                                          <  2:2119407/150‑1 (MQ=255)

                                                                                                             tGTTGATGATTATGTGTGTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGTTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                                        <  1:1185749/150‑1 (MQ=255)

                                                                                                              gTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCa                                                                                                 <  1:2169787/92‑1 (MQ=255)

                                                                                                              gTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCa                                                                                                 >  2:2169787/1‑92 (MQ=255)

                                                                                                               ttGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                                      <  1:1868813/150‑1 (MQ=255)

                                                                                                                   tgatTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTccc                                  <  1:1039508/150‑1 (MQ=255)

                                                                                                                       tATGTGGCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTcc                              <  1:617755/150‑1 (MQ=255)

                                                                                                                        ttgtgTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTccc                             <  1:130139/149‑1 (MQ=255)

                                                                                                                        aTGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTCTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTccc                             <  1:338461/150‑1 (MQ=255)

                                                                                                                         tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCa                            >  2:904760/1‑150 (MQ=255)

                                                                                                                         tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCa                            <  2:1490544/150‑1 (MQ=255)

                                                                                                                            gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCatta                         >  1:176754/1‑150 (MQ=255)

                                                                                                                            gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCatta                         >  2:2057730/1‑150 (MQ=255)

                                                                                                                             tCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTGCCTTCCCattat                        <  1:1685650/150‑1 (MQ=255)

                                                                                                                              cTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattatt                       >  2:389676/1‑150 (MQ=255)

                                                                                                                               tGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattattt                      >  2:2345080/1‑149 (MQ=255)

                                                                                                                                  aTTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa                   <  1:2254351/150‑1 (MQ=255)

                                                                                                                                  aTTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa                   <  1:1333605/150‑1 (MQ=255)

                                                                                                                                      gacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa               >  2:1419693/4‑150 (MQ=255)

                                                                                                                                      gacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa               >  2:1343500/4‑150 (MQ=255)

                                                                                                                                      gATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa               <  2:300204/150‑1 (MQ=255)

                                                                                                                                      gATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCGGACGGGCATTACCGTGGCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa               <  2:642827/150‑1 (MQ=255)

                                                                                                                                              cAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctgc       >  2:863025/1‑150 (MQ=255)

                                                                                                                                              cAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctgc       <  2:1688732/150‑1 (MQ=255)

                                                                                                                                               tggcggcCCCCGTCTGTTTGTCGCGGTTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctgct      <  1:240216/149‑1 (MQ=255)

                                                                                                                                                ggcggcCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctgctg     <  2:2432461/150‑1 (MQ=255)

                                                                                                                                                   ggcCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGttt  >  1:1119467/1‑150 (MQ=255)

AGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTT  >  NC_000913/403970‑404266

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