BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,307,294	A→T	8.3%	L169Q (CTG→CAG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,294	0	A	T	8.3%	256.6 / 11.1	105	L169Q (CTG→CAG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand): ref base A (51/45); new base T (3/6); total (54/51)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 105 total aligned reads displayed.

GAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGG  >  NC_000913/4307146‑4307438

gAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCAATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                                                                                                 <  2:1525361/150‑1 (MQ=255)

   aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                              >  2:1697511/1‑150 (MQ=255)

     gggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                            >  2:1467300/1‑150 (MQ=255)

       gCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGata                                                                                                                                          <  1:796515/150‑1 (MQ=255)

       gCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATAGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGata                                                                                                                                          >  2:49700/1‑150 (MQ=255)

         gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTAATATCTCCCAGGGGGATATg                                                                                                                                        <  1:1191636/150‑1 (MQ=255)

           tttGCGCATGGACTAAAAGATCGCGCAATGGGTAATTTCGGGGCGGTTGTTCGTACCAGAGTCTTAGCGCCATTCCATAGCAATTGGTTTCAAGGCACCCAGGATTTCCACACGCGCAGTGTTGGGTCTTATCTCCCAGTGGGATATGAc                                                                                                                                      <  2:1375137/150‑1 (MQ=255)

           tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATGTCCCCCTGGGAGATATGAc                                                                                                                                      <  1:724112/150‑1 (MQ=255)

           tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                      <  1:194953/150‑1 (MQ=255)

           tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                      <  2:817263/150‑1 (MQ=255)

           tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACAAGCTTCTTAGCGCTATTCCAGAGCTATTGGTTTGCAGGCACCCTGGATTGCCACACGCGTAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                      <  1:106509/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTTCCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                     <  2:1160071/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                     <  2:123299/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                     <  2:39903/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                     <  2:1262582/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                     <  1:1368549/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                     <  1:1375345/150‑1 (MQ=255)

            ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCGGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                     <  1:584297/150‑1 (MQ=255)

             tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                    <  2:1032335/150‑1 (MQ=255)

             tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                    >  2:689188/1‑150 (MQ=255)

                   tGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTA‑CAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                             >  1:1541975/1‑150 (MQ=255)

                     gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                            >  1:1297502/1‑150 (MQ=255)

                      aCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCt                                                                                                                           >  2:698650/1‑150 (MQ=255)

                      aCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAAATGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcca                                                                                                                           >  1:282619/1‑149 (MQ=255)

                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                          >  2:1439336/1‑150 (MQ=255)

                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                          <  1:1187446/150‑1 (MQ=255)

                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCACCTg                                                                                                                          >  1:1381305/1‑150 (MQ=255)

                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCAAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCATGGGGATATGACCCAGTTCGCCTg                                                                                                                          >  2:1139634/1‑150 (MQ=255)

                        gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                                                                                                         <  2:1104945/150‑1 (MQ=255)

                         aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCTGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                        <  1:1697737/150‑1 (MQ=255)

                         aaaaGATCGCGCAATGGGTAATTTCGGGGCGGTTGGTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAAATGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                        <  2:1104794/150‑1 (MQ=255)

                             gATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCacac                                                                                                                    <  1:1467300/150‑1 (MQ=255)

                              aTCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACAcc                                                                                                                   >  1:1298237/1‑150 (MQ=255)

                                cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCCGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                 >  1:937525/1‑150 (MQ=255)

                                cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                 >  2:858156/1‑150 (MQ=255)

                                  cgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                                                                                               >  1:638137/1‑150 (MQ=255)

                                   gctaTGGGTAATTTGGGGGCTGATGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGtgtg                                                                                                              <  2:344971/147‑1 (MQ=255)

                                       tGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAcc                                                                                                          <  1:174972/150‑1 (MQ=255)

                                        gggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAccc                                                                                                         >  2:145312/1‑150 (MQ=255)

                                        gggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAccc                                                                                                         >  1:417377/1‑150 (MQ=255)

                                         ggTAATTTCGGGGCTGTTGTTCGTACCTGCGTGTTAGCGGCATTGCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCg                                                                                                        >  1:625025/1‑150 (MQ=255)

                                          gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGt                                                                                                       >  2:1235904/1‑150 (MQ=255)

                                          gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACACAGTTCGCCTGCCACACCGTGTGCACCCGt                                                                                                       >  2:1229884/1‑150 (MQ=255)

                                                 cGGGGCTGTTGTTCGTACCAGGGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCATGCACCCTGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCAAGGTGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgc                                                                                                <  2:1224355/150‑1 (MQ=255)

                                                 cGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgc                                                                                                >  1:1051971/1‑150 (MQ=255)

                                                 cGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgc                                                                                                >  2:363501/1‑150 (MQ=255)

                                                   gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                                                                                              <  2:191034/150‑1 (MQ=255)

                                                    ggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCa                                                                                             >  2:565159/1‑150 (MQ=255)

                                                         ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGtt                                                                                        <  1:698650/150‑1 (MQ=255)

                                                          tgttCGTACCAGCGTCTTAGCGCCATTCCAGATCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                       <  2:1648895/150‑1 (MQ=255)

                                                          tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                       >  1:2571/1‑150 (MQ=255)

                                                          tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                       >  1:593209/1‑150 (MQ=255)

                                                          tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTACCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACACGTCCACGGCGCACCGTTc                                                                                       >  1:935998/1‑150 (MQ=255)

                                                           gttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCA‑GGGGATATTACCCAGTTCGCCTGCCACCCCGTGTGCACCCGTCCCCGGCGCACCGGTCAt                                                                                     >  2:510874/1‑150 (MQ=255)

                                                            ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                     <  1:565159/150‑1 (MQ=255)

                                                             tCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATc                                                                                    <  2:1183314/150‑1 (MQ=255)

                                                              cGTACCATCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGTCACCCAGGCTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGTGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCAACGGCGCACCGTTCATcc                                                                                   <  2:1298546/150‑1 (MQ=255)

                                                                   cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCCGTTCGCCTGCCCCACCGTGTTCACCCGTCCACGGCCCACCGTTCATCCACACt                                                                              >  1:1085070/1‑150 (MQ=255)

                                                                    aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTg                                                                             >  1:1040366/1‑150 (MQ=255)

                                                                      cGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCg                                                                           <  2:1051971/150‑1 (MQ=255)

                                                                            aGCGCCATTCCAGAGCAATTGGTTTCCAGGCACACAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAACCGTCCACGGCGCACAGTTCATCCACACTGGGAAcccc                                                                     >  2:533382/1‑150 (MQ=255)

                                                                            aGCGCCATTCCAGAGCAATTGGTTTCCAGACACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAcccc                                                                     >  1:219499/1‑150 (MQ=255)

                                                                              cgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTAGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccccca                                                                   >  2:1317359/1‑148 (MQ=255)

                                                                               gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACTTTTTGTAAATGTTAAGGGGTCAACGTTGATCTACACTGCGAACCCCATc                                                                  <  2:1298237/150‑1 (MQ=255)

                                                                               gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCACCCCATc                                                                  >  2:1620155/1‑150 (MQ=255)

                                                                                   ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCCCACTGCGAACCCCATCCCCg                                                              >  1:1433997/1‑150 (MQ=255)

                                                                                    tCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTACGAACCCCATCCCCGt                                                             >  1:1606453/1‑150 (MQ=255)

                                                                                     ccAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa                                                            >  2:1652913/1‑150 (MQ=255)

                                                                                      ccgatctaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTACGAACCCCATCCCCGTAc                                                           <  2:1705765/143‑1 (MQ=255)

                                                                                      ccgatctaTTGGATTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGTGGGATATGACCCAGTTCGCCTGCCACACCTTGTGCACCCGTCCTCGGTGCACCGTTCATCCACACTACGAACCCCATCCCCGTAc                                                           <  2:357414/143‑1 (MQ=255)

                                                                                       agagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc                                                          >  2:803078/1‑150 (MQ=255)

                                                                                       agagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCAAATCTCCCAGGGGGATATGACCCATTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCCCACTGCGCACCCCATCCCCTTAcc                                                          >  2:1347761/1‑150 (MQ=255)

                                                                                          gCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgag                                                       <  1:489942/150‑1 (MQ=255)

                                                                                           cAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                      >  2:321652/1‑150 (MQ=255)

                                                                                           cAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGAGTTGGGTCATATCTCCCAGGGGGATATGACCCAGATCGCCAGACACACCGTGACCACCCATCCACGGCGCACCGTTCATACACACTACGAACCCCATCCCCGTACCcaga                                                      >  2:716964/1‑150 (MQ=255)

                                                                                                ggTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgc                                                 >  1:1386060/1‑150 (MQ=255)

                                                                                                  tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                               <  1:145312/150‑1 (MQ=255)

                                                                                                    tCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccgcc                                             <  2:625025/150‑1 (MQ=255)

                                                                                                      cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATTACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACAGCGCACCCCAACCCCGTACCGAGATAGGCCGCCAg                                           >  2:1259560/1‑150 (MQ=255)

                                                                                                       aGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGa                                          <  1:689188/150‑1 (MQ=255)

                                                                                                                       cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTaa                          >  1:47710/1‑150 (MQ=255)

                                                                                                                        acacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTACACCCGTCCACGGCGCACGGTTCATCCACACTGCGAACCCCTTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAg                         <  1:732741/150‑1 (MQ=255)

                                                                                                                             cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCCCCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGGAAGGCGGt                    >  1:868590/1‑150 (MQ=255)

                                                                                                                              ccAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtt                   >  1:1113834/2‑150 (MQ=255)

                                                                                                                                  tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCATATTGTGTGTAACTGTCCAGGGCGCAACGTTCATCAACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct               <  2:1297502/150‑1 (MQ=255)

                                                                                                                                  tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct               <  2:1040366/150‑1 (MQ=255)

                                                                                                                                   gtTGGGTCATCTCTCCCAGGTGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGTCGCACCGTTCATCCACACTTCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta              <  1:1413970/150‑1 (MQ=255)

                                                                                                                                     tGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTACGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact            <  2:1606453/150‑1 (MQ=255)

                                                                                                                                      gggTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGACCCCCATCCCCGTACCGAGATAGGCCGCCAGCACCAGTTGGTGCGTAAGGCGGTTTTatacta           >  2:125941/1‑150 (MQ=255)

                                                                                                                                       ggTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTCtactac          >  2:428666/1‑150 (MQ=255)

                                                                                                                                       ggTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTGGCGTAAGGCGGTTTTCtactac          >  2:1720378/1‑150 (MQ=255)

                                                                                                                                        gTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg         <  1:808292/150‑1 (MQ=255)

                                                                                                                                          cATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTc       <  1:1229884/150‑1 (MQ=255)

                                                                                                                                          cATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTc       <  1:1235904/150‑1 (MQ=255)

                                                                                                                                          cATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCACATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTc       >  1:421165/1‑150 (MQ=255)

                                                                                                                                            tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc     <  2:638137/150‑1 (MQ=255)

                                                                                                                                             atCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa    <  1:321652/150‑1 (MQ=255)

                                                                                                                                             atCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa    >  2:1208102/1‑150 (MQ=255)

                                                                                                                                             atCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGAGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGATCCCCATCCGCGTACCGAGATAGGCCGCGAGAACCAGTTGTTGCGTAAGGCGGTGTTCTACGACGACCCa    >  1:133022/1‑150 (MQ=255)

                                                                                                                                               ctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg  <  2:707332/150‑1 (MQ=255)

GAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGG  >  NC_000913/4307146‑4307438

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