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.4%	4070.6 / NA	1068	M85I (ATG→ATT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (1/2); major base A (504/554); minor base T (3/3); total (509/559)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1068 total aligned reads displayed.

GGTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAACATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGATG  >  AM406671/1853326‑1853617

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

  tAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATAc                                                                                                                                               <  2:4201104/149‑1 (MQ=255)

   aGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATAcc                                                                                                                                              <  1:3106017/149‑1 (MQ=255)

   aGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATAcc                                                                                                                                              >  2:2287047/1‑149 (MQ=255)

   aGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTAAAAATAcc                                                                                                                                              >  1:5308326/1‑149 (MQ=255)

      cccGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAGACCACCATTACCAATCGCCCTGTCTGCTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAACATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTaa                                                                                                                                           >  2:3620753/1‑149 (MQ=255)

         gTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAAtt                                                                                                                                        >  2:5423865/1‑149 (MQ=255)

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

            gACGCTAAGGAATCCATAAAACAACTTGTCAAGCGCCCTAAACCACCATTACCATTCGCCATGTCTGGTTCAACTGTAGCTACTTCATCAAGTTCAATGCCAAAATCTGCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAg                                                                                                                                     <  1:2842449/149‑1 (MQ=255)

             aCGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGt                                                                                                                                    >  2:6079231/1‑149 (MQ=255)

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

               gCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTaa                                                                                                                                  <  1:2605572/149‑1 (MQ=255)

               gCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTaa                                                                                                                                  <  1:2960386/149‑1 (MQ=255)

                cTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTGAATGCCAAAATCTTGCAGTCCTTCGCGAAGGGTTTTTAAAATACCTAAATTAAGTaaa                                                                                                                                 >  1:6637335/1‑149 (MQ=255)

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

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

                   aGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTg                                                                                                                              >  2:7336623/1‑149 (MQ=255)

                   aGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTg                                                                                                                              >  2:2935044/1‑149 (MQ=255)

                      aaTACATAAATCAACTTGCCAAGCGCCCCAAGCCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCGTCCAGACCTTCGCGAACGTTTTTTAAAATACCTAAATTCAGTAAATTGGAt                                                                                                                           <  2:140710/149‑1 (MQ=255)

                       aTCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCACTCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGAtt                                                                                                                          >  1:5894110/1‑149 (MQ=255)

                       aTCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGAtt                                                                                                                          >  1:5163176/1‑149 (MQ=255)

                         ccATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTa                                                                                                                        >  1:1984069/1‑149 (MQ=255)

                          cATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTaa                                                                                                                       <  2:7544809/149‑1 (MQ=255)

                           aTAAAACAACTTGCCAAGCGCCCTAAACTACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTCTTTAAAATACCTAAATTAAGTCAATTGGATTTaat                                                                                                                      >  2:481734/1‑148 (MQ=255)

                            tAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGc                                                                                                                     >  1:1857155/1‑149 (MQ=255)

                              aaaCAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCAt                                                                                                                   <  2:7039651/149‑1 (MQ=255)

                               aacaacTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCAtt                                                                                                                  <  2:1810752/149‑1 (MQ=255)

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

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

                                      tGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCg                                                                                                           <  1:1759966/149‑1 (MQ=255)

                                        ccAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGt                                                                                                         <  2:6454330/149‑1 (MQ=255)

                                         cAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTa                                                                                                        <  2:8577816/149‑1 (MQ=255)

                                           aGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaaa                                                                                                      >  2:6808620/1‑149 (MQ=255)

                                                ccTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTc                                                                                                 >  1:6886660/1‑149 (MQ=255)

                                                  tAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCaa                                                                                               >  1:1146831/1‑149 (MQ=255)

                                                  tAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCaa                                                                                               <  2:2381495/149‑1 (MQ=255)

                                                   aaaCCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAAt                                                                                              <  1:9050131/149‑1 (MQ=255)

                                                    aaCCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAAtt                                                                                             <  2:4357748/149‑1 (MQ=255)

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

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

                                                       caccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaa                                                                                          <  2:7910598/149‑1 (MQ=255)

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

                                                         ccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATAACTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaa                                                                                        >  2:4517975/1‑149 (MQ=255)

                                                          caTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGGTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGGTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaaa                                                                                       >  1:1447541/1‑149 (MQ=255)

                                                          caTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCACCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaaa                                                                                       >  2:6194749/1‑149 (MQ=255)

                                                           aTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAAtt                                                                                             <  2:1511300/142‑1 (MQ=255)

                                                            ttACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaata                                                                                     <  1:1005384/149‑1 (MQ=255)

                                                               ccAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACAGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAataata                                                                                  >  1:5810005/1‑149 (MQ=255)

                                                                cAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAAtaataa                                                                                 >  2:8222467/1‑149 (MQ=255)

                                                                cAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAAtaataa                                                                                 >  1:7017799/1‑149 (MQ=255)

                                                                 aaTCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAg                                                                                <  1:3863897/149‑1 (MQ=255)

                                                                   tCGCGATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCt                                                                              <  2:5605262/149‑1 (MQ=255)

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

                                                                       cATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtt                                                                          >  1:929020/1‑149 (MQ=255)

                                                                        aTGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGttt                                                                         <  2:6569907/149‑1 (MQ=255)

                                                                        aTGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAAATAAGTAAATTGAATTTAAGCATTTTCCCCGGTAAAAAATAAATAAAAAAATAATAAGCTTGttt                                                                         <  2:8938068/149‑1 (MQ=255)

                                                                          gTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGttttt                                                                       >  2:5890261/1‑149 (MQ=255)

                                                                           tCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtttttt                                                                      >  1:395477/1‑149 (MQ=255)

                                                                             gggTGCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTcc                                                                    <  2:6339191/148‑1 (MQ=255)

                                                                               gTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCa                                                                  <  2:3285442/149‑1 (MQ=255)

                                                                               gTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCa                                                                  <  1:7146927/149‑1 (MQ=255)

                                                                                ttCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAg                                                                 <  1:2756173/149‑1 (MQ=255)

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

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

                                                                                       ttAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATc                                                          <  1:8908036/149‑1 (MQ=255)

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

                                                                                         aGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCt                                                        <  2:7195192/149‑1 (MQ=255)

                                                                                         aGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCt                                                        <  2:777112/149‑1 (MQ=255)

                                                                                          gCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCtt                                                       >  1:628051/1‑149 (MQ=255)

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

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

                                                                                          gCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCtt                                                       >  2:4659159/1‑149 (MQ=255)

                                                                                            tACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTtat                                                     >  1:6421600/1‑149 (MQ=255)

                                                                                             aCTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATAGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTtata                                                    <  2:9200849/149‑1 (MQ=255)

                                                                                              cTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTatat                                                   >  2:3488571/1‑149 (MQ=255)

                                                                                               ttCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATAtt                                                  <  1:8044863/149‑1 (MQ=255)

                                                                                                 catcaAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGa                                                <  1:7778592/149‑1 (MQ=255)

                                                                                                 catcaAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGa                                                >  1:8521420/1‑149 (MQ=255)

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

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

                                                                                                         tCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGa                                        >  1:7393967/1‑149 (MQ=255)

                                                                                                         tCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGa                                        <  1:6087530/149‑1 (MQ=255)

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

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

                                                                                                            aTGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAAc                                     <  1:3712759/149‑1 (MQ=255)

                                                                                                              gCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAAccc                                   >  2:1167278/1‑149 (MQ=255)

                                                                                                                     tCTTCCAGACCTTCGCGCACGGTTTTTAAAATACGTACATTAAGTAAATTGGATTTAGGCTTTTTCCCCGGTAATAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCCTTTTCATCCTTATATTGATTATCTGAAACCTACATAtt                            >  2:2761091/1‑149 (MQ=255)

                                                                                                                       ttCCAGGCCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCATGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAg                          <  2:5234537/149‑1 (MQ=255)

                                                                                                                         ccAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTCAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCa                        >  1:7590868/1‑149 (MQ=255)

                                                                                                                           aGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAta                      >  1:7898824/1‑149 (MQ=255)

                                                                                                                            gACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAACAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATGTTAGCAtag                     <  1:4756201/149‑1 (MQ=255)

                                                                                                                             aCCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAtagt                    >  2:2134682/1‑149 (MQ=255)

                                                                                                                             aCCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAtagt                    <  1:8252901/149‑1 (MQ=255)

                                                                                                                                ttCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATagtagt                 <  2:1340612/149‑1 (MQ=255)

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

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

                                                                                                                                        cGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcaca         >  2:8072381/1‑149 (MQ=255)

                                                                                                                                        cGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcaca         >  1:3491620/1‑149 (MQ=255)

                                                                                                                                        cGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcaca         >  1:6764687/1‑149 (MQ=255)

                                                                                                                                             tttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTTTTTCACAACAGa    >  2:3144956/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` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`