BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,520,729	C→T	5.8%	W97* (TGG→TGA) ‡	hofP ←	DNA catabolic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,520,729	0	C	T	5.8%	324.2 / 5.3	102	W97* (TGG→TGA) ‡	hofP	DNA catabolic protein
Reads supporting (aligned to +/- strand): ref base C (53/43); new base T (2/4); total (55/47)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 104 total aligned reads displayed.

CAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGC  >  NC_000913/3520580‑3520876

cagccagcGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAAc                                                                                                                                                     <  1:1288116/150‑1 (MQ=255)

  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)

         ccTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGATATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGc                                                                                                                                            <  2:1528923/150‑1 (MQ=255)

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

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

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

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

              aaCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAATGGTTGTGGcaacgg                                                                                                                                       <  2:1388742/150‑6 (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:892500/150‑1 (MQ=255)

                             ggAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGt                                                                                                                        <  2:1474322/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:1531262/1‑150 (MQ=255)

                                                  gTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTgg                                                                                                   >  1:1408597/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                                                                        >  2:355036/1‑150 (MQ=255)

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

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

                                                                                  gggTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCa                                                                   <  2:912652/150‑1 (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:676001/1‑150 (MQ=255)

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

                                                                                          gCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCg                                                           >  2:86275/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:1568515/1‑150 (MQ=255)

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

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

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

                                                                                                                           tttGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCg                          <  1:1654259/150‑1 (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)

                                                                                                                                         ccGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCATCTGTAAAATTGTCCAGCCGTTTTCCAGCAAATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACAcc            <  1:220506/150‑1 (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)

CAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGCACCCGTCGCCATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGC  >  NC_000913/3520580‑3520876

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