BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM406671	1,210,283	(A)_5→6	pseudogene (515/1046 nt)	llmg_pseudo_39 ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM406671	1,210,278	1	.	A	98.0%	3090.6 / 36.9	784	pseudogene (520/1046 nt)	llmg_pseudo_39	pseudogene
Reads supporting (aligned to +/- strand): ref base . (12/4); new base A (398/370); total (410/374)
Fisher's exact test for biased strand distribution p-value = 7.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 843 total aligned reads displayed.

TTCAGTCAAAATTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑A‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCATCAA  >  AM406671/1210130‑1210421

 tCAGTCAAAATTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATCGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAgag‑a                                                                                                                                                  >  2:6750492/1‑149 (MQ=255)

   aGTCAAAATTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑a‑aa                                                                                                                                               <  1:6916312/149‑1 (MQ=255)

     tCAAAATTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑a‑aaa‑a                                                                                                                                            <  2:6919774/149‑1 (MQ=255)

        aaaTTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑a‑aaa‑aata                                                                                                                                         >  1:4656434/1‑146 (MQ=255)

         aaTTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑a‑aaa‑aatat                                                                                                                                        >  2:8068377/1‑145 (MQ=255)

         aaTTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑a‑aaa‑aatat                                                                                                                                        >  1:2699417/1‑145 (MQ=255)

                   aaaCTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAg                                                                                                                               >  1:5526403/1‑149 (MQ=255)

                         tCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGt                                                                                                                         >  1:4241406/1‑149 (MQ=255)

                            gTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTaa                                                                                                                      >  2:6438500/1‑149 (MQ=255)

                              tctctcCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAAt                                                                                                                    >  2:1067112/1‑149 (MQ=255)

                               ctctcCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATc                                                                                                                   >  1:6569118/1‑149 (MQ=255)

                                tctcCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCa                                                                                                                  >  1:1245657/1‑149 (MQ=255)

                                     tGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATc                                                                                                                   <  1:826928/143‑1 (MQ=255)

                                     tGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATc                                                                                                                   >  2:826955/1‑143 (MQ=255)

                                      gAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCa                                                                                                            <  1:4856103/149‑1 (MQ=255)

                                      gAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCa                                                                                                            >  1:7235949/1‑149 (MQ=255)

                                             ccATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGaa                                                                                                     >  2:4582795/1‑149 (MQ=255)

                                              cATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAAt                                                                                                    <  2:4239254/149‑1 (MQ=255)

                                               aTGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATa                                                                                                   <  2:6879187/149‑1 (MQ=255)

                                                tGTATTGGACGAACCACCAATGAAAGCTNATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATaa                                                                                                  >  1:5138012/1‑149 (MQ=255)

                                                 gTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAg                                                                                                 >  1:2250638/1‑149 (MQ=255)

                                                  tATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAgg                                                                                                >  1:5848274/1‑149 (MQ=255)

                                                     tGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGa                                                                                             >  2:2131795/1‑149 (MQ=255)

                                                     tGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGa                                                                                             <  2:7146866/149‑1 (MQ=255)

                                                        aCGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAAtt                                                                                          <  1:4199167/149‑1 (MQ=255)

                                                        aCGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAAtt                                                                                          >  1:7719650/1‑149 (MQ=255)

                                                         cGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAAttt                                                                                         <  1:83157/149‑1 (MQ=255)

                                                          gAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTc                                                                                        <  2:3411293/149‑1 (MQ=255)

                                                           aaCCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCt                                                                                       <  2:747946/149‑1 (MQ=255)

                                                            accaccAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCtt                                                                                      <  1:1388654/149‑1 (MQ=255)

                                                               accaATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTgg                                                                                   >  2:8428253/1‑149 (MQ=255)

                                                                ccaATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGc                                                                                  >  2:5430599/1‑149 (MQ=255)

                                                                 caATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCt                                                                                 >  2:5214279/1‑149 (MQ=255)

                                                                 caATGAAAGCTAATTTCTGATAAGGCTTTTGAATCATGTGAGAAACAAGAAGCTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGAGAGCTGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCt                                                                                 >  2:5012085/1‑149 (MQ=255)

                                                                   aTGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTaa                                                                               <  1:709243/149‑1 (MQ=255)

                                                                           tAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTaa                                                                       <  2:1384324/149‑1 (MQ=255)

                                                                            aaTTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAag                                                                      <  1:3130758/149‑1 (MQ=255)

                                                                             aTTTCTGATAGCCTTTTTGTTTCATGTTAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAaga                                                                     >  2:7664633/1‑149 (MQ=255)

                                                                             aTTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAaga                                                                     >  1:2943020/1‑149 (MQ=255)

                                                                              tttCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTACATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAagag                                                                    <  1:113662/149‑1 (MQ=255)

                                                                               ttCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGagc                                                                   >  2:6272353/1‑149 (MQ=255)

                                                                                 cTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGagcag                                                                 >  2:3254721/1‑149 (MQ=255)

                                                                                 cTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGagcag                                                                 <  1:5457926/149‑1 (MQ=255)

                                                                                  tGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGagcagc                                                                <  2:405827/149‑1 (MQ=255)

                                                                                  tGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGagcagc                                                                >  1:3257160/1‑149 (MQ=255)

                                                                                   gATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCt                                                               <  1:2006734/149‑1 (MQ=255)

                                                                                    aTAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCtt                                                              >  1:6049322/1‑149 (MQ=255)

                                                                                     tAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTg                                                             <  1:7849670/149‑1 (MQ=255)

                                                                                      aGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTgg                                                            >  1:4420959/1‑149 (MQ=255)

                                                                                      aGCCTTTTTGAATCATGTGAGAAACAAGTAGTGGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATGTCTTTGGCTAAATTGCTAAGAGCAGCTTgg                                                            <  2:908154/149‑1 (MQ=255)

                                                                                        ccTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTa                                                          >  1:6208074/1‑149 (MQ=255)

                                                                                        ccTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTa                                                          <  1:8951699/149‑1 (MQ=255)

                                                                                          tttttGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGt                                                        <  1:3806007/149‑1 (MQ=255)

                                                                                           ttttGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTc                                                       <  1:9226558/149‑1 (MQ=255)

                                                                                           ttttGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTc                                                       <  1:2376468/149‑1 (MQ=255)

                                                                                            tttGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCa                                                      >  2:1804624/1‑149 (MQ=255)

                                                                                             ttGAATCATGTGAGAAACAAGTAGTTGANTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGGTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAg                                                     >  1:5059181/1‑149 (MQ=255)

                                                                                             ttGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAGATTGCTAAGAGCAGCTTGGTAGTCAg                                                     >  1:4116860/1‑149 (MQ=255)

                                                                                              tGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAgg                                                    <  1:7989860/149‑1 (MQ=255)

                                                                                              tGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCCGACAGATGAGTTAAATCACGTCCAAGTAGACTAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAgg                                                    <  2:271748/149‑1 (MQ=255)

                                                                                                 aTCATGTGAGAAACAAGTAGTTGCGTGGGCGCGCGTTGGGCTCTTTTTAGAG‑AGAAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCACGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGAAGCTTGGTAGTAAGCAtt                                                 >  2:7014002/1‑149 (MQ=255)

                                                                                                 aTCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGAtt                                                 >  2:3570691/1‑149 (MQ=255)

                                                                                                   cATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTaa                                               <  1:3955173/149‑1 (MQ=255)

                                                                                                    aTGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAat                                              <  1:1732479/149‑1 (MQ=255)

                                                                                                    aTGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAat                                              >  2:5926141/1‑149 (MQ=255)

                                                                                                      gtgAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAatat                                            <  2:3075553/149‑1 (MQ=255)

                                                                                                           aaaCAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑A‑AAACATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCCGtt                                       >  2:243989/1‑149 (MQ=255)

                                                                                                            aaCAAGTAATTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGttt                                      <  1:5292491/149‑1 (MQ=255)

                                                                                                              cAAGTAGTAGAGTGGGTGTGCGTTGGTCACTTTTTAGAGAA‑AAA‑ATATTGGCTGACAGATGAGTTAAATCACGTCCACGGAGAGTAAGGGGAATTTCTTTGGCTAAATTGCTAAGAACAGCTTGGTAGTCAGGATTAATATAAGGTTCa                                    >  2:1985147/1‑149 (MQ=255)

                                                                                                                     ttGAGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGt                                                           >  2:2313178/1‑119 (MQ=255)

                                                                                                                        aGTGGGAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtc                          <  2:1266975/149‑1 (MQ=255)

                                                                                                                          tGGGAGTGCGTTGGTCTCTTTTTAGAG‑ACAAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtcct                        <  1:6348715/149‑1 (MQ=255)

                                                                                                                          tGGGAGTGCGTTGGTCTCTGTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtcct                        <  1:4512049/149‑1 (MQ=255)

                                                                                                                          ggggAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtcct                        <  2:5139634/148‑1 (MQ=255)

                                                                                                                           gggAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtcctc                       >  2:8817982/1‑149 (MQ=255)

                                                                                                                           gggAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtcctc                       >  1:6809389/1‑149 (MQ=255)

                                                                                                                           gggAGTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAAtcctc                       <  2:5601470/149‑1 (MQ=255)

                                                                                                                               gTGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTCTGGCTAAATTGCTAAGACCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTaa                   <  2:3792931/149‑1 (MQ=255)

                                                                                                                                tGCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCATaaa                  <  1:3570691/149‑1 (MQ=255)

                                                                                                                                 gCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTaaaa                 >  1:7569684/1‑149 (MQ=255)

                                                                                                                                 gCGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGAAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTaaaa                 <  1:1022564/149‑1 (MQ=255)

                                                                                                                                  cGTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAAt                <  2:4818521/149‑1 (MQ=255)

                                                                                                                                   gTTGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAAtt               <  2:1593071/149‑1 (MQ=255)

                                                                                                                                    ttGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTa              <  2:6096840/149‑1 (MQ=255)

                                                                                                                                    ttGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATCTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTa              <  2:2948336/149‑1 (MQ=255)

                                                                                                                                     tGGTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTaa             >  1:815251/1‑149 (MQ=255)

                                                                                                                                       gTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTaaaa           <  1:2834375/149‑1 (MQ=255)

                                                                                                                                       gTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTaaaa           <  1:1909374/149‑1 (MQ=255)

                                                                                                                                       gTCTCTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTaaaa           <  1:5740890/149‑1 (MQ=255)

                                                                                                                                        tctcTTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAAc          <  1:2149560/149‑1 (MQ=255)

                                                                                                                                         ctctTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACt         >  1:3388466/1‑149 (MQ=255)

                                                                                                                                           ctTTTTAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGc       >  1:152877/1‑149 (MQ=255)

                                                                                                                                            tttttAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCa      >  2:3928138/1‑149 (MQ=255)

                                                                                                                                             ttttAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGGTAAAGTGCTAAGCGCAGCTTGGTAGCCAGGGTTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCAt     <  2:4440446/149‑1 (MQ=255)

                                                                                                                                             ttttAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCAt     >  1:8721481/1‑149 (MQ=255)

                                                                                                                                             ttttAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCAt     >  1:254697/1‑149 (MQ=255)

                                                                                                                                              tttAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCATc    >  1:8513387/1‑149 (MQ=255)

                                                                                                                                               ttAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCATCa   <  1:6459438/149‑1 (MQ=255)

                                                                                                                                                tAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCATCaa  >  2:7830185/1‑149 (MQ=255)

                                                                                                                                                tAGAGAA‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCATCaa  >  2:7380152/1‑149 (MQ=255)

TTCAGTCAAAATTTCTCGAAAACTATCTGTTCTCTCCTGAGTCATCCATGTATTGGACGAACCACCAATGAAAGCTAATTTCTGATAGCCTTTTTGAATCATGTGAGAAACAAGTAGTTGAGTGGGAGTGCGTTGGTCTCTTTTTAGAG‑A‑AAA‑ATATTTGCTGACAGATGAGTTAAATCACGTCCAAGTAGAATAAGGGGAATTTCTTTGGCTAAATTGCTAAGAGCAGCTTGGTAGTCAGGATTAATATCAGTTTCATCTAATAATCCTCCTAAAATTAAAACTGCATCAA  >  AM406671/1210130‑1210421

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