BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,241,716	C→G	6.4%	C180S (TGC→TCC) ‡	xylE ←	D‑xylose transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,241,716	0	C	G	6.4%	441.5 / 16.8	125	C180S (TGC→TCC) ‡	xylE	D‑xylose transporter
Reads supporting (aligned to +/- strand): ref base C (56/61); new base G (3/5); total (59/66)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 125 total aligned reads displayed.

CATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAA  >  NC_000913/4241571‑4241863

cATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGt                                                                                                                                                  >  3:8910/1‑149 (MQ=255)

  tCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTaa                                                                                                                                                >  1:387770/1‑149 (MQ=255)

   cagccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTaaa                                                                                                                                               <  1:1146226/149‑1 (MQ=255)

    agccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAAc                                                                                                                                              >  2:1058540/1‑149 (MQ=255)

       cagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTat                                                                                                                                           <  1:857239/149‑2 (MQ=255)

       cagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTaa                                                                                                                                           >  1:78487/1‑149 (MQ=255)

        agcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTatt                                                                                                                                          <  3:134016/149‑3 (MQ=255)

        agcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTaaa                                                                                                                                          >  2:101876/1‑149 (MQ=255)

         gcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                         <  1:763865/149‑4 (MQ=255)

             ggACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTg                                                                                                                                     <  4:569820/149‑1 (MQ=255)

             ggACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTg                                                                                                                                     <  1:101876/149‑1 (MQ=255)

              gACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGc                                                                                                                                    >  1:822297/1‑149 (MQ=255)

              gACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGc                                                                                                                                    <  1:1013283/149‑1 (MQ=255)

                 tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTaaaa                                                                                                                                         >  3:417585/1‑141 (MQ=255)

                 tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTaaaa                                                                                                                                         <  4:417585/141‑1 (MQ=255)

                 tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCg                                                                                                                                 >  3:944674/1‑149 (MQ=255)

                 tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCg                                                                                                                                 >  4:950051/1‑149 (MQ=255)

                 tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCg                                                                                                                                 >  2:351878/1‑149 (MQ=255)

                 tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCg                                                                                                                                 >  1:391946/1‑149 (MQ=255)

                     tGGCACGGTATACAGCAGCATTAAGAACGGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGaaaa                                                                                                                             <  2:773138/149‑1 (MQ=255)

                        cACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAaataa                                                                                                                          <  3:364009/149‑1 (MQ=255)

                         aCGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAaataat                                                                                                                         >  2:954940/1‑149 (MQ=255)

                          cGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATc                                                                                                                        <  2:78487/149‑1 (MQ=255)

                             tataCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCa                                                                                                                     >  3:895052/1‑149 (MQ=255)

                                aCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGGCAGCTGGCATCACCGGAACGGGCAATAGAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAAc                                                                                                                  >  3:969394/1‑149 (MQ=255)

                                    cagcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGt                                                                                                              <  1:550316/149‑1 (MQ=255)

                                    cagcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGt                                                                                                              <  1:351878/149‑1 (MQ=255)

                                      gcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                         <  3:773199/120‑4 (MQ=255)

                                      gcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                         >  4:773199/1‑117 (MQ=255)

                                      gcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGc                                                                                                                                    <  3:706571/125‑1 (MQ=255)

                                      gcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGc                                                                                                                                    >  4:706571/1‑125 (MQ=255)

                                        aTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGc                                                                                                                      >  1:545488/1‑137 (MQ=255)

                                        aTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTaaa                                                                                                          >  2:456638/1‑149 (MQ=255)

                                           aaGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAaga                                                                                                       <  2:299004/149‑1 (MQ=255)

                                           aaGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAaga                                                                                                       <  2:7131/149‑1 (MQ=255)

                                            aGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAagag                                                                                                      <  2:391946/149‑1 (MQ=255)

                                            aGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAagag                                                                                                      <  2:1088772/149‑1 (MQ=255)

                                                  gcagTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGt                                                                                                >  4:700833/1‑149 (MQ=255)

                                                       gCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTccc                                                                                           <  2:577851/149‑1 (MQ=255)

                                                         aGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCcgc                                                                                         >  4:670539/1‑149 (MQ=255)

                                                         aGGGATACATGCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCcgc                                                                                         >  3:776229/1‑149 (MQ=255)

                                                           ggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                       >  3:826462/1‑149 (MQ=255)

                                                           ggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                       >  1:237657/1‑149 (MQ=255)

                                                           ggATACATTCCGAGCCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCCCCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAAACGCAAACTGGTTAAAAGAGACCAGTTTCCCCCGa                                                                                       >  2:965659/1‑149 (MQ=255)

                                                             aTACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTc                                                                                             >  4:816755/1‑141 (MQ=255)

                                                               aCATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAatat                                                                                   <  4:895052/149‑1 (MQ=255)

                                                                 aTTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGa                                                                                 >  4:1105915/1‑149 (MQ=255)

                                                                  ttCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAagag                                                                                                      >  3:617246/1‑127 (MQ=255)

                                                                   tCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGc                                                                               <  3:182671/149‑1 (MQ=255)

                                                                      gAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTgg                                                                            >  4:367479/1‑149 (MQ=255)

                                                                           aaaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCt                                                                              <  4:1031469/142‑1 (MQ=255)

                                                                           aaaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCt                                                                              >  3:1031469/1‑142 (MQ=255)

                                                                           aaaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCa                                                                       >  1:163063/1‑149 (MQ=255)

                                                                           aaaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCa                                                                       >  2:112201/1‑149 (MQ=255)

                                                                                      gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                         <  1:1075467/72‑4 (MQ=38)

                                                                                      gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                         >  2:1075467/1‑69 (MQ=38)

                                                                                      gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAata                                                            >  3:580773/1‑149 (MQ=255)

                                                                                      gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAata                                                            <  3:1105915/149‑1 (MQ=255)

                                                                                       ccAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAatat                                                           <  3:225275/149‑1 (MQ=255)

                                                                                         aGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATAc                                                         <  1:456638/149‑1 (MQ=255)

                                                                                         aGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATAc                                                         <  3:670539/149‑1 (MQ=255)

                                                                                          gCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTaaaa                                                                                                         <  3:259131/100‑1 (MQ=255)

                                                                                          gCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTaaaa                                                                                                         >  4:259131/1‑100 (MQ=255)

                                                                                           ccGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACAt                                                       <  3:623448/149‑1 (MQ=255)

                                                                                            cGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACAtt                                                      <  2:172907/149‑1 (MQ=255)

                                                                                             gTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTg                                                     <  3:700833/149‑1 (MQ=255)

                                                                                              tCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATAc                                                         <  1:53986/144‑1 (MQ=255)

                                                                                              tCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATAc                                                         >  2:53986/1‑144 (MQ=255)

                                                                                              tCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTgg                                                    >  3:66740/1‑149 (MQ=255)

                                                                                               cAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGc                                                   <  3:367479/149‑1 (MQ=255)

                                                                                                      cagccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCa                                            <  4:776229/149‑1 (MQ=255)

                                                                                                       agccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTc                                                                   <  4:213196/125‑1 (MQ=255)

                                                                                                       agccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTc                                                                   >  3:213196/1‑125 (MQ=255)

                                                                                                       agccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAt                                           >  3:833965/1‑149 (MQ=255)

                                                                                                        gccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGa                                                                           >  1:1001189/1‑116 (MQ=255)

                                                                                                        gccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGa                                                                           <  2:1001189/116‑1 (MQ=255)

                                                                                                        gccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAtt                                          <  3:950051/149‑1 (MQ=255)

                                                                                                          cagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGa                                        <  4:944674/149‑1 (MQ=255)

                                                                                                                gCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTaa                                  >  4:656829/1‑149 (MQ=255)

                                                                                                                   tCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAAcc                               >  3:454453/1‑149 (MQ=255)

                                                                                                                     aCCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCaa                             >  4:1111492/1‑149 (MQ=255)

                                                                                                                      ccGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAAc                            <  4:451392/149‑1 (MQ=255)

                                                                                                                      ccGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAAc                            <  3:1111492/149‑1 (MQ=255)

                                                                                                                       cGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACg                           >  2:489892/1‑149 (MQ=255)

                                                                                                                             gggCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATa                     >  4:1138363/1‑149 (MQ=255)

                                                                                                                              ggCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATAc                    <  4:826462/149‑1 (MQ=255)

                                                                                                                              ggCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATAc                    >  3:1023867/1‑149 (MQ=255)

                                                                                                                                cAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCg                  >  4:275099/1‑149 (MQ=255)

                                                                                                                                 aaTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATa                     >  4:606485/1‑145 (MQ=255)

                                                                                                                                 aaTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATa                     <  3:606485/145‑1 (MQ=255)

                                                                                                                                 aaTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGc                 >  1:1061726/1‑149 (MQ=255)

                                                                                                                                 aaTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGc                 >  2:153882/1‑149 (MQ=255)

                                                                                                                                   tAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTcc                                                                  >  4:556170/1‑98 (MQ=255)

                                                                                                                                   tAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTcc                                                                  <  3:556170/98‑1 (MQ=255)

                                                                                                                                    aaaaTAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaa              <  3:45010/149‑1 (MQ=255)

                                                                                                                                            ttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGa      <  3:1138363/149‑1 (MQ=255)

                                                                                                                                            ttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGa      <  1:1066097/149‑1 (MQ=255)

                                                                                                                                             tACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGat     <  2:1061726/149‑1 (MQ=255)

                                                                                                                                                gCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGataaa  <  2:237657/149‑1 (MQ=255)

                                                                                                                                                gCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGataaa  <  4:919669/149‑1 (MQ=255)

CATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAA  >  NC_000913/4241571‑4241863

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