BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	698,817	T→C	5.3%	M121V (ATG→GTG)	asnB ←	asparagine synthetase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	698,817	0	T	C	5.3%	2535.4 / 83.4	766	M121V (ATG→GTG)	asnB	asparagine synthetase B
Reads supporting (aligned to +/- strand): ref base T (381/344); new base C (24/17); total (405/361)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

Only 100 of 770 total aligned reads displayed.

CGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCACCGTTTACCGCCAGT  >  NC_000913/698669‑698965

 ggCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATg                                                                                                                                                     >  4:2129880/1‑149 (MQ=255)

  gCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATTAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGc                                                                                                                                                    >  1:5354649/1‑149 (MQ=255)

   cAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGcc                                                                                                                                                   <  3:2701035/149‑1 (MQ=255)

   cAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGcc                                                                                                                                                   <  2:768358/149‑1 (MQ=255)

   cAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGcc                                                                                                                                                   <  2:2487039/149‑1 (MQ=255)

       cTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGc                                                                                                                                               >  3:2147826/1‑149 (MQ=255)

         ggCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCaa                                                                                                                                             >  3:6542300/1‑149 (MQ=255)

         ggCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCaa                                                                                                                                             <  1:7201581/149‑1 (MQ=255)

             ccAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtcg                                                                                                                                         >  3:3762443/1‑149 (MQ=255)

               aGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtcgtc                                                                                                                                       <  4:6252777/149‑1 (MQ=255)

               aGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtcgtc                                                                                                                                       <  1:6662167/149‑1 (MQ=255)

                   cTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtcgtc                                                                                                                                       <  4:7061742/145‑1 (MQ=255)

                   cTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtcgtc                                                                                                                                       >  3:7061742/1‑145 (MQ=255)

                   cTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCgaga                                                                                                                                   >  4:2432590/1‑149 (MQ=255)

                    tttCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCaa                                                                                                                                             >  2:3414401/1‑138 (MQ=255)

                    tttCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGaa                                                                                                                                  <  1:1985696/149‑1 (MQ=255)

                        aTTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTc                                                                                                                              >  2:175146/1‑149 (MQ=255)

                         tttCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTcc                                                                                                                             >  1:3624948/1‑149 (MQ=255)

                             tGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGcc                                                                                                                         <  1:391502/149‑1 (MQ=255)

                              gAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGccc                                                                                                                        >  4:3524699/1‑149 (MQ=255)

                              gAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGccc                                                                                                                        >  2:3761483/1‑149 (MQ=255)

                               aGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCt                                                                                                                       <  4:141362/149‑1 (MQ=255)

                               aGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGAAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCt                                                                                                                       <  1:6809825/149‑1 (MQ=255)

                                    acaTACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtcgt                                                                                                                                        >  4:6098926/1‑127 (MQ=255)

                                      aTACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCt                                                                                                                >  1:5289648/1‑149 (MQ=255)

                                      aTACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCt                                                                                                                >  2:6536250/1‑149 (MQ=255)

                                        aCAGCTGACCGTGTTCGTCCTACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGa                                                                                                              <  2:4192777/149‑1 (MQ=255)

                                        aCAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGa                                                                                                              >  1:5670777/1‑149 (MQ=255)

                                          aGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATa                                                                                                            >  4:1192738/1‑149 (MQ=255)

                                           gCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGAACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATAc                                                                                                           >  4:1701324/1‑149 (MQ=255)

                                              gACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAgc                                                                                                        >  2:3929063/1‑149 (MQ=255)

                                               aCCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAgcg                                                                                                       >  2:317740/1‑149 (MQ=255)

                                                ccGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGc                                                                                                                          >  2:5323076/1‑129 (MQ=255)

                                                 cGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGcgcg                                                                                                     >  1:3331324/1‑149 (MQ=255)

                                                  gtgtTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCt                                                                                                                >  4:2993742/1‑137 (MQ=255)

                                                  gtgtTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGa                                                                                                    <  4:5097489/149‑1 (MQ=255)

                                                     ttCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGa                                                                                                 >  2:4351303/1‑149 (MQ=255)

                                                        gtccttccccATATACAGTGGGATGATCCCCAGATGGTCGCCACCAATCAGGTAGCCATCTTTTTCGCTGTCGTACAGTGCTAAGGCAAACATTCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGTTACAGCGCGAGGATCa                                                                                              <  2:5417029/143‑1 (MQ=255)

                                                        gTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCa                                                                                              >  1:1423416/1‑149 (MQ=255)

                                                        gTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACACGCCCTGCAAGTCGTCGGGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCa                                                                                              >  4:1880611/1‑149 (MQ=255)

                                                        gTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACACGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCa                                                                                              >  2:2808287/1‑149 (MQ=255)

                                                          cATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACt                                                                                            <  1:4295248/149‑1 (MQ=255)

                                                            tACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTc                                                                                          <  3:4221448/149‑1 (MQ=255)

                                                                ccATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtc                                                                                                                                          <  3:6361596/97‑1 (MQ=255)

                                                                ccATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAgtc                                                                                                                                          <  3:2951160/97‑1 (MQ=255)

                                                                ccAGATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACACGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCt                                                                                                                >  2:1182072/1‑123 (MQ=255)

                                                                 cATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGt                                                                                     <  4:2824965/149‑1 (MQ=255)

                                                                   tataCAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCa                                                                                   <  2:3001376/149‑1 (MQ=255)

                                                                   tataCAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCa                                                                                   <  3:4959059/149‑1 (MQ=255)

                                                                     taCAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGa                                                                                 <  4:5450318/149‑1 (MQ=255)

                                                                        aGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGAccc                                                                              >  4:4245871/1‑149 (MQ=255)

                                                                         gTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCt                                                                                                                <  4:4910499/114‑1 (MQ=255)

                                                                            ggATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTc                                                                          <  4:3869883/149‑1 (MQ=255)

                                                                            ggATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTc                                                                          <  3:3914406/149‑1 (MQ=255)

                                                                            ggATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTc                                                                          >  3:1947908/1‑149 (MQ=255)

                                                                             gatgatCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAg                                                                                                   >  1:893439/1‑123 (MQ=255)

                                                                             gatgatCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCt                                                                         <  4:1012593/149‑1 (MQ=255)

                                                                               tgatCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCttt                                                                                                                     >  2:1482620/1‑103 (MQ=255)

                                                                                  tCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAAc                                                                    >  3:1220246/1‑149 (MQ=255)

                                                                                    cccAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTg                                                                  <  3:5138799/149‑1 (MQ=255)

                                                                                     ccAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTgg                                                                 <  4:5595782/149‑1 (MQ=255)

                                                                                     ccAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTgg                                                                 >  1:4727681/1‑149 (MQ=255)

                                                                                     ccAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTgg                                                                 >  1:1217728/1‑149 (MQ=255)

                                                                                       aGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCCGTCTGGAACTGGTa                                                               <  2:4705426/149‑1 (MQ=255)

                                                                                         aTGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAAc                                                             <  3:238908/149‑1 (MQ=255)

                                                                                            gTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGAt                                                          <  4:5771448/149‑1 (MQ=255)

                                                                                            gTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGAt                                                          <  3:4366347/149‑1 (MQ=255)

                                                                                             tCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATc                                                         <  4:4538197/149‑1 (MQ=255)

                                                                                                cgACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCGCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTAATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGcc                                                      <  2:2242327/149‑1 (MQ=255)

                                                                                                  aCCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAGGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCAGTCTGGAACTGGTAACGATCGCCat                                                    >  3:4659109/1‑149 (MQ=255)

                                                                                                     aaTCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATAtt                                                 >  1:6195462/1‑149 (MQ=255)

                                                                                                       tCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATa                                                                                                            >  1:2749453/1‑88 (MQ=255)

                                                                                                         aGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGgc                                             >  3:3452503/1‑149 (MQ=255)

                                                                                                          ggTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTcc                                                                                                                 <  4:6494949/80‑1 (MQ=255)

                                                                                                           gTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGgcgc                                           >  4:6394771/1‑149 (MQ=255)

                                                                                                            tAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGgcgcg                                          >  4:3168979/1‑149 (MQ=255)

                                                                                                               gCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCaa                                       >  2:5616087/1‑149 (MQ=255)

                                                                                                                 aTCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATg                                     >  1:4175601/1‑149 (MQ=255)

                                                                                                                   cTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGcc                                   <  4:6816447/149‑1 (MQ=255)

                                                                                                                    tttttCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCt                                  >  2:7035009/1‑149 (MQ=255)

                                                                                                                    tttttCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCt                                  <  3:6149803/149‑1 (MQ=255)

                                                                                                                      tttCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCtgg                                <  3:5572704/149‑1 (MQ=255)

                                                                                                                       ttCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCtggt                               >  1:1849257/1‑149 (MQ=255)

                                                                                                                       ttCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCtggt                               <  2:6259801/149‑1 (MQ=255)

                                                                                                                          gCTGTCGTACAGTGCAAAGGCAATCATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAGCTGGTAACGATCGCCATAGTCGGCGCGCAACGCCTggtggt                            <  4:5296644/149‑1 (MQ=255)

                                                                                                                            tGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAAc                                                                    <  1:2368047/107‑1 (MQ=255)

                                                                                                                            tGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTg                          >  3:1596294/1‑149 (MQ=255)

                                                                                                                             gTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGt                         >  1:1057799/1‑149 (MQ=255)

                                                                                                                              tCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCTTGGTGGTTGTa                        <  1:3802243/149‑1 (MQ=255)

                                                                                                                               cGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGa                                                                                                    <  3:101524/72‑1 (MQ=255)

                                                                                                                               cGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAg                       <  3:3146243/149‑1 (MQ=255)

                                                                                                                                   cAGTGCAAAGGCAAACACGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGAttt                   >  3:4819617/1‑149 (MQ=255)

                                                                                                                                      tGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTATGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCAc                >  3:6741582/1‑149 (MQ=255)

                                                                                                                                       gCAAAGGCAAAAATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCAcc               <  2:5033082/149‑1 (MQ=255)

                                                                                                                                        cAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCt                                  <  1:5050468/129‑1 (MQ=255)

                                                                                                                                         aaaGGCAAACATGCCCTGCGAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCACCGt             >  3:5890465/1‑149 (MQ=255)

                                                                                                                                         aaaGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCAAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCACCGt             >  4:2648314/1‑149 (MQ=255)

                                                                                                                                            ggCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCt                                                                         >  3:576392/1‑86 (MQ=255)

                                                                                                                                                 aCATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCACCGTTTACCGcc     <  1:2984378/149‑1 (MQ=255)

                                                                                                                                                  cATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCACCGTTTACCGCCa    <  3:5418362/149‑1 (MQ=255)

CGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGCAAACATGCCCTGCAAGTCGTCGAGAAATTCCGGCCCTTTTTCCTGATACAGCGCGAGGATCACTTCACAGTCAGACCCGGTCTGGAACTGGTAACGATCGCCATATTCGGCGCGCAATGCCTGGTGGTTGTAGATTTCACCGTTTACCGCCAGT  >  NC_000913/698669‑698965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`