BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,520,731	A→G	6.0%	W97R (TGG→CGG) ‡	hofP ←	DNA catabolic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,520,731	0	A	G	6.0%	289.0 / 6.2	100	W97R (TGG→CGG) ‡	hofP	DNA catabolic protein
Reads supporting (aligned to +/- strand): ref base A (49/44); new base G (2/4); total (52/48)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 103 total aligned reads displayed.

GCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGC  >  NC_000913/3520582‑3520876

gccagcGCCTGCACCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCa                                                                                                                                                   >  2:1229958/1‑150 (MQ=255)

gccagcGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCa                                                                                                                                                   >  1:1645627/1‑150 (MQ=255)

gccagcGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCa                                                                                                                                                   >  2:1534663/1‑150 (MQ=255)

 ccagcGCCTGCAGCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCAt                                                                                                                                                  <  1:567649/150‑1 (MQ=255)

    gcgcCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATtgt                                                                                                                                               <  1:1116519/150‑1 (MQ=255)

     cgcCTGCAACAGCTGAGCTACCGGAACGTCATCCTCCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTTGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATtgtg                                                                                                                                              <  1:1344671/150‑1 (MQ=255)

       ccTGCAACACCTGAGCTACTGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGc                                                                                                                                            <  1:1473466/150‑1 (MQ=255)

           cAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAATGGTTGTGGcaacg                                                                                                                                        <  2:1355091/150‑5 (MQ=255)

           cAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAATGGTTGTGGcaacg                                                                                                                                        <  2:1669553/150‑5 (MQ=255)

           cAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGtt                                                                                                                                        <  2:147125/150‑1 (MQ=255)

           cAACACCTGAGCGACCGGAACGTACTCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGtt                                                                                                                                        <  2:152532/150‑1 (MQ=255)

            aaCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAATGGTTGTGGcaacgg                                                                                                                                       <  2:1388742/150‑6 (MQ=255)

            aaCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAATGGTTGTGGcaacgg                                                                                                                                       <  1:1673247/150‑6 (MQ=255)

            aaCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTc                                                                                                                                       >  1:892500/1‑150 (MQ=255)

             acacCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCg                                                                                                                                      <  1:1601704/150‑1 (MQ=255)

              cacCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGc                                                                                                                                     >  1:1474322/1‑150 (MQ=255)

                 cTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGACGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATATCTATCTCCTTGCCGTTGACACAACCTATGAGGCGGTTCGCAGt                                                                                                                                  >  2:953123/1‑150 (MQ=255)

                  tGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGtt                                                                                                                                 <  1:1033356/150‑1 (MQ=255)

                        aCCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTccc                                                                                                                           >  1:724036/1‑150 (MQ=255)

                          cGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCgg                                                                                                                         <  2:503962/150‑1 (MQ=255)

                           ggcaCGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGt                                                                                                                        <  1:610321/147‑1 (MQ=255)

                           ggAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGt                                                                                                                        <  2:1474322/150‑1 (MQ=255)

                           ggAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGt                                                                                                                        <  2:892500/150‑1 (MQ=255)

                             aaCGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTAc                                                                                                                      >  1:1653003/1‑150 (MQ=255)

                               catCATCCACCCTCAGCGTCACTATTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTAccc                                                                                                                    <  2:219958/148‑1 (MQ=255)

                               cGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTAccc                                                                                                                    <  2:428585/150‑1 (MQ=255)

                                gTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGTGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGAGGCGGTTCGCAGTTTGTCCCGGTACCCa                                                                                                                   >  1:288845/1‑150 (MQ=255)

                                 tCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAg                                                                                                                  <  1:1534663/150‑1 (MQ=255)

                                    tCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGt                                                                                                               <  2:617856/150‑1 (MQ=255)

                                    tCCACCATCAGCGGCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGGTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGt                                                                                                               <  2:590516/150‑1 (MQ=255)

                                     ccaccaTCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGtt                                                                                                              >  1:800172/1‑150 (MQ=255)

                                        ccaTCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGGTTCATTTGTTTCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGGTTGTCCCGGTACCCAGCGTTAGt                                                                                                           <  1:45401/150‑1 (MQ=255)

                                         caTCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTa                                                                                                          <  2:1645627/150‑1 (MQ=255)

                                         caTCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTa                                                                                                          <  1:718453/150‑1 (MQ=255)

                                                gTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTgg                                                                                                   >  1:1100240/1‑150 (MQ=255)

                                                gTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTgg                                                                                                   >  1:1408597/1‑150 (MQ=255)

                                                gTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTgg                                                                                                   >  1:1531262/1‑150 (MQ=255)

                                                 tCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGc                                                                                                  >  1:1431878/1‑150 (MQ=255)

                                                    cTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCCCAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTc                                                                                               >  1:132044/1‑150 (MQ=255)

                                                    cTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCAAAACCATTGTGGCGGTTCGCAGTTTGTCCAGGTACCCAGCGTAAGTACGTCTGGCGTc                                                                                               >  2:159007/1‑150 (MQ=255)

                                                         tGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTg                                                                                          >  2:1644578/1‑150 (MQ=255)

                                                           cGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTa                                                                                        >  1:703369/1‑150 (MQ=255)

                                                             gCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTaaa                                                                                      >  2:267077/1‑150 (MQ=255)

                                                               tttGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAAt                                                                                    <  2:1419281/150‑1 (MQ=255)

                                                                 tGCCGCCTGTACGCCGGGTATCCGCATCAACAGTAGTGCGTCTATCCATTGCTTCATAAGAATCTCCTTGCCGTTGCCACAACCATTCTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGGTAGTACGTCTGGCGTCAGCTGTAAATTTg                                                                                  <  2:781326/150‑1 (MQ=255)

                                                                 tGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTg                                                                                  <  2:750342/150‑1 (MQ=255)

                                                                 tGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACATCTGGCGGCAGCTGTAAAATTg                                                                                  <  2:842931/150‑1 (MQ=255)

                                                                  gccgccTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTTTCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGt                                                                                 >  2:1602709/1‑150 (MQ=255)

                                                                    cgccTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTcc                                                                               <  1:1624859/150‑1 (MQ=255)

                                                                    cgccTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGCCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTcc                                                                               >  1:238188/1‑150 (MQ=255)

                                                                      ccTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAg                                                                             >  2:1178835/1‑150 (MQ=255)

                                                                        tGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGcc                                                                           >  1:1554254/1‑150 (MQ=255)

                                                                         gTACGCCGGGTATCAGCATCAACAGTATTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCAGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCg                                                                          <  1:953123/150‑1 (MQ=255)

                                                                         gTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTGCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCg                                                                          >  2:1502467/1‑150 (MQ=255)

                                                                           aCGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGtt                                                                        >  1:1364261/1‑150 (MQ=255)

                                                                           aCGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGtt                                                                        >  2:355036/1‑150 (MQ=255)

                                                                              ccGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCTTTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACACAGCGATAGTACGTCTGGCGTCAGCGGTAAAATTGTCCAGCCGtttta                                                                     >  2:132145/1‑149 (MQ=255)

                                                                                gggTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCa                                                                   <  2:912652/150‑1 (MQ=255)

                                                                                 ggTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAg                                                                  >  1:1567537/1‑150 (MQ=255)

                                                                                 ggTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAg                                                                  <  2:187826/150‑1 (MQ=255)

                                                                                  gTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGc                                                                 <  1:169207/150‑1 (MQ=255)

                                                                                    aTCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGcac                                                               >  1:1048228/1‑150 (MQ=255)

                                                                                     tCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGcaca                                                              >  2:1654259/1‑150 (MQ=255)

                                                                                      cAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTGCGCAGTTTGTCCCGGTACCCAGCGTTAGTCCGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGGTTTCCAGCACAt                                                             >  1:1165862/1‑150 (MQ=255)

                                                                                       aGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATc                                                            >  1:1467118/1‑150 (MQ=255)

                                                                                        gCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCg                                                           >  2:86275/1‑150 (MQ=255)

                                                                                        gCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCg                                                           >  1:1599447/1‑150 (MQ=255)

                                                                                        gCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCg                                                           >  2:676001/1‑150 (MQ=255)

                                                                                         cATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTATGGCGTCAGCTGTAAAATTGTCCAGCAGTTTTCCAGCACATCGt                                                          >  1:1551617/1‑150 (MQ=255)

                                                                                          aTCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTCCGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGtt                                                         >  2:858796/1‑150 (MQ=255)

                                                                                          aTCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGtt                                                         >  2:1162698/1‑150 (MQ=255)

                                                                                                   agtGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAACTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCAc                                                >  1:1105822/1‑150 (MQ=255)

                                                                                                    gtGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCAcc                                               <  1:8234/150‑1 (MQ=255)

                                                                                                      gCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCg                                             >  2:732842/1‑150 (MQ=255)

                                                                                                        ggCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTc                                           <  1:355036/150‑1 (MQ=255)

                                                                                                              ccATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCAtt                                     <  1:1367026/150‑1 (MQ=255)

                                                                                                                 ttGCTTCTTTTGTATCTCCTTGCAGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGGCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCt                                  <  2:288845/150‑1 (MQ=255)

                                                                                                                   gCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCttt                                >  1:1435875/1‑150 (MQ=255)

                                                                                                                     ttCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTg                              >  1:1534892/1‑150 (MQ=255)

                                                                                                                     ttCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTg                              >  1:1568515/1‑150 (MQ=255)

                                                                                                                       cATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGcc                            <  2:1100240/150‑1 (MQ=255)

                                                                                                                         tttGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCg                          <  1:1654259/150‑1 (MQ=255)

                                                                                                                         tttGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCg                          >  2:1395570/1‑150 (MQ=255)

                                                                                                                          ttGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGt                         <  1:267077/150‑1 (MQ=255)

                                                                                                                               tctcCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCtttt                    >  2:1096518/1‑150 (MQ=255)

                                                                                                                                ctcCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTtta                   <  1:86275/150‑1 (MQ=255)

                                                                                                                                  ccTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTttatt                 >  2:966253/1‑150 (MQ=255)

                                                                                                                                   cTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTttatta                <  1:1229958/150‑1 (MQ=255)

                                                                                                                                   cTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTttatta                <  2:1531262/150‑1 (MQ=255)

                                                                                                                                      gCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGCCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTacac             <  2:238188/150‑1 (MQ=255)

                                                                                                                                      gCCGTTGCCACAACCATTGTGGCGGTTCGCAGGTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTacac             >  2:1232280/1‑150 (MQ=255)

                                                                                                                                         accgCCACAATGGTTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCga          >  2:359945/4‑150 (MQ=255)

                                                                                                                                         accgCCACAATGGTTGTGGCGGTTCGCAGTTTGTCACGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACACGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCga          >  2:366184/4‑150 (MQ=255)

                                                                                                                                            gCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGcaca                                                              >  2:469780/1‑95 (MQ=255)

                                                                                                                                            gCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGcaca                                                              <  1:469780/95‑1 (MQ=255)

                                                                                                                                             ccACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCgatgaa      >  2:1244751/1‑149 (MQ=255)

                                                                                                                                              cacaACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGGCGCCATTTCTTTTGCCCGTCTTTTATTACACCGatgatg     >  2:950018/1‑150 (MQ=255)

                                                                                                                                                caACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGGTGCACCCGTCGCCATTTCTGTTGGCCGTCATTTATTACACCGATGATgcg   >  1:127021/1‑150 (MQ=255)

                                                                                                                                                 acccATTGGGGCGGTTCGCAGTTTCACCCGGTCCCCAGCGTTCGTACGTCTGGCGCCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGAACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCGATGATgcgc  <  1:818588/148‑1 (MQ=255)

                                                                                                                                                 aaCCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCGATGATgcgc  >  1:1208662/1‑150 (MQ=255)

GCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGC  >  NC_000913/3520582‑3520876

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