BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM406671	1,853,472	C→A	M85I (ATG→ATT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM406671	1,853,472	0	C	A	99.6%	3838.3 / NA	990	M85I (ATG→ATT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (0/1); major base A (472/513); minor base T (3/1); total (475/515)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 990 total aligned reads displayed.

GGTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAACATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGATG  >  AM406671/1853326‑1853617

ggTCGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAGATCTTCCAGACCGTCGCGAACGGTTTTTAAAAt                                                                                                                                                 <  2:8980045/149‑1 (MQ=255)

ggTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAAt                                                                                                                                                 <  2:5303867/149‑1 (MQ=255)

 gTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATa                                                                                                                                                >  2:2796516/1‑149 (MQ=255)

 gTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATa                                                                                                                                                <  1:615548/149‑1 (MQ=255)

  tAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATAc                                                                                                                                               >  1:4788751/1‑149 (MQ=255)

    ggCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCt                                                                                                                                             <  1:5881371/149‑1 (MQ=255)

     gCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACCTCATCAAGTTCAATGCCATAATCTTCCAGTCCTTCGCGAACGGTTTTTAAAATACCTa                                                                                                                                            >  1:6951051/1‑149 (MQ=255)

           aGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTaa                                                                                                                                      >  2:809536/1‑149 (MQ=255)

           aGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTaa                                                                                                                                      >  1:1376166/1‑149 (MQ=255)

           aGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTaa                                                                                                                                      <  2:3062813/149‑1 (MQ=255)

            gACGCTAAGGAATCCAGAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAg                                                                                                                                     <  1:6450618/149‑1 (MQ=255)

             aCGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGt                                                                                                                                    <  2:3478516/149‑1 (MQ=255)

                cTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTaaa                                                                                                                                 <  1:926234/149‑1 (MQ=255)

                cTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTaaa                                                                                                                                 <  2:1142099/149‑1 (MQ=255)

                  aaGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAAtt                                                                                                                               >  2:7825288/1‑149 (MQ=255)

                  aaGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAAtt                                                                                                                               <  2:4417775/149‑1 (MQ=255)

                   aGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTg                                                                                                                              <  1:5714952/149‑1 (MQ=255)

                     gAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGa                                                                                                                            <  2:1593792/149‑1 (MQ=255)

                       aTCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGAtt                                                                                                                          <  1:3416281/149‑1 (MQ=255)

                         ccATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTa                                                                                                                        <  2:9331619/149‑1 (MQ=255)

                           aTAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCGAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAg                                                                                                                      >  1:382618/1‑149 (MQ=255)

                           aTAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAg                                                                                                                      <  1:4661387/149‑1 (MQ=255)

                            tAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGc                                                                                                                     <  1:1128366/149‑1 (MQ=255)

                              aaaCAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCAt                                                                                                                   >  1:154588/1‑149 (MQ=255)

                                acaacTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCAttt                                                                                                                 >  1:2891715/1‑149 (MQ=255)

                                  aacTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTc                                                                                                               >  1:5152241/1‑149 (MQ=255)

                                    cTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTccc                                                                                                             >  2:7592135/1‑149 (MQ=255)

                                    cTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTccc                                                                                                             <  1:9390165/149‑1 (MQ=255)

                                      tGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCg                                                                                                           >  1:3995069/1‑149 (MQ=255)

                                        ccAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGt                                                                                                         >  1:8840313/1‑149 (MQ=255)

                                          aaGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaa                                                                                                       >  1:3052327/1‑149 (MQ=255)

                                          aaGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaa                                                                                                       >  1:777940/1‑149 (MQ=255)

                                           aGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaaa                                                                                                      <  2:1969200/149‑1 (MQ=255)

                                           aGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaaa                                                                                                      <  2:2932481/149‑1 (MQ=255)

                                               cccTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACt                                                                                                  >  2:1794759/1‑149 (MQ=255)

                                               cccTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACt                                                                                                  <  1:9128151/149‑1 (MQ=255)

                                                 cTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCa                                                                                                >  2:4108863/1‑149 (MQ=255)

                                                  tAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAGATTGGATTTAAGCATTTTCCCCGGTAAAAACTCaa                                                                                               <  1:4145865/149‑1 (MQ=255)

                                                     accaccATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTATTATACCTATATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTg                                                                                            >  2:4235482/1‑149 (MQ=255)

                                                     accaccATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTg                                                                                            <  2:2114759/149‑1 (MQ=255)

                                                        accaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaa                                                                                         <  2:7978570/149‑1 (MQ=255)

                                                         ccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaa                                                                                        >  2:4646471/1‑149 (MQ=255)

                                                         ccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaa                                                                                        >  1:759342/1‑149 (MQ=255)

                                                         ccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaa                                                                                        <  1:4974880/149‑1 (MQ=255)

                                                           aTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaat                                                                                      >  1:255270/1‑149 (MQ=255)

                                                              aCCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaataat                                                                                   <  2:9135375/149‑1 (MQ=255)

                                                              aCCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaataat                                                                                   <  1:1259744/149‑1 (MQ=255)

                                                                 aaTCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAg                                                                                >  1:794510/1‑149 (MQ=255)

                                                                  aTCGCCATGTCTGGTTCAACTTTAGCTANTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGGAAAAACTCAATTGAAAAATAATAAGc                                                                               >  1:5374808/1‑149 (MQ=255)

                                                                   tCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCt                                                                              >  2:508916/1‑149 (MQ=255)

                                                                     gCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTAAATAAAAAAATAATAAGCTTg                                                                            <  2:8171538/149‑1 (MQ=255)

                                                                      ccATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCTGGTAAAAACTCAATTGAAAAATAATAAGCTTGt                                                                           <  2:5570604/149‑1 (MQ=255)

                                                                      ccATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGt                                                                           >  2:3216311/1‑149 (MQ=255)

                                                                      ccATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGt                                                                           >  1:3607132/1‑149 (MQ=255)

                                                                       cATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtt                                                                          <  1:3278351/149‑1 (MQ=255)

                                                                        aTGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGttt                                                                         >  1:1859610/1‑149 (MQ=255)

                                                                         tGTCTGGTTCAACTTTAGCAACTTCATCAAGTTCAATGCCAAAATATTACAGAACTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCACAGTAAAAACTCAATTGAAAAATAATAAGCTTGtttt                                                                        <  2:6590327/149‑1 (MQ=255)

                                                                          gTCTGGATCAACTTTAGCTACTTCATCCAGTTCAATGCCAAAATCTTCCAGACGTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGAATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTGGttttt                                                                       <  2:6365454/149‑1 (MQ=255)

                                                                               gTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCa                                                                  >  1:2307148/1‑149 (MQ=255)

                                                                               gTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCa                                                                  >  2:8515153/1‑149 (MQ=255)

                                                                                 tCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGt                                                                <  2:687970/149‑1 (MQ=255)

                                                                                 tCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGt                                                                <  1:8310824/149‑1 (MQ=255)

                                                                                  cAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGtt                                                               >  2:7387119/1‑149 (MQ=255)

                                                                                  cAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGtt                                                               <  1:5507150/149‑1 (MQ=255)

                                                                                    aCTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTc                                                             >  2:489068/1‑149 (MQ=255)

                                                                                    aCTTTAGCTACTTCATAAAGTTCAATTCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTAAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTc                                                             <  2:3933391/149‑1 (MQ=255)

                                                                                       ttAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATc                                                          >  1:625716/1‑149 (MQ=255)

                                                                                        tAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATcc                                                         >  1:728209/1‑149 (MQ=255)

                                                                                        tAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTCTTTCCCAGTTTCATcc                                                         >  2:1020216/1‑149 (MQ=255)

                                                                                          gCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCtt                                                       <  2:1424009/149‑1 (MQ=255)

                                                                                           cTACTTCCTCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTta                                                      <  2:1313268/149‑1 (MQ=255)

                                                                                           cTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTta                                                      >  1:6713127/1‑149 (MQ=255)

                                                                                             aCTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTtata                                                    >  2:6181137/1‑149 (MQ=255)

                                                                                                tcatcaAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTg                                                 <  1:6032728/149‑1 (MQ=255)

                                                                                                    caAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATTCTTATATTGATTa                                             >  2:5808793/1‑149 (MQ=255)

                                                                                                    caAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTa                                             <  1:2699302/149‑1 (MQ=255)

                                                                                                       gTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCa                                          >  1:6586368/1‑149 (MQ=255)

                                                                                                       gTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCa                                          >  1:7389837/1‑149 (MQ=255)

                                                                                                       gTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCa                                          >  2:5995463/1‑149 (MQ=255)

                                                                                                          cAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGaa                                       <  1:1007281/149‑1 (MQ=255)

                                                                                                           aaTGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGaaa                                      >  1:213376/1‑149 (MQ=255)

                                                                                                            aTGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAAc                                     >  1:9576162/1‑149 (MQ=255)

                                                                                                            aTGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAAc                                     >  2:8859603/1‑149 (MQ=255)

                                                                                                             tGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAAcc                                    >  1:9010137/1‑149 (MQ=255)

                                                                                                                 aaaaTCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGa                                <  2:7146221/149‑1 (MQ=255)

                                                                                                                   aaTCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGata                              <  1:809536/149‑1 (MQ=255)

                                                                                                                       ttCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAg                          >  1:8497298/1‑149 (MQ=255)

                                                                                                                          cAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAt                       <  1:1981702/149‑1 (MQ=255)

                                                                                                                            gACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAtag                     <  2:9499998/149‑1 (MQ=255)

                                                                                                                            gACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTCTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAtag                     >  2:3669713/1‑149 (MQ=255)

                                                                                                                                 tCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCTTTGTTTTTTCCTTATATTGATTATCATAAACCCAGATATTCGCATagtagtt                >  2:8334310/1‑148 (MQ=255)

                                                                                                                                   gcgAACGGTTTTTAAAATACCAAAATTAAGTAAATTAGATTTAAGCAGTTTCCCCGGTAAAAACTCAATTAAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGtt              <  2:8675294/149‑1 (MQ=255)

                                                                                                                                    cgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGttt             <  2:3316497/149‑1 (MQ=255)

                                                                                                                                      aaCGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTca           >  2:267686/1‑149 (MQ=255)

                                                                                                                                       aCGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcac          <  1:7098552/149‑1 (MQ=255)

                                                                                                                                        cGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcaca         <  1:8670872/149‑1 (MQ=255)

                                                                                                                                          gTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCacaac       >  2:237606/1‑149 (MQ=255)

                                                                                                                                           tttttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCATAAACCCAGATATTAGCATAGTAGTGTTTCacaaca      >  2:6520263/1‑149 (MQ=255)

                                                                                                                                            ttttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAg     >  1:1639525/1‑149 (MQ=255)

                                                                                                                                             tttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGa    >  1:4383684/1‑149 (MQ=255)

GGTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAACATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGATG  >  AM406671/1853326‑1853617

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