BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,799	C→A	44.0%	intergenic (‑882/‑286)	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,799	0	C	A	44.0%	84.4 / 117.1	122	intergenic (‑882/‑286)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (33/35); new base A (23/31); total (56/66)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 123 total aligned reads displayed.

GTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGC  >  CP000730/1929660‑1929939

gTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCg                                                                                                                                                  <  1:43883/141‑1 (MQ=255)

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

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

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

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

              gCCTACTTTCTTAAATTAACAATATCTATTCTCCTTACACTTGTCTGGTTGACTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                    <  1:295617/141‑1 (MQ=38)

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

              gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                    <  2:461480/141‑1 (MQ=255)

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

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

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

                      tCTTAAATTAACAATATCTATTCTCATAGAATGTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc                                                                                                                            <  2:386968/141‑1 (MQ=255)

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

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                           >  1:714853/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)

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

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

                                                 aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATCTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                 <  1:837743/141‑1 (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)

                                                        gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga                                                                                          >  2:180522/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)

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

                                                                   tgtAGACGATTCAAGCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAGCTCTCCGGCGTCAATTTAGACGCAGAGGGGAGGTG‑TA‑TAA‑g‑g                                                                           >  1:504572/1‑141 (MQ=255)

                                                                         cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGc                                                                    <  1:79055/141‑1 (MQ=255)

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

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

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

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

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

                                                                               aaTCTATCAGCTAGAATCATATAACTTATCAGCAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTACCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAttt                                                              <  1:607812/141‑1 (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)

                                                                                                              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                             >  1:259999/14‑141 (MQ=255)

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

                                                                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCGT‑TAC‑T‑T‑TCCTCATCGCTTTCGTTCACATCATAGCACCAGTCATCAGCGGCt                             >  1:616651/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)

                                                                                                                aGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                             >  1:559205/1‑141 (MQ=255)

                                                                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                           <  1:343615/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)

                                                                                                                      tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                      <  2:665036/135‑4 (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)

                                                                                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                     <  2:245269/136‑5 (MQ=255)

                                                                                                                        gatatgCTTGTATTTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                    >  1:572585/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)

                                                                                                                          tatgCTTGTATATTCACATATTGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                  <  2:318809/139‑1 (MQ=255)

                                                                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                  >  1:208309/3‑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)

                                                                                                                          gatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTGTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTCAGCATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                  <  2:56404/139‑1 (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)

                                                                                                                              cTTGAATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt              <  1:698023/141‑4 (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)

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

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

                                                                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac          <  2:307152/141‑5 (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‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGC  >  CP000730/1929660‑1929939

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: `‑`