BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,241,715	G→T	6.9%	C180* (TGC→TGA) ‡	xylE ←	D‑xylose transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,241,715	0	G	T	6.9%	342.8 / 11.9	116	C180* (TGC→TGA) ‡	xylE	D‑xylose transporter
Reads supporting (aligned to +/- strand): ref base G (59/49); new base T (5/3); total (64/52)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 116 total aligned reads displayed.

GCGACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAA  >  NC_000913/4241567‑4241862

gcgACATCAGCCAGGGAGGACTTTCTGGCACGGTATACAGCAGCACTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGACCGGGCAATAAAATCGTTTACGc                                                                                                                                                    <  2:827602/150‑1 (MQ=255)

 cgACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCa                                                                                                                                                   <  2:10138/150‑1 (MQ=255)

   aCATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGt                                                                                                                                                 >  1:1017237/1‑150 (MQ=255)

     aTCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTTCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTaa                                                                                                                                               <  1:578044/150‑1 (MQ=255)

     aTCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGTAACGGGCAATAAAATAGTTTACGCAGTaa                                                                                                                                               <  2:663185/150‑1 (MQ=255)

       cagccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCCGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAAc                                                                                                                                             >  1:882042/1‑150 (MQ=255)

       cagccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTATAc                                                                                                                                             >  2:338982/1‑150 (MQ=255)

          ccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTaa                                                                                                                                          >  2:1413155/1‑150 (MQ=255)

            agcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                        <  2:524047/150‑4 (MQ=255)

             gcGAGGACTTTCTGGGTCTGTATACAGGAGCATTAAGAACAGCATTGCAGGGATACATTCCGAGGCAAACATATACCGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAg                                                                                                                                       <  2:318600/150‑1 (MQ=255)

             gcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAg                                                                                                                                       <  1:1447316/150‑1 (MQ=255)

             gcGAGGACTTTCTGGCACGGTAGACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCAAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTatttt                                                                                                                                       <  2:1083111/150‑5 (MQ=255)

               gAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattttat                                                                                                                                     <  1:133614/150‑7 (MQ=255)

                 ggACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGc                                                                                                                                   >  2:224058/1‑150 (MQ=255)

                 ggACTATCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGc                                                                                                                                   <  1:1568570/150‑1 (MQ=255)

                  gACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGcc                                                                                                                                  <  1:1209515/150‑1 (MQ=255)

                             aCGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATc                                                                                                                       <  1:1643386/150‑1 (MQ=255)

                              cGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCg                                                                                                                      <  2:553988/150‑1 (MQ=255)

                               ggTATACAGCAGCATTTAGAACAGCAGTGCAGGGATACATTCCGAGTCAAACATCTAACGCCAGCAGTCATTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAATATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGc                                                                                                                     <  1:170195/150‑1 (MQ=255)

                                gTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCa                                                                                                                    <  1:1223816/150‑1 (MQ=255)

                                  ataCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCaaa                                                                                                                  >  1:1708009/1‑150 (MQ=255)

                                   taCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAAc                                                                                                                 <  2:193368/150‑1 (MQ=255)

                                    aCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCCAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACt                                                                                                                >  1:1117337/1‑150 (MQ=255)

                                     cagcagCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCATCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTg                                                                                                               <  2:121511/150‑1 (MQ=255)

                                     cagcagCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTg                                                                                                               <  1:421852/150‑1 (MQ=255)

                                       gcagcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGt                                                                                                             <  2:87528/150‑1 (MQ=255)

                                       gcagcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGt                                                                                                             >  2:735291/1‑150 (MQ=255)

                                        ccccactAACAACAGCAGCGCACGGCTACCTTCCGACGCAAACATATACCCCCACCCGTCAGTACTCACCCCGCTCGCATCACCGGAACGGGCAATAAAATAGTTAACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGtt                                                                                                            <  2:1132716/144‑1 (MQ=255)

                                        cagcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGtt                                                                                                            <  2:1229010/150‑1 (MQ=255)

                                             ttAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAag                                                                                                       <  2:724845/150‑1 (MQ=255)

                                             ttAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAag                                                                                                       <  2:895875/150‑1 (MQ=255)

                                             ttAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGAAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAag                                                                                                       <  1:113089/150‑1 (MQ=255)

                                                 gAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGAc                                                                                                   >  1:855238/1‑150 (MQ=255)

                                                 gAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGAc                                                                                                   >  2:807849/1‑150 (MQ=255)

                                                 gAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATAGCAAACTGGTTAAAAGAGAc                                                                                                   >  1:860676/1‑150 (MQ=255)

                                                  aaCAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGAcc                                                                                                  <  1:208782/150‑1 (MQ=255)

                                                  aaCAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGAcc                                                                                                  <  1:1476222/150‑1 (MQ=255)

                                                   aCAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCa                                                                                                 >  2:1448717/1‑150 (MQ=255)

                                                     agcagTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGt                                                                                               <  1:245061/150‑1 (MQ=255)

                                                              gggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                      >  1:147523/1‑150 (MQ=255)

                                                              gggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                      >  1:671099/1‑150 (MQ=255)

                                                              gggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                      >  1:170402/1‑150 (MQ=255)

                                                              gggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                      >  2:1481072/1‑150 (MQ=255)

                                                                 aTACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAATGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAata                                                                                   >  2:1365385/1‑150 (MQ=255)

                                                                    cATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGa                                                                                >  1:771377/1‑150 (MQ=255)

                                                                     aTTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAg                                                                               <  2:742058/150‑1 (MQ=255)

                                                                       tCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAATAGAGACCAGTTTCCCGCGAATATGAGCt                                                                             >  2:195484/1‑150 (MQ=255)

                                                                       tCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCt                                                                             <  2:533760/150‑1 (MQ=255)

                                                                        ccGAGGCAAACTTATAACGCCAGCCGGCACTTTTCCGCCAGATGGCATCACCGGAAATGGCAAATCAATAGTTTACGCAGTAAACTAAAAGGTGCCCGAAAATAATCGCAAACTGGGTAAACGAGACCAGTTTCCCGCGAATATGACCTg                                                                            <  2:41976/150‑1 (MQ=255)

                                                                         cGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCATCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTgg                                                                           >  2:1552803/1‑150 (MQ=255)

                                                                          gAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGa                                                                          >  2:1341737/1‑150 (MQ=255)

                                                                             gCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGATTATGAGCTGGAGcc                                                                       >  1:1069034/1‑150 (MQ=255)

                                                                             gCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGcc                                                                       >  1:732661/1‑150 (MQ=255)

                                                                              cAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCa                                                                      >  2:225230/1‑150 (MQ=255)

                                                                               aaaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAg                                                                     <  2:54067/150‑1 (MQ=255)

                                                                                aaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGt                                                                    >  1:1464424/1‑150 (MQ=255)

                                                                                  cATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTc                                                                  >  1:988702/1‑150 (MQ=255)

                                                                                    tataACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCGATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCg                                                                >  1:1072000/1‑150 (MQ=255)

                                                                                    tataACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCg                                                                <  1:511117/150‑1 (MQ=255)

                                                                                       aaCGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCaa                                                             >  2:164901/1‑150 (MQ=255)

                                                                                       aaCGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCaa                                                             >  2:1061249/1‑150 (MQ=255)

                                                                                          gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAatat                                                          <  2:1647717/150‑1 (MQ=255)

                                                                                              gCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACAt                                                      >  2:498165/1‑150 (MQ=255)

                                                                                                 gccAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTgg                                                   >  2:1063134/3‑150 (MQ=255)

                                                                                                 gTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTgg                                                   >  2:1373471/1‑150 (MQ=255)

                                                                                                 gTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTgg                                                   <  1:1522857/150‑1 (MQ=255)

                                                                                                    aGTATTCCGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCga                                                >  1:281438/1‑150 (MQ=255)

                                                                                                     gTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCgag                                               >  1:496963/1‑150 (MQ=255)

                                                                                                       aTTCAGCCAGCTGGCATCACCTGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACAAGTTTAACGCGAATATGAGAATCAGCCAGTTCCGCAATATACAATGgagagaa                                             >  1:333667/1‑149 (MQ=255)

                                                                                                       aTTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGagag                                             <  1:1639224/150‑1 (MQ=255)

                                                                                                         tCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCG‑AAATAATACCAAACTTGTTAAAAAAGACCCGTTTACCGCGAATATTAACTTGAGCCAGTTACGCAATATACATTTGCGACAGCaa                                          >  2:1713094/1‑149 (MQ=255)

                                                                                                          cagccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTATAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAt                                          >  1:1642991/1‑150 (MQ=255)

                                                                                                            gccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTg                                        >  1:826150/1‑150 (MQ=255)

                                                                                                            gccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTg                                        >  1:562036/1‑150 (MQ=255)

                                                                                                              cagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAg                                      >  1:972133/1‑150 (MQ=255)

                                                                                                               agcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAgg                                     >  2:990544/1‑150 (MQ=255)

                                                                                                               agcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAgg                                     >  2:333147/1‑150 (MQ=255)

                                                                                                                gcTTGCATCCCAGGAACGGTCAATAAAATAGTTTACGCATTAAACTAAAAGTTGCCCGAATATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCTCGAATATGAGCTGGAGCCAGTTCCGTAATATACATTGGCGAGAGCATTGAGGc                                    <  2:1204789/150‑1 (MQ=255)

                                                                                                                   ggCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTaa                                 >  1:898160/1‑150 (MQ=255)

                                                                                                                    gCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTTGCGAGAGCAGTGAGGCTaaa                                >  1:95216/1‑150 (MQ=255)

                                                                                                                           cGGAACGGGCAATAAAATAGTTTACGCACTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGc                         <  1:1328624/150‑1 (MQ=255)

                                                                                                                            ggAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGcc                        <  2:1708009/150‑1 (MQ=255)

                                                                                                                             gAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTg                                                                            >  2:88323/1‑97 (MQ=255)

                                                                                                                                  ggCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATAcc                  >  2:1618759/1‑150 (MQ=255)

                                                                                                                                   gCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCg                 <  1:1341737/150‑1 (MQ=255)

                                                                                                                                     aaTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGcc               <  2:29639/150‑1 (MQ=255)

                                                                                                                                      aTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCa              <  2:935056/150‑1 (MQ=255)

                                                                                                                                       acttttAGTTTACTGCGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaa             >  2:805702/6‑150 (MQ=255)

                                                                                                                                        cttttAGTTTACTGCGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaat            >  1:918303/5‑150 (MQ=255)

                                                                                                                                         ttttAGTTTACTGCGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaata           >  1:365684/4‑150 (MQ=255)

                                                                                                                                         ttttAGTTTACTGCGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaata           >  1:1021169/4‑150 (MQ=255)

                                                                                                                                          tttAGTTTACTGCGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaataa          >  1:220689/3‑150 (MQ=255)

                                                                                                                                          aaTAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaataa          <  2:196787/150‑1 (MQ=255)

                                                                                                                                            tAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATg        <  2:42143/150‑1 (MQ=255)

                                                                                                                                             aGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGc       >  2:217499/1‑150 (MQ=255)

                                                                                                                                             aGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGc       >  1:1004601/1‑150 (MQ=255)

                                                                                                                                                ttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGGCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGat    <  1:131102/150‑1 (MQ=255)

                                                                                                                                                ttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTACACCAACGCCAATACCGCCAATAATGCGat    >  1:313147/1‑150 (MQ=255)

                                                                                                                                                ttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGat    <  2:771377/150‑1 (MQ=255)

                                                                                                                                                  aCGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGataa  >  2:568908/1‑150 (MQ=255)

GCGACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAA  >  NC_000913/4241567‑4241862

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`