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.2%	4408.3 / NA	1153	M85I (ATG→ATT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (563/579); minor base T (7/2); total (570/583)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1153 total aligned reads displayed.

GGTAGGCCCGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAACATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGATG  >  AM406671/1853326‑1853617

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

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

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

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

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

        cGTAGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAAt                                                                                                                                         >  2:23490/1‑149 (MQ=255)

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

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

           aGACGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATCCCAAAATCTTCCAGACCCTCGCGAACGGTTTTTAAAATACCTAAATTaa                                                                                                                                      >  2:3584045/1‑149 (MQ=255)

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

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

              cGCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTa                                                                                                                                   <  2:9721041/149‑1 (MQ=255)

               gCTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTaa                                                                                                                                  >  1:4425481/1‑149 (MQ=255)

                cTAAGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCAGCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTaaa                                                                                                                                 <  2:10274983/149‑1 (MQ=255)

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

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

                   aGGAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGATCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTg                                                                                                                              <  2:7707788/149‑1 (MQ=255)

                    ggAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTgg                                                                                                                             <  1:7854150/149‑1 (MQ=255)

                    ggAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTgg                                                                                                                             <  1:9458653/149‑1 (MQ=255)

                    ggAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTgg                                                                                                                             >  1:9800868/1‑149 (MQ=255)

                    ggAATCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTgg                                                                                                                             <  1:4175209/149‑1 (MQ=255)

                      aaTCCATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGAt                                                                                                                           <  1:6218125/149‑1 (MQ=255)

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

                          cATAAAACAACTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTaa                                                                                                                       >  2:2173288/1‑149 (MQ=255)

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

                               aacaacTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCAtt                                                                                                                  >  2:8509374/1‑149 (MQ=255)

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

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

                                 caacTTGCCAAGCGGCCTAAACCTCCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATGTTCCAGACCTTGGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCAtttt                                                                                                                <  1:1205820/149‑1 (MQ=255)

                                   acatGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTcc                                                                                                              <  1:8639511/146‑1 (MQ=255)

                                   acTTGCCAAGCGCCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTcc                                                                                                              <  1:9163847/149‑1 (MQ=255)

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

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

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

                                            gcgcCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaaaa                                                                                                     >  1:2151877/1‑149 (MQ=255)

                                            gcgcCCTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTaaaa                                                                                                     >  1:2445951/1‑149 (MQ=255)

                                               cccTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAAGCTTCCAGACCTTCGCGAACGGTTTCTAAAATACCTAAAGTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACt                                                                                                  <  2:8214746/149‑1 (MQ=255)

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

                                                 cTAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATGAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCa                                                                                                <  2:1844135/149‑1 (MQ=255)

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

                                                  tAAACCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAAGCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAAGTGGATTGAAGCATTTTCCCCGGTAAAAACTCaa                                                                                               <  2:1523691/149‑1 (MQ=255)

                                                  tAAACCACCATTACCAATCGCCATGTCTGGATCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCaa                                                                                               >  1:1055255/1‑149 (MQ=255)

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

                                                    aaCCACCATTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAAtt                                                                                             >  1:1000865/1‑149 (MQ=255)

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

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

                                                      ccaccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATATCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTATAAACTCAATTGa                                                                                           >  2:476134/1‑149 (MQ=255)

                                                      ccaccaTTACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGa                                                                                           >  1:7482913/1‑149 (MQ=255)

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

                                                             tACCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTaaaataaaaaataa                                                                                    <  2:8418288/149‑1 (MQ=255)

                                                              aCCAATCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAaataat                                                                                   >  2:7445853/1‑149 (MQ=255)

                                                                 aaTCGCCATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAACACTCAATTGAATAATAATAAg                                                                                >  1:7642958/1‑149 (MQ=255)

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

                                                                      ccATGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTTAATTTCTTAAATTAAGTAAATTGGATTAAAGCATATTCCCCGTTAAAAACTTAATTAAAAAATAATAATCTGGt                                                                           >  2:3968815/1‑149 (MQ=255)

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

                                                                         tGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtttt                                                                        >  1:1600399/1‑149 (MQ=255)

                                                                         tGTCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtttt                                                                        >  1:5097678/1‑149 (MQ=255)

                                                                           tCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtttttt                                                                      <  1:5695233/149‑1 (MQ=255)

                                                                           tCTGGTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGtttttt                                                                      <  2:1728864/149‑1 (MQ=255)

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

                                                                              ggTTCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTccc                                                                   <  1:1989020/149‑1 (MQ=255)

                                                                                 tCAACTTTAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGAGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGt                                                                <  1:5560321/149‑1 (MQ=255)

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

                                                                                      tttAGCTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCAt                                                           <  2:8764545/149‑1 (MQ=255)

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

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

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

                                                                                           cTACTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTTTTTTTTTTTAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAAGAATATGCTTGTTTTTTCCCAGTTTCATCCTta                                                      >  2:5745750/1‑149 (MQ=255)

                                                                                              cTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATcc                                                         >  2:9228160/1‑143 (MQ=255)

                                                                                              cTTCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTta                                                      <  2:8007568/146‑1 (MQ=255)

                                                                                               ttCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATAtt                                                  >  1:6468086/1‑149 (MQ=255)

                                                                                               ttCATCAAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTGGTTTTTTCCCAGTTTCATCCTtatact                                                  <  2:4417282/149‑3 (MQ=255)

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

                                                                                                  atcaAGTTGAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGAt                                               <  1:6614532/149‑1 (MQ=255)

                                                                                                   tcaAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGAtt                                              <  2:9599923/149‑1 (MQ=255)

                                                                                                    caAGTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGAttt                                             >  1:3788634/1‑148 (MQ=255)

                                                                                                       gTTCAATGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCa                                          <  2:10426004/149‑1 (MQ=255)

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

                                                                                                         tCAATGCCAAAATCTTCCAGACCTTCGCGAACTGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGTTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTGCCAGTTGCATCCTTATATTGTGTATCaaa                                        <  2:6861725/149‑3 (MQ=255)

                                                                                                           aaTGCCAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGGTTCATCCTTATATTGATTATCAGaaa                                      <  2:5135102/149‑1 (MQ=255)

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

                                                                                                               ccAAAATCTTCCAGACCTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCa                                  <  1:3598233/149‑1 (MQ=255)

                                                                                                                aaaaaGCTTCCAGACATTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTAAATTTAAGCATTTTCCCCGGTAAAAACAAAAAAGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATACTTATATTGATTATCAGAAACCCAg                                 <  2:6272449/148‑1 (MQ=255)

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

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

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

                                                                                                                              ccTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGa                                <  1:7242064/136‑1 (MQ=255)

                                                                                                                              ccTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAtagta                   <  1:10751087/149‑1 (MQ=255)

                                                                                                                              ccTTCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCAtagta                   <  2:3326527/149‑1 (MQ=255)

                                                                                                                                ttCGCGAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATagtagt                 >  2:760587/1‑149 (MQ=255)

                                                                                                                                 tCGCGAACGGTTTTTAAAATACCTCAATTAATTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTAGATTGAGTATCAGAAACCCAGACATTAGCATAGTAgtg                <  1:5987638/149‑1 (MQ=255)

                                                                                                                                   gcgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGtt              >  2:10733193/1‑149 (MQ=255)

                                                                                                                                   gcgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGtt              <  1:8103912/149‑1 (MQ=255)

                                                                                                                                   gcgAACGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGtt              <  2:5557283/149‑1 (MQ=255)

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

                                                                                                                                       aCGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAAGTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATTGTAGTGTTTcac          >  1:10147278/1‑149 (MQ=255)

                                                                                                                                        cGGTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCGTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcaca         <  2:7186555/149‑1 (MQ=255)

                                                                                                                                         ggTTTTTAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCacaa        <  2:1974643/149‑1 (MQ=255)

                                                                                                                                             tttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTcac          <  1:8425683/143‑1 (MQ=255)

                                                                                                                                              ttAAAATACCTAAATTAAGTAAATTGGATTTAAGCATTTTCCCCGGTAAAAACTCAATTGAAAAATAATAAGCTTGTTTTTTCCCAGTTTCATCCTTATATTGATTATCAGAAACCCAGATATTAGCATAGTAGTGTTTCACAACAGAt   <  1:9635075/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: `‑`