BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,797	G→A	43.4%	intergenic (‑880/‑288)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,797	0	G	A	43.4%	69.2 / 128.3	121	intergenic (‑880/‑288)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (32/36); new base A (22/30); total (55/66)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTC  >  CP000730/1929660‑1929931

 tCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCAga                                                                                                                                         >  2:688181/1‑141 (MQ=255)

   ataataAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGTGAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgatg                                                                                                                                       <  2:829607/141‑1 (MQ=255)

          aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                >  1:208213/1‑141 (MQ=255)

            tCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtag                                                                                                                              >  2:582428/1‑141 (MQ=255)

            tCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtag                                                                                                                              >  1:371496/1‑141 (MQ=255)

            tCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTAACTGTGCTTGTATATTTGCCGATGATATAAaag                                                                                                                              >  1:371350/1‑141 (MQ=255)

              gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAGGCAAGTGACTGTGCTTGTATATTTGCCGATGAATTAAtagtg                                                                                                                            >  2:301731/1‑140 (MQ=255)

              gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                            >  1:820416/1‑141 (MQ=255)

              gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                            >  2:504010/1‑141 (MQ=255)

               ccTAGTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTGTCAGAAGCAAGTGACTGTGCCTGTATATTTGCCGATAAGAGAAtagtag                                                                                                                           >  2:551202/1‑141 (MQ=25)

                  aCTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGt                                                                                                                        <  1:270173/141‑1 (MQ=255)

                  aCTTTCTTAAACTAACAATATCTATTCTCATATAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTAAGCTTAGATATTTGCCGATGATATAATAGTAGAGt                                                                                                                        <  2:377837/141‑1 (MQ=37)

                      tCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc                                                                                                                    <  2:255910/141‑1 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                   >  1:714853/1‑141 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                   >  1:714883/1‑141 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAACTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCAACTGATATAATAGTAGAGTCGCCt                                                                                                                   <  2:457781/141‑1 (MQ=255)

                           aaTTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTAGCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtct                                                                                                               >  1:555796/1‑141 (MQ=255)

                               aaCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTCGAGTCGCCTATCTCTCa                                                                                                           <  2:603149/141‑1 (MQ=255)

                                  aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                        >  2:753867/1‑141 (MQ=255)

                                    tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGt                                                                                                      <  1:219668/141‑1 (MQ=255)

                                          ttCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAttt                                                                                                <  1:322815/141‑1 (MQ=255)

                                          ttCCCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAttt                                                                                                <  2:258342/141‑1 (MQ=255)

                                               aTAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACg                                                                                           <  1:726022/141‑1 (MQ=255)

                                                tAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGATTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGc                                                                                          >  2:239109/1‑141 (MQ=255)

                                                  gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                        <  1:57996/141‑1 (MQ=255)

                                                  gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                        <  1:502200/141‑1 (MQ=255)

                                                  gAATTTGTCCAATTAAGTGTAGACGATTCAATCGATCAGGTAGAATCATATAACTTATCAGCAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGGGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                        <  2:435133/141‑1 (MQ=255)

                                                   actttGTGCAATTACGCGTAGACGCGCCAATCCATCAGCTAGAATCATATAACTTACCCGAAACACGTGCCTGTGCTAGTCTATACCCCGATCATATAATAGTAGAGTCGCATACTGCTCCAGCGTCAGTTCAGCCGCaga                                                                                       >  2:34265/3‑141 (MQ=255)

                                                   aaattGTCCAATTAAGTGTAGACGATTAAATCTATCAGCTTGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGGCAATTTAGACGCaga                                                                                       <  2:785526/138‑1 (MQ=255)

                                                    aTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagag                                                                                      <  1:794918/141‑1 (MQ=255)

                                                      ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                    >  2:16187/1‑141 (MQ=255)

                                                          ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGgggagg                                                                                >  2:220951/1‑141 (MQ=255)

                                                          ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGaggagg                                                                                >  1:536532/1‑141 (MQ=255)

                                                          ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGACGCCTATCTCTCGGGCGTCAATGTAGACGCcgagagacgg                                                                                >  1:233739/1‑137 (MQ=255)

                                                               tAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑ta‑t                                                                         >  1:332680/1‑141 (MQ=255)

                                                                  gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑g                                                                     <  1:233854/141‑1 (MQ=255)

                                                                           aTTCAAGCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCtt                                                          <  1:19522/141‑1 (MQ=255)

                                                                             tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                        <  2:820416/141‑1 (MQ=255)

                                                                             tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                        <  2:857778/141‑1 (MQ=255)

                                                                             tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                        <  1:504010/141‑1 (MQ=255)

                                                                             tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                        <  1:301731/141‑1 (MQ=255)

                                                                              cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAtt                                                       >  1:828033/1‑141 (MQ=255)

                                                                                aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGTGGTG‑TA‑TAA‑G‑G‑TGATGCTTAtttt                                                     <  1:582428/141‑1 (MQ=255)

                                                                                aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGGTGCGTAtttt                                                     >  1:324975/1‑141 (MQ=255)

                                                                                   tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGt                                                  >  1:140705/1‑141 (MQ=255)

                                                                                            gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCActtcat                                         <  2:555796/141‑1 (MQ=255)

                                                                                             aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTGTATTTGCCGATGATATAATAGTAGAGTAGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGGTCACATCATa                                        >  1:196502/1‑141 (MQ=255)

                                                                                                   ataACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                  >  1:815784/1‑141 (MQ=255)

                                                                                                     aaCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGt                                <  2:283593/141‑1 (MQ=255)

                                                                                                           tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                          >  2:459763/1‑141 (MQ=255)

                                                                                                              atcggaattagatatgCATGTACATTCACAGATGCTATAATAATAGAGTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAACCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                      <  1:458780/127‑1 (MQ=255)

                                                                                                              atcggaagtatatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTT‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                      <  1:435092/127‑1 (MQ=255)

                                                                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATATAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                      <  1:115700/127‑1 (MQ=255)

                                                                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                      <  1:242005/127‑1 (MQ=255)

                                                                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                      <  2:616381/127‑1 (MQ=255)

                                                                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCATATCATAGCACCAGTCATCAGTGGc                      >  2:568844/15‑141 (MQ=255)

                                                                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGCTCAAGACATAGCACCAGTCATCAGCGGc                      >  2:560697/15‑141 (MQ=255)

                                                                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                     >  2:567492/14‑141 (MQ=255)

                                                                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                     >  1:112339/14‑141 (MQ=255)

                                                                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                     >  1:259999/14‑141 (MQ=255)

                                                                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCGT‑TAC‑T‑T‑TCCTCATCGCTTTCGTTCACATCATAGCACCAGTCATCAGCGGCt                     >  1:616651/14‑141 (MQ=255)

                                                                                                               tcggaagtagatatgCTTGTATATTCACAGATCGTATAATAATAGAGTCGCCTATCTCTCAGACGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                     >  1:653607/14‑141 (MQ=255)

                                                                                                               tcggaaggagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCt                     >  2:199155/14‑141 (MQ=255)

                                                                                                                gggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                    <  2:436160/129‑1 (MQ=255)

                                                                                                                aggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGACGCCTATCTCTCAGGCGTCAA‑‑TGGACGAAGGGAGGAGGTG‑CATTAC‑T‑C‑TGCTAAGCATCTCCGTTCACCTCATAGCCCCAGTGATCAGTGCCtg                    >  2:154804/13‑141 (MQ=255)

                                                                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                   <  1:606532/130‑1 (MQ=255)

                                                                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                   <  2:336666/130‑1 (MQ=255)

                                                                                                                   aagtagatatgCTTGTTTATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGc                 <  1:356595/132‑1 (MQ=255)

                                                                                                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                 <  2:638240/132‑1 (MQ=255)

                                                                                                                   aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                  >  1:189934/1‑141 (MQ=255)

                                                                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGGCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATGATAGCACCAGTCATCATTGGCTGTGctg               >  1:515352/8‑139 (MQ=255)

                                                                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTATGctg               >  1:725748/8‑139 (MQ=255)

                                                                                                                       gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATATTACGAGGCATCAGTGGCTGTGCCAt              <  2:646762/141‑1 (MQ=255)

                                                                                                                       gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt              <  1:691294/141‑1 (MQ=255)

                                                                                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt             >  1:559027/6‑137 (MQ=255)

                                                                                                                        gatatgCTTGTATTTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg            >  1:572585/5‑136 (MQ=255)

                                                                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg            >  1:556839/5‑136 (MQ=255)

                                                                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg            <  1:882442/137‑6 (MQ=255)

                                                                                                                          tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCATTGc           <  2:779727/141‑1 (MQ=255)

                                                                                                                          tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGc           >  1:122779/1‑141 (MQ=255)

                                                                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCg          >  1:208487/3‑141 (MQ=255)

                                                                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACAAGTCATCATTGGCTGTGCTGTTGCg          >  2:226741/3‑141 (MQ=255)

                                                                                                                           atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt         <  1:637536/140‑1 (MQ=255)

                                                                                                                           atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt         <  1:633230/140‑1 (MQ=255)

                                                                                                                            tgCTTGTATATTCAAAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑GATTAC‑T‑T‑TGCGAATCATTTTAGTTCACATCATAGCACCAGTCATCAGTTGCTGtgcttttgcttt        >  2:378795/1‑141 (MQ=255)

                                                                                                                              cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt      <  1:110514/141‑4 (MQ=255)

                                                                                                                              cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt      <  1:383506/141‑4 (MQ=255)

                                                                                                                              cTTGTATATTCACAGATGGTATAATAATAGAGGCGCCTATCTTTGAGGCGTCAAGTT‑‑ACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCTGTCATCAATGGCTGTGCTGTTGCGtatt      >  1:447839/1‑141 (MQ=255)

                                                                                                                               ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtccga     <  2:8910/141‑5 (MQ=255)

                                                                                                                               ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta     >  1:821819/1‑140 (MQ=255)

                                                                                                                               ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta     <  1:251866/141‑5 (MQ=255)

                                                                                                                               ttGTATATTCACAGAAGGTATAATATGAGAGTCGCCTAGCTCTCAGGCTTCAAGTT‑‑ACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGGCAGCAGTGGCTGTGCTGTTGCGtagta     <  2:166106/141‑5 (MQ=255)

                                                                                                                                tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGttttat    <  2:631480/141‑3 (MQ=255)

                                                                                                                                tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat    <  2:138253/141‑6 (MQ=255)

                                                                                                                                  tataTTTGCCGAGGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt   <  2:718397/141‑1 (MQ=255)

                                                                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac  <  2:307152/141‑5 (MQ=255)

                                                                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac  <  1:194973/141‑8 (MQ=255)

                                                                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac  <  2:429289/141‑8 (MQ=255)

                                                                                                                                   atatTTGCCGATGATATAATAGTAGAGTGGCATATCTCTCAGGCGTCAAGTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGTTTTTTTc  >  1:264965/1‑141 (MQ=255)

                                                                                                                                   atagtagCCTATTATATAATAGTGGAGTCGCCGATGTCTCAGGCGTCAATTTAGACGCGGAGAAGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTGTTTAAGTTCACCAAAAACAGCCAGTCATCCGTGGTTGTGCGATTGCGATttcttc  >  2:639315/7‑141 (MQ=255)

GTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTC  >  CP000730/1929660‑1929931

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