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%	3937.7 / NA	1031	M85I (ATG→ATT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (0/3); major base A (490/533); minor base T (3/1); total (493/538)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1031 total aligned reads displayed.

GGTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAACATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGATG  >  AM406671/1853326‑1853617

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

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

     gCCCGTAGACGCTAAGGAATCCATAAAANAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTa                                                                                                                                            >  1:5764273/1‑149 (MQ=255)

     gCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTa                                                                                                                                            <  2:7776799/149‑1 (MQ=255)

     gCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTa                                                                                                                                            >  1:167522/1‑149 (MQ=255)

     gCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTa                                                                                                                                            <  2:1139305/149‑1 (MQ=255)

       ccGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTaaa                                                                                                                                          <  2:1346598/149‑1 (MQ=255)

        cGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAAt                                                                                                                                         <  2:7669799/149‑1 (MQ=255)

         gTAGACCCTAATGAATCCATAAAACAACTTGCCAATCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAAtt                                                                                                                                        <  2:6847/149‑1 (MQ=255)

          tAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTa                                                                                                                                       <  1:2262039/149‑1 (MQ=255)

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

            gACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAg                                                                                                                                     <  1:4977592/149‑1 (MQ=255)

            aaCGCTAAGGAATCCATAAAGCAACTTGCCAAGCGCCCTAAACCCCCATTAACAATCGCCATGTCTGGTTCAACTTTAGCTACATCAGCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATAAAg                                                                                                                                     <  2:419609/148‑1 (MQ=255)

             aCGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAATATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGt                                                                                                                                    <  1:1495810/149‑1 (MQ=255)

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

              cGCTAAAGAATCCATAAAACAACTTGCCCAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCCATGCCAAAATCTTCCAAACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTa                                                                                                                                   <  2:1905958/149‑1 (MQ=255)

                 tAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAAt                                                                                                                                >  2:4296042/1‑149 (MQ=255)

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

                    ggAATCCATAAAACAACTTGCCAAGCGCCATAAAACAACATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTgg                                                                                                                             <  2:847157/149‑1 (MQ=255)

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

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

                        tCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGCTTTTAAAATACCTAAATTAAGTAAATTGGAttt                                                                                                                         <  2:6045472/149‑1 (MQ=255)

                             aaaaCAACTTGCCAAGCGCACTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGCTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCa                                                                                                                    >  2:2466853/1‑149 (MQ=255)

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

                                     ttGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTcccc                                                                                                            <  1:7017053/149‑1 (MQ=255)

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

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

                                         cAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAATTTGGATTTAAGCATTTTCCCCGGTa                                                                                                        >  1:6515589/1‑149 (MQ=255)

                                         cAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTa                                                                                                        >  2:1571347/1‑149 (MQ=255)

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

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

                                                ccTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTCGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTc                                                                                                 >  2:7045957/1‑149 (MQ=255)

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

                                                   aaaCCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAAt                                                                                              >  1:4329274/1‑149 (MQ=255)

                                                   aaaCCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAAt                                                                                              >  1:3585060/1‑149 (MQ=255)

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

                                                     accaccATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTg                                                                                            >  2:9796968/1‑149 (MQ=255)

                                                      ccaccaTTCCCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGa                                                                                           <  1:6020375/149‑1 (MQ=255)

                                                      ccaccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGa                                                                                           <  1:2496852/149‑1 (MQ=255)

                                                       caccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaa                                                                                          >  1:7736571/1‑149 (MQ=255)

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

                                                          caTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGaaaaa                                                                                       >  2:7433010/1‑149 (MQ=255)

                                                            ttACCAATTGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaata                                                                                     <  2:8735522/149‑1 (MQ=255)

                                                             tACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaataa                                                                                    <  1:7914432/149‑1 (MQ=255)

                                                             tACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaataa                                                                                    >  1:3134447/1‑149 (MQ=255)

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

                                                               ccAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAataata                                                                                  <  1:6848438/149‑1 (MQ=255)

                                                               ccAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAataata                                                                                  <  1:6367431/149‑1 (MQ=255)

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

                                                                   tCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCt                                                                              <  1:2060202/149‑1 (MQ=255)

                                                                     gCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTg                                                                            <  1:2984275/149‑1 (MQ=255)

                                                                     gCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTg                                                                            <  2:10413601/149‑1 (MQ=255)

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

                                                                        aTGTCTGGTTCAACTTTAGCTACTTCATCAATTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAAAACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGttt                                                                         <  2:7197031/149‑1 (MQ=255)

                                                                        aTGTCTGGTTCAACTTTAGCTACTTCATAAAGTACAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGAATATAAGCATTTTCCCCGGTAAAAAAGAAAATAAAAAATAATAAGCTTGttt                                                                         <  2:5706190/149‑1 (MQ=255)

                                                                                ttCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAg                                                                 >  1:4793382/1‑149 (MQ=255)

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

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

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

                                                                                   aaCTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGttt                                                              <  1:4182168/149‑1 (MQ=255)

                                                                                   aaCTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGttt                                                              >  2:5481065/1‑149 (MQ=255)

                                                                                     cTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTTCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCa                                                            <  1:919943/149‑1 (MQ=255)

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

                                                                                        tAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATcc                                                         <  1:3798244/149‑1 (MQ=255)

                                                                                          gCTAATTCATCAAGTTCAATGCCAAAATCTTCCACACCTTCGCGAACGGTTTTTAAATTACTTAAATTATTTAAATTGAATTTAAGCACTTTCCCCGGTAAAAACACAAAAGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCtt                                                       <  2:6166612/149‑1 (MQ=255)

                                                                                           cTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTta                                                      <  1:2732900/149‑1 (MQ=255)

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

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

                                                                                                tcatcaAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTg                                                 >  2:4684397/1‑149 (MQ=255)

                                                                                                 catcaAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCt                                                        <  1:5887992/141‑1 (MQ=255)

                                                                                                    caAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTa                                             >  1:9339220/1‑149 (MQ=255)

                                                                                                     aaGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATCGGAGTTAAGCATTTTCTCCGGTAAAAACTCAATTGAAAAATAATAAGGTTGTTTTTTCTCAGTTTCATCTTTATATTtattat                                            >  2:6628472/1‑149 (MQ=255)

                                                                                                      aGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTAAAATAAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATc                                           <  2:5234803/149‑1 (MQ=255)

                                                                                                       gTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAATTTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCa                                          >  2:5887817/1‑149 (MQ=255)

                                                                                                        ttCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAg                                         <  1:3342920/149‑1 (MQ=255)

                                                                                                         tCAATGCCAAAATCTTCCAGACCTTCGCGAACGATTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGa                                        <  2:2943156/149‑1 (MQ=255)

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

                                                                                                              gCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAAccc                                   <  1:7617306/149‑1 (MQ=255)

                                                                                                                 aaaaTCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATGAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGa                                <  2:3302613/149‑1 (MQ=255)

                                                                                                                  aaaTCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGat                               <  2:3818672/149‑1 (MQ=255)

                                                                                                                  aaaTCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGat                               <  2:9061282/149‑1 (MQ=255)

                                                                                                                  aaaTCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGat                               <  1:4023573/149‑1 (MQ=255)

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

                                                                                                                     tCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATAtt                            <  1:10008819/149‑1 (MQ=255)

                                                                                                                        tCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGc                         >  2:3608034/1‑149 (MQ=255)

                                                                                                                         ccAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAATCCAGATATTAGCa                        >  2:3432565/1‑149 (MQ=255)

                                                                                                                         ccAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCa                        <  1:5386069/149‑1 (MQ=255)

                                                                                                                         ccAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCa                        >  1:6725321/1‑149 (MQ=255)

                                                                                                                           aGACCCTCGCGAACAGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCCGAAACCCAGATACTAGCAta                      <  2:2095108/149‑1 (MQ=255)

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

                                                                                                                                  cgcgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAgtgt               >  2:4027636/1‑149 (MQ=255)

                                                                                                                                  cgcgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAgtgt               <  1:9088167/149‑1 (MQ=255)

                                                                                                                                    cgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCATGTTTTTTCCCAGTTTCATCCGTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGttt             >  2:2686570/1‑149 (MQ=255)

                                                                                                                                     gtactgTTTTTAACATCCCTAAATTAAGTAAATTGGATTTAAGCAGTTTCCCCGTTAAAAACTCAATTGCAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTCTATTGTTCATCAGAAACACAGATATTAGCATAGTAGTGTTTc            <  2:1784115/144‑1 (MQ=255)

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

                                                                                                                                      aaCGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTca           <  2:9601873/149‑1 (MQ=255)

                                                                                                                                       aCGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcac          >  1:8755589/1‑149 (MQ=255)

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

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

                                                                                                                                             tttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGa    <  1:8280853/149‑1 (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: `‑`