BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,251,392	C→A	10.3%	A123A (GCC→GCA)	yqjG →	putative S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,251,392	0	C	A	10.3%	304.0 / 9.7	123	A123A (GCC→GCA)	yqjG	putative S‑transferase
Reads supporting (aligned to +/- strand): ref base C (60/48); new base A (2/13); total (62/61)
Fisher's exact test for biased strand distribution p-value = 2.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCC  >  NC_000913/3251246‑3251541

cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGAt                                                                                                                                                    <  1:1100729/150‑1 (MQ=255)

cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGAt                                                                                                                                                    <  2:1380977/150‑1 (MQ=255)

atCAAAGGACTGTAACCGTTTATTTCCGTTTCCGTAGTGAACCAGCTGATGCTGGAAAACGGTTGGACATTTGATGACAGTTTTCCGGTAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGAt                                                                                                                                                    <  1:1008456/149‑1 (MQ=255)

 tCATAGGACTTGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATc                                                                                                                                                   <  1:477973/150‑1 (MQ=255)

 tCAAATGACTGGAACCGTTTATTTCCGTATCCTTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATc                                                                                                                                                   <  1:783408/150‑1 (MQ=255)

 tCAAAGGACTGGCACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATc                                                                                                                                                   <  2:1768676/150‑1 (MQ=255)

 tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATc                                                                                                                                                   >  2:1762685/1‑150 (MQ=255)

  cAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATcc                                                                                                                                                  <  1:920841/150‑1 (MQ=255)

  aaaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATcc                                                                                                                                                  <  2:1663761/149‑1 (MQ=255)

    aaGGACTGGAACCGTTTCTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACTCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCcac                                                                                                                                                <  2:1569432/150‑1 (MQ=255)

    aaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCcac                                                                                                                                                <  1:1952743/150‑1 (MQ=255)

     aGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCcaca                                                                                                                                               >  2:1147947/1‑150 (MQ=255)

     aGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCcaca                                                                                                                                               <  1:364717/150‑1 (MQ=255)

           ggAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTTTCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAg                                                                                                                                         >  1:1269364/1‑150 (MQ=255)

                cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATTTTTTTTTTTATTTTTTTTTTGACCACGGAGATCTACACGACAGCGGAc                                                                                                                                    <  2:649250/150‑1 (MQ=255)

                cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGAc                                                                                                                                    <  1:577317/150‑1 (MQ=255)

                cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGG‑AAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACg                                                                                                                                   >  2:978476/1‑150 (MQ=255)

                   ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTTGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTAAAGCGGACGAg                                                                                                                                 >  2:1813987/1‑150 (MQ=255)

                   ttATTTCCGGTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATTTATTTTTTTATTTTTTTTTTGTCCACGGAGATCTACACTACAGCGGACGAg                                                                                                                                 <  2:596309/150‑1 (MQ=255)

                    tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGt                                                                                                                                >  1:1068521/1‑150 (MQ=255)

                     aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGCACGAGtt                                                                                                                               >  1:886488/1‑150 (MQ=255)

                      tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTa                                                                                                                              >  1:613307/1‑150 (MQ=255)

                      tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTa                                                                                                                              >  2:1475940/1‑150 (MQ=255)

                         ccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTg                                                                                                                           >  2:1167772/1‑150 (MQ=255)

                          cgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTCTATCTCCACGCCGATCCCCACTACAGCGGACGAGTTACTGt                                                                                                                          >  1:756377/1‑150 (MQ=255)

                              tccgtCGTGCACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACACTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATTAATTTATGTATATTTTTTTTGTCCACGGAGATCTACACTACAGCGGACGAGTTACTGTTccc                                                                                                                      <  2:369254/150‑1 (MQ=255)

                               ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCg                                                                                                                     <  1:1963674/150‑1 (MQ=255)

                               ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCg                                                                                                                     >  2:1541149/1‑150 (MQ=255)

                               ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCATTATATCAACTCCGATCCACACTACAGCGGACGAGTTACTGTTCCCg                                                                                                                     >  2:1006937/1‑150 (MQ=255)

                                 gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTg                                                                                                                   <  1:867671/150‑1 (MQ=255)

                                  tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATTATTTTTTTTATTTTTTTTATGTCCACGGAGATCAACACTACAGCGGACGAGTTACTGTTCCCGTGc                                                                                                                  <  2:81719/150‑1 (MQ=255)

                                  tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGc                                                                                                                  <  1:1271547/150‑1 (MQ=255)

                                   aGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCt                                                                                                                 >  2:550085/1‑150 (MQ=255)

                                    gTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCtg                                                                                                                >  2:769507/1‑150 (MQ=255)

                                      gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCtgtg                                                                                                              >  2:2050212/1‑150 (MQ=255)

                                       aaCCCGCTGATGCTGGATAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATTATTTTTTGTATCTGTTTTTTGTCCACGCAGATGTACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTgg                                                                                                             <  2:158152/150‑1 (MQ=255)

                                             cTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACaaa                                                                                                       <  1:686450/150‑1 (MQ=255)

                                             cTGATGCTGGAAAACGGTTGGACCTTAGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCTTGCTGTGGGACaaa                                                                                                       >  2:1944031/1‑150 (MQ=255)

                                              tGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGTACGAGTTACTGTTCCCGTGATGTGGGACaaaa                                                                                                      >  2:347676/1‑150 (MQ=255)

                                                 tGCTGGAAAACGGTTGCACCAATGATGACACTTTTCCGCGAGCACCCGGAGACACGCACTATCAAAATGAATTTCTGTATAAGCTTTATCTCCACGCCGATCCACACAACAGCGGACGAGTTACTGTTCCCGTGCTCTGGGAAAAAAAGa                                                                                                   >  2:2088962/1‑150 (MQ=255)

                                                   cTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCAGATCCACAAAACAGCTGACGAGTTACTGTTCCCGTGCTGTGGGGCAAAAAGAAc                                                                                                 >  1:1459058/1‑150 (MQ=255)

                                                    tGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAAcc                                                                                                >  1:1812581/1‑150 (MQ=255)

                                                    tGGAAAACGGTTGGACCTTTGATGACAGGTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAAcc                                                                                                <  1:51244/150‑1 (MQ=255)

                                                      gAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATATCAACGCCGATCCACACTACAGCGGACGAGTTAATGTTCCCGTGCTGTGGGACAAAAAGAACcac                                                                                              >  1:1052897/1‑150 (MQ=255)

                                                       aaaaCGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTTTATCATCTTTATCTCCACGCAGATCTACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACcaca                                                                                             <  2:211688/150‑1 (MQ=255)

                                                         aaCGGTTGGACCTTTGATGACAGGTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGAAGCAGCTTTCTCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGGGGGACAAAAAGAACCACAcc                                                                                           <  1:991620/150‑1 (MQ=255)

                                                                ggACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCa                                                                                    <  1:1569592/150‑1 (MQ=255)

                                                                 gACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAg                                                                                   <  1:1348367/150‑1 (MQ=255)

                                                                   ccTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCTAACTACAGCGGACGAGTTACAGTTCCCTTGCTGTGCGACAAAAAGAACCACACCATCGTCAGCa                                                                                 >  2:1260641/1‑150 (MQ=255)

                                                                   ccTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCa                                                                                 >  1:1456232/1‑150 (MQ=255)

                                                                     tttGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCAGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAAc                                                                               <  2:347038/150‑1 (MQ=255)

                                                                       tgatgaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGa                                                                             >  2:1375313/1‑150 (MQ=255)

                                                                        gatgaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAATTTTTTTTTTTTATTTTTTTTTTGTCCACGGAGATCTACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGaa                                                                            <  2:926188/150‑1 (MQ=255)

                                                                        gatgaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGaa                                                                            >  1:1671118/1‑150 (MQ=255)

                                                                          tgaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATc                                                                          >  2:2038116/1‑150 (MQ=255)

                                                                          tgaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATc                                                                          >  1:1235249/1‑150 (MQ=255)

                                                                          tgaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATc                                                                          >  2:1104525/1‑150 (MQ=255)

                                                                           gaCAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATAATCTTTTTCTCCACGCAGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATca                                                                         <  2:1810171/150‑1 (MQ=255)

                                                                             cAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATcagc                                                                       >  2:1848965/1‑150 (MQ=255)

                                                                             cAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATcagc                                                                       >  1:2034434/1‑150 (MQ=255)

                                                                             cAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATcagc                                                                       <  1:1937520/150‑1 (MQ=255)

                                                                                ttttCCGGGATCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACTCCGATCCAAACTACAGCGGACGAGGTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGa                                                                    <  1:1852925/150‑1 (MQ=255)

                                                                                     cGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCCGCAACGAATCAGCAGAAatca                                                               >  2:926552/1‑150 (MQ=255)

                                                                                         aGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                           >  1:284666/1‑150 (MQ=255)

                                                                                           cAACCGGAGACACGCTCTATCTAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCTCACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCa                                                         <  1:715637/150‑1 (MQ=255)

                                                                                             aCCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATg                                                       >  1:867122/1‑150 (MQ=255)

                                                                                             aCCCGAGACACGCTCTATCAAAATCAATATCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACCAGTAACTGTTCCCGTCCTGTGCGACAAAAAGCACCACACCAACGTCAGCAACGAAACAGCACAAATCATCCCCAtt                                                       >  2:253987/1‑149 (MQ=255)

                                                                                                 gagaCACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCa                                                   <  2:561344/150‑1 (MQ=255)

                                                                                                     cacGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATAc                                               >  2:1666082/1‑150 (MQ=255)

                                                                                                     cacGCTCTATCAAAATGAATGTCTGTATCAGCTTTATCTCCACCCCGAACCACACTACAGCGGAGGAGTGTCTGTTCCCGTGCTGTGGGACAAACAGAACCAGACCATCGTCAGCAACGAATCAGCAGAAATCAACCGGAGGTTCAATAc                                               >  2:111725/1‑150 (MQ=255)

                                                                                                       cGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACcg                                             >  1:513126/1‑150 (MQ=255)

                                                                                                       cGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACcg                                             >  2:1375459/1‑150 (MQ=255)

                                                                                                       agCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACcg                                             <  1:1495227/149‑1 (MQ=255)

                                                                                                        gCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACcgc                                            >  1:509989/1‑150 (MQ=255)

                                                                                                            tATCAAAATGTTTTTTTTTATAAGTTTTTTGTCCACGGAGATCTACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGttt                                        <  2:1763239/150‑1 (MQ=255)

                                                                                                            aaTCAAAATGAATTTCTGTATCAGCTTAATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGtat                                        >  1:809906/2‑148 (MQ=255)

                                                                                                             aTCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTg                                       >  1:2061936/1‑150 (MQ=255)

                                                                                                              tCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGa                                      >  1:2094549/1‑150 (MQ=255)

                                                                                                               cAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCAGCCGCAGGTTCAATACcgcggtggtt                                     >  1:2008540/1‑144 (MQ=255)

                                                                                                                         ttCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGgcg                           <  1:1734291/150‑1 (MQ=255)

                                                                                                                           cTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGcgcg                         <  1:769820/150‑1 (MQ=255)

                                                                                                                             gTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGaa                       >  2:1531922/1‑150 (MQ=255)

                                                                                                                             gTATCAGCTTTATCACCACGCCGATCCACACTACAGCGGACGAGGTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGaa                       <  1:1337417/150‑1 (MQ=255)

                                                                                                                               aTCAGCTTTATCTCCACGCCGCTCCACTCTACAGCGGCCGAGTTAATGTTCCCGTGCTGTGTGACAAAAAGCACCACACCATCGTCAGCAACGAATCAGCAGAAATCAACCTCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAg                     >  1:2091648/1‑150 (MQ=255)

                                                                                                                               aTCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAACTCGGCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGaaat                     >  1:1282707/1‑149 (MQ=255)

                                                                                                                               aTCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGaacg                     >  2:1069802/1‑148 (MQ=255)

                                                                                                                                 cAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCg                   >  1:1919571/1‑150 (MQ=255)

                                                                                                                                   tcTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCggg                 <  1:416911/149‑1 (MQ=255)

                                                                                                                                    cTTTATCTCCACGCCGAACCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGt                <  2:1941760/150‑1 (MQ=255)

                                                                                                                                       tttctcCACGGCGATCTACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGGCAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGAt             <  2:1459058/148‑1 (MQ=255)

                                                                                                                                       tagctcCACGCCGATCCACACTCCAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGAt             <  1:960266/147‑1 (MQ=255)

                                                                                                                                         tctcCACGCCGATCCACACTACAGCGGACGAGTTACTGATCCCGTGCTGTGGGACAAAAAGAACCACACCATCGACAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTAGATGCGCTCGGCGCGACAGCGGGTGATta           >  1:373063/1‑150 (MQ=255)

                                                                                                                                         tatcCACGCCGATCCACACTACATCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGAAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGAGAACGCGGGTGATta           >  2:1766235/3‑150 (MQ=255)

                                                                                                                                           tcCACTCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtact         <  2:1052897/150‑1 (MQ=255)

                                                                                                                                            ccACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGATTCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtacta        >  2:68536/1‑150 (MQ=255)

                                                                                                                                            ccACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtacta        <  1:1762685/150‑1 (MQ=255)

                                                                                                                                            ccACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtacta        <  1:1597824/150‑1 (MQ=255)

                                                                                                                                             cACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGGGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtactac       <  1:1167772/150‑1 (MQ=255)

                                                                                                                                              aCGCCGATACACACTACAGCTGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTAcc      <  2:523561/150‑1 (MQ=255)

                                                                                                                                                  cGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTcc  <  2:642826/150‑1 (MQ=255)

CTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCAAAATGAATTTCTGTATCAGCTTTATCTCCACGCCGATCCACACTACAGCGGACGAGTTACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCC  >  NC_000913/3251246‑3251541

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