BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,410	A→T	7.3%	P56P (CCT→CCA)	yieL ←	putative hydrolase YieL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,410	0	A	T	7.3%	449.0 / 17.2	146	P56P (CCT→CCA)	yieL	putative hydrolase YieL
Reads supporting (aligned to +/- strand): ref base A (79/56); new base T (6/5); total (85/61)
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.40e-01

Only 100 of 150 total aligned reads displayed.

TACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGC  >  NC_000913/3900263‑3900559

tACACCATCAACTTTGAAATAATACTCCTACAGGTTGCCTTTCAGGATAGGTGTGCGCGACTACCAGAGTCCTGCTTCGTCTTTGGTCATCGGTTGAATATTGTCCGTAACCGGAACACCCACCACTACAGTGACATTTTTTGCGCCTgg                                                                                                                                                     <  2:1316222/150‑1 (MQ=255)

tACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgg                                                                                                                                                     <  1:552257/150‑1 (MQ=255)

 accccATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAggg                                                                                                                                                    >  1:442170/4‑150 (MQ=255)

 acacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAggg                                                                                                                                                    >  1:125768/1‑150 (MQ=255)

 acacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCAACGGGTGAATATTGTCCGCAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAggg                                                                                                                                                    >  1:1319182/1‑150 (MQ=255)

   acCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGc                                                                                                                                                  >  1:753319/1‑150 (MQ=255)

    ccATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGTAACCGGAACACCCACCACTACAGAGACATTTTTTGCACCAGGGGCa                                                                                                                                                 >  1:607500/1‑150 (MQ=255)

     cATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaa                                                                                                                                                >  1:561080/1‑150 (MQ=255)

       tCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCACCAGGGGCAAAg                                                                                                                                              >  1:951149/1‑150 (MQ=255)

        cAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGt                                                                                                                                             >  2:738106/1‑150 (MQ=255)

        cAACATAGAAAAAATACACCTACAGGATCCCTTTCACGATAGGTGTGCGCCACCACCAGACTCCACCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGt                                                                                                                                             >  2:20551/1‑150 (MQ=255)

           cATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaatg                                                                                                                                          <  2:389237/150‑6 (MQ=255)

            aTTGTAAAAATACTCGTACAGGTTGCCTTTCATGATTGGGGTGCTCCAGGTTCAGACTCCTGCTTCGTCTTTTGTAATGTTGTGAATATTGTCCTGAACCTTAACACCCACCACTACTTAGAAATTTTTTGCCCCAGGTTTAAAGTAgag                                                                                                                                         <  2:983373/150‑3 (MQ=255)

            aTTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTTGTCATCGGGTGAATATTGTCCTGAACCGGAACACCCACCACTACAGAGACATTTTTTGTCCCTGGGGCaaaaaatgt                                                                                                                                         <  2:108212/150‑7 (MQ=255)

               gAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGaaa                                                                                                                                      <  1:467998/150‑1 (MQ=255)

                 aaaaaTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACTACCAGACTCCTGCTTCGTCTTTTGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAgg                                                                                                                                    <  2:902287/150‑1 (MQ=255)

                 aaaaaTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaatg                                                                                                                                          <  1:1063151/144‑6 (MQ=255)

                 aaaaaTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAgg                                                                                                                                    <  2:1052844/150‑1 (MQ=255)

                       aCTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGACACGACCAGACTACTGCTTCGTCTTTGGTCATCGGGTGAATATTGTACGGAACCGGAAAACCCACCACTACAGAGAAATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACga                                                                                                                              >  2:808915/1‑149 (MQ=255)

                          cGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGt                                                                                                                           <  1:1423663/150‑1 (MQ=255)

                          cGTACAGGTTGCCTTTCAGGATAGGGGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGt                                                                                                                           <  1:874401/150‑1 (MQ=255)

                            tACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTa                                                                                                                         <  1:483531/150‑1 (MQ=255)

                                ggTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCgg                                                                                                                     >  2:1051716/1‑150 (MQ=255)

                                ggTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCGGTTATCgg                                                                                                                     >  2:606300/1‑150 (MQ=255)

                                ggTTGCCATTCAGGATAGGTGTGCGCCAAGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGCAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCgg                                                                                                                     >  2:427818/1‑150 (MQ=255)

                                 gTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGc                                                                                                                    >  1:835039/1‑150 (MQ=255)

                                  ttGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCa                                                                                                                   >  2:673639/1‑150 (MQ=255)

                                     ccTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTg                                                                                                                >  1:830270/1‑150 (MQ=255)

                                      cTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                                                                               >  1:603059/1‑150 (MQ=255)

                                         tCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCt                                                                                                            <  2:208514/150‑1 (MQ=255)

                                          cAGGATAGGTGGGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTg                                                                                                           <  1:1007121/150‑1 (MQ=255)

                                            ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                                                                         >  1:219212/1‑150 (MQ=255)

                                            ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                                                                         <  1:15950/150‑1 (MQ=255)

                                              aTAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTc                                                                                                       <  2:125768/150‑1 (MQ=255)

                                               tAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTca                                                                                                      <  2:942236/150‑1 (MQ=255)

                                                 ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTcaca                                                                                                    >  2:126512/1‑150 (MQ=255)

                                                  gtgtGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACAt                                                                                                   <  1:152253/150‑1 (MQ=255)

                                                     tgCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACt                                                                                                <  1:1292423/150‑1 (MQ=255)

                                                        gcCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGaaa                                                                                                                                      <  1:416792/109‑1 (MQ=255)

                                                                cAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCggg                                                                                     <  2:244910/150‑1 (MQ=255)

                                                                  gACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGAt                                                                                   >  1:523822/1‑150 (MQ=255)

                                                                  caCTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGAt                                                                                   <  1:616421/149‑1 (MQ=255)

                                                                   aCTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATa                                                                                  >  1:51523/1‑150 (MQ=255)

                                                                      ccTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTg                                                                               <  2:753319/150‑1 (MQ=255)

                                                                           ttCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGc                                                                          <  1:114531/150‑1 (MQ=255)

                                                                              gTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCgg                                                                       >  2:501664/1‑150 (MQ=255)

                                                                                cTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGGGGGTGCCGggc                                                                     <  2:84059/150‑1 (MQ=255)

                                                                                  ttGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctg                                                                   <  2:219212/150‑1 (MQ=255)

                                                                                      gcATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGATAGTTCACGCTGTTATCGTCATTGACCTTAGTCACATACTGTTTCACCGGTATAGTGGGTGCCTGgatggctg                                                               <  2:609897/149‑1 (MQ=255)

                                                                                       cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg                                                              >  1:488687/1‑150 (MQ=255)

                                                                                           gggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCata                                                          >  2:254966/1‑150 (MQ=255)

                                                                                            ggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCatat                                                         >  2:1083775/1‑150 (MQ=255)

                                                                                              tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc                                                       <  1:1158917/150‑1 (MQ=255)

                                                                                              tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc                                                       >  1:20882/1‑150 (MQ=255)

                                                                                               gAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCg                                                      >  1:249642/1‑150 (MQ=255)

                                                                                                 atatTGTCCGGAACCGGATCACCCACCACTACAGATACATTTTTTGCCACAGGGGCAAAGTAGCGAAACGTCACGCTGTTATAGGAAATGACCTGAGTCACATACTGTTTCACAGGGAAAGTGGCTGCCGGGCTGGCTGGCATATCCgct                                                    >  2:725966/1‑150 (MQ=255)

                                                                                                   atTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCgctgc                                                  >  1:323578/1‑150 (MQ=255)

                                                                                                    ttGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCgctgct                                                 >  1:713746/1‑150 (MQ=255)

                                                                                                       tCCGGAACCGGAACACCCACCACTACAGAGACATTCTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATg                                              >  1:1249443/1‑150 (MQ=255)

                                                                                                        ccggaaccggaaCACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATgg                                             >  1:233783/1‑150 (MQ=255)

                                                                                                         cggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGc                                            <  1:20551/150‑1 (MQ=255)

                                                                                                          ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGcc                                           >  1:1132869/1‑150 (MQ=255)

                                                                                                           gaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGTCTGGCATATCCGCTGCTATGGccc                                          <  2:592359/150‑1 (MQ=255)

                                                                                                              ccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCGCCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACt                                       >  1:964249/1‑150 (MQ=255)

                                                                                                              ccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACt                                       >  1:1060056/1‑150 (MQ=255)

                                                                                                              ccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACt                                       >  2:587579/1‑150 (MQ=255)

                                                                                                                ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGa                                     >  2:647327/1‑150 (MQ=255)

                                                                                                                ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGa                                     >  1:652962/1‑150 (MQ=255)

                                                                                                                  aacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGc                                   <  2:442170/150‑1 (MQ=255)

                                                                                                                  aacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGc                                   <  1:442445/150‑1 (MQ=255)

                                                                                                                   acacCCACCACTACAGAGACATTTTTTGCACCAGGGCCAAAATATCGAAATGTAACGCTGTAATAAGCATTGAACAGACTCACATACTTTATAAAAGTGAAAGTGGGTGCCGGGCTGCCTGGCATATACGCTGCTATGTCCCACTGAGCa                                  >  2:290964/1‑150 (MQ=255)

                                                                                                                     acccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGa                                >  2:1199857/1‑150 (MQ=255)

                                                                                                                      cccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaa                               >  2:1427042/1‑150 (MQ=255)

                                                                                                                        caccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAAt                             >  1:1015616/1‑150 (MQ=255)

                                                                                                                        caccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAAt                             >  1:791305/1‑150 (MQ=255)

                                                                                                                            acTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCg                         >  2:269081/1‑150 (MQ=255)

                                                                                                                               aCAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTg                      >  1:1196814/1‑150 (MQ=255)

                                                                                                                                  gagaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCa                   >  1:1095455/1‑150 (MQ=255)

                                                                                                                                  gagaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCa                   >  1:26549/1‑150 (MQ=255)

                                                                                                                                  gagaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCa                   >  2:380364/1‑150 (MQ=255)

                                                                                                                                  gagaCATTTTTTGCCACAGGAGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCa                   >  2:1164560/1‑150 (MQ=255)

                                                                                                                                      gctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg               >  1:1332063/6‑150 (MQ=255)

                                                                                                                                      cATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAg                                    <  2:489735/129‑1 (MQ=255)

                                                                                                                                       ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc              >  2:290110/5‑150 (MQ=255)

                                                                                                                                       ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc              >  1:1428397/5‑150 (MQ=255)

                                                                                                                                       ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc              >  1:288105/5‑150 (MQ=255)

                                                                                                                                        tactttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc             >  1:717975/4‑150 (MQ=255)

                                                                                                                                          ttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAg           >  2:488705/1‑150 (MQ=255)

                                                                                                                                            ttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCTCCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCGGGCAATAGCCAGCg         <  2:964249/150‑1 (MQ=255)

                                                                                                                                            ttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCg         <  1:1367726/150‑1 (MQ=255)

                                                                                                                                              gCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGtt       <  2:233783/150‑1 (MQ=255)

                                                                                                                                               ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCTTTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCGGTa      <  2:1319182/150‑1 (MQ=255)

                                                                                                                                               ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCGGGGCTGGCTGGCATATCCGCTGGTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCGGTa      <  2:51523/150‑1 (MQ=255)

                                                                                                                                               ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTa      >  1:386914/1‑150 (MQ=255)

                                                                                                                                               ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGGGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTa      <  1:1177494/150‑1 (MQ=255)

                                                                                                                                               ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGGTTCACCGGGATAGTGGGTGCCGggct                                                                    <  2:834769/88‑1 (MQ=255)

                                                                                                                                                 ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa    <  1:235075/150‑1 (MQ=255)

                                                                                                                                                 ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa    <  1:1303487/150‑1 (MQ=255)

                                                                                                                                                 ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATAACGCTGGCTATTGCCAGCGGTatt    <  1:476827/150‑3 (MQ=255)

                                                                                                                                                   aGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc  <  1:808915/150‑1 (MQ=255)

TACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGC  >  NC_000913/3900263‑3900559

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