BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,241,715	0	G	T	5.4%	377.0 / 7.5	111	C180* (TGC→TGA) ‡	xylE	D‑xylose:H(+) symporter
Reads supporting (aligned to +/- strand): ref base G (56/49); new base T (4/2); total (60/51)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 112 total aligned reads displayed.

GCGACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAAT  >  NC_000913/4241567‑4241864

gcgACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGc                                                                                                                                                      <  1:26374/150‑1 (MQ=255)

gcgACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGc                                                                                                                                                      <  2:1343202/150‑1 (MQ=255)

gcgACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGc                                                                                                                                                      >  2:351762/1‑150 (MQ=255)

 cgACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCa                                                                                                                                                     >  2:955105/1‑150 (MQ=255)

    cATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTa                                                                                                                                                  <  1:520687/150‑1 (MQ=255)

    cATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTa                                                                                                                                                  >  1:289201/1‑150 (MQ=255)

      tCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTaaa                                                                                                                                                >  1:1249973/1‑150 (MQ=255)

      tCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTaaa                                                                                                                                                >  1:1277689/1‑150 (MQ=255)

       cagccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAAc                                                                                                                                               >  2:1536804/1‑150 (MQ=255)

       cagccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAAc                                                                                                                                               >  2:1175154/1‑150 (MQ=255)

        agccagcGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACt                                                                                                                                              <  1:1536804/150‑1 (MQ=255)

               gAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGAAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACTGGCAATAAAATAGTTTACGCAGTAAACTAAAAGtt                                                                                                                                       >  2:1549536/1‑150 (MQ=255)

                aGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                          >  2:560408/1‑143 (MQ=255)

                aGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACTGCGTAAACTattt                                                                                                                                          <  1:560408/146‑4 (MQ=255)

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

                     tttCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGa                                                                                                                                 <  2:418095/150‑1 (MQ=255)

                           gCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAaataa                                                                                                                           >  1:1128634/1‑150 (MQ=255)

                              cGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATATTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCg                                                                                                                        >  2:1300276/1‑150 (MQ=255)

                                   taCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAAc                                                                                                                   >  1:1393443/1‑150 (MQ=255)

                                    aCAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGATGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACt                                                                                                                  >  2:257202/1‑150 (MQ=255)

                                      agcagcATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTgg                                                                                                                <  2:1019008/150‑1 (MQ=255)

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

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

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

                                         agcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTa                                                                                                             <  2:1157062/150‑1 (MQ=255)

                                         agcaTTAAGAACAGCAGTGCAGGGATAAATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTAAGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTa                                                                                                             >  2:560500/1‑150 (MQ=255)

                                          gcaTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGccc                                                                                                                                   <  1:1216323/127‑1 (MQ=255)

                                            aTTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAAT‑AAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAag                                                                                                         <  1:986607/150‑1 (MQ=255)

                                             ttAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAag                                                                                                         >  2:922918/1‑150 (MQ=255)

                                             ttAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAAAATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAag                                                                                                         <  2:711320/150‑1 (MQ=255)

                                               aaGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAagag                                                                                                       <  2:1463709/150‑1 (MQ=255)

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

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

                                                      gcagTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGtt                                                                                                >  2:1585644/1‑150 (MQ=255)

                                                         gTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTcc                                                                                             >  1:188825/1‑150 (MQ=255)

                                                             aGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCcgcg                                                                                         <  2:1500989/150‑1 (MQ=255)

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

                                                              gggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                        <  1:380050/150‑1 (MQ=255)

                                                              gggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                        <  1:1175154/150‑1 (MQ=255)

                                                                gATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAat                                                                                      <  2:869530/150‑1 (MQ=255)

                                                                gATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACT‑AAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAata                                                                                     <  2:134037/150‑1 (MQ=255)

                                                                   aCATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATg                                                                                   >  2:1036554/1‑150 (MQ=255)

                                                                   aCATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACCAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATg                                                                                   <  2:1059574/150‑1 (MQ=255)

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

                                                                    cATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGa                                                                                  <  2:664021/150‑1 (MQ=255)

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

                                                                          gAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGa                                                                            <  1:351762/150‑1 (MQ=255)

                                                                           aGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAg                                                                           <  1:1588617/150‑1 (MQ=255)

                                                                             gCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGcc                                                                         <  2:27665/150‑1 (MQ=255)

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

                                                                                 aCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGtt                                                                     <  2:579349/150‑1 (MQ=255)

                                                                                   atatAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATATTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTcc                                                                   <  1:1300276/150‑1 (MQ=255)

                                                                                    tataACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCg                                                                  >  2:242746/1‑150 (MQ=255)

                                                                                    tataACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCg                                                                  >  1:242426/1‑150 (MQ=255)

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

                                                                                      taACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCa                                                                >  1:700235/1‑150 (MQ=255)

                                                                                       aaCGCCAGCAGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCaa                                                               <  1:688249/150‑1 (MQ=255)

                                                                                        aCGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCagt                                                              >  1:77171/1‑148 (MQ=255)

                                                                                        aCGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAat                                                              <  1:242746/150‑1 (MQ=255)

                                                                                        aCGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAat                                                              <  1:777768/150‑1 (MQ=255)

                                                                                         cGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGGGAATATGAGCTGGAGCCAGTTCCGCAata                                                             >  2:325127/1‑150 (MQ=255)

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

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

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

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

                                                                                                   cAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCg                                                   <  2:188825/150‑1 (MQ=255)

                                                                                                    aGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCga                                                  <  2:1169762/150‑1 (MQ=255)

                                                                                                    aGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCga                                                  >  2:1306796/1‑150 (MQ=255)

                                                                                                     gTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCgag                                                 <  1:955105/150‑1 (MQ=255)

                                                                                                     gTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCgag                                                 <  2:1125910/150‑1 (MQ=255)

                                                                                                      tATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAGTATACATTGGCgaga                                                <  2:77171/150‑1 (MQ=255)

                                                                                                       aTTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGagag                                               >  1:1202057/1‑150 (MQ=255)

                                                                                                        ttCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGc                                              <  2:1241139/150‑1 (MQ=255)

                                                                                                         tCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCa                                             >  2:929851/1‑150 (MQ=255)

                                                                                                          cagccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAt                                            >  1:1518803/1‑150 (MQ=255)

                                                                                                          cagccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAt                                            >  1:890574/1‑150 (MQ=255)

                                                                                                           agccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAtt                                           <  1:563506/150‑1 (MQ=255)

                                                                                                                gcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACAtt                                                       <  1:1093009/133‑1 (MQ=255)

                                                                                                                gcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACAtt                                                       >  2:1093009/1‑133 (MQ=255)

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

                                                                                                                     cATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAAc                                 <  1:1036554/150‑1 (MQ=255)

                                                                                                                      aTCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAAcc                                <  1:972137/150‑1 (MQ=255)

                                                                                                                        cACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCaa                              >  1:1157158/1‑150 (MQ=255)

                                                                                                                         aCCGGAACGGGCATTAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCacc                             >  2:244032/1‑148 (MQ=255)

                                                                                                                          ccGGAACTGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACg                            <  1:1549536/150‑1 (MQ=255)

                                                                                                                           cGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGc                           >  1:967106/1‑150 (MQ=255)

                                                                                                                                 gggCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATAc                     >  2:1403087/1‑150 (MQ=255)

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

                                                                                                                                  ggCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATAcc                    <  1:1588655/150‑1 (MQ=255)

                                                                                                                                    cAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGc                  <  2:25500/150‑1 (MQ=255)

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

                                                                                                                                         aaaTAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaata             >  1:1496582/1‑150 (MQ=255)

                                                                                                                                         aaaTAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGAACAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCTCCaata             >  2:1509993/1‑150 (MQ=255)

                                                                                                                                          tttAGTTTACTGCGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaataa            >  2:291327/3‑150 (MQ=255)

                                                                                                                                             aGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGc         <  1:1604031/150‑1 (MQ=255)

                                                                                                                                             aGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGc         <  2:1084853/150‑1 (MQ=255)

                                                                                                                                              gTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCg        >  1:670802/1‑150 (MQ=255)

                                                                                                                                               tttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGa       <  1:1306796/150‑1 (MQ=255)

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

                                                                                                                                                    gCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGataaat  <  1:1403087/150‑1 (MQ=255)

GCGACATCAGCCAGCGAGGACTTTCTGGCACGGTATACAGCAGCATTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAAT  >  NC_000913/4241567‑4241864

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