BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,848	C→A	71.6%	intergenic (‑931/‑237)	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,848	0	C	A	71.6%	307.3 / 176.3	235	intergenic (‑931/‑237)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (32/36); new base A (88/79); total (120/115)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

Only 100 of 236 total aligned reads displayed.

AGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAA  >  CP000730/1929709‑1929986

  aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                                                               <  2:234262/141‑1 (MQ=255)

       gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga                                                                                                                                          >  2:677448/1‑141 (MQ=255)

       gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga                                                                                                                                          >  2:829688/1‑141 (MQ=255)

                   gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGGGGAGGGG‑TA‑TAA‑G‑G‑t                                                                                                                         >  1:491417/1‑141 (MQ=255)

                     aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATATCTCAGGCGTCAATTTAGTCGCAGTGAGGAGGTG‑TA‑TAA‑G‑G‑T‑Ga                                                                                                                      <  1:516640/141‑1 (MQ=255)

                                                         aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                                  >  1:500526/1‑141 (MQ=255)

                                                         aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                                  <  1:95442/141‑1 (MQ=255)

                                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                             <  1:385809/127‑1 (MQ=255)

                                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                             >  1:403493/15‑141 (MQ=255)

                                                               cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                           <  1:273132/129‑1 (MQ=255)

                                                               cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                           >  2:552933/13‑141 (MQ=255)

                                                                ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                          <  1:698084/130‑1 (MQ=255)

                                                                 gaagtagatatgCTTGTATATTCACAGGTGGTACAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                         <  1:229031/131‑1 (MQ=255)

                                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                         <  1:86919/131‑1 (MQ=255)

                                                                 gaagtagagtttcTTGGATATTCACAGTTGGTATAATCATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                         <  1:142484/129‑1 (MQ=255)

                                                                  aatgagagtagcgtGTATATTCACAGATGGTATAATAATAGGGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                        <  2:409890/128‑1 (MQ=255)

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

                                                                  aaggagagatgCTTGTATATTCACAGATTGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                        <  2:762431/132‑1 (MQ=255)

                                                                   aGTGACTGTGCTTGTATATTTGCCGATGATATAACAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                        <  2:406851/141‑1 (MQ=255)

                                                                   aGTGAATATGCTTGTATATTCATAGATGGTATAATAATAGTGTCGTCTATCTCTCAGGCGTCAATT‑‑GACGAAGTGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                       <  2:213878/141‑2 (MQ=255)

                                                                    gtagatatgattGTATCTTCACCGATGGTATAGGAATAGAGTCGCTTATCGCTCAGGGGGCAATT‑‑GACGAAGAGAGGAGGGG‑CA‑TTA‑C‑T‑TGGCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                      <  2:890004/131‑3 (MQ=255)

                                                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                      <  1:583158/134‑3 (MQ=255)

                                                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                      >  2:69433/8‑139 (MQ=255)

                                                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGACTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CAGTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTCGCcgttcgt                                                                      >  2:278690/8‑137 (MQ=255)

                                                                    gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                                                       >  2:588576/1‑141 (MQ=255)

                                                                     tGACTGTGCTTGTATATTTGCGGATGATATAATAGGAGAGTCGCCTATCTCGCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                                                      <  2:625450/141‑1 (MQ=255)

                                                                      agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                    >  1:745034/6‑137 (MQ=255)

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

                                                                       gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                   >  2:651107/5‑136 (MQ=255)

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

                                                                       gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                   >  1:128562/5‑136 (MQ=255)

                                                                       aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                    >  1:185008/1‑141 (MQ=255)

                                                                       aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                    <  2:529880/141‑1 (MQ=255)

                                                                       aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                    >  2:723604/1‑141 (MQ=255)

                                                                         tgtgCTTGTATATTTGCCGATGATATAATCGTAGCGTCGACTATCAAAAGGGAGTAAGATGATGCGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                  <  2:198155/141‑1 (MQ=255)

                                                                          atgCTTGTATATTCACAGATGGTATAATAGTAGAGTCACCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑G‑T‑T‑GCTAATCATTTTGGTTCGCATCATAGAAACCACCATCAGTGGCTGTGctggtgcgt                                                                >  2:704367/2‑133 (MQ=255)

                                                                          atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                >  2:936563/2‑141 (MQ=255)

                                                                          atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                >  1:567940/2‑141 (MQ=255)

                                                                            gCTTGTCTATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                              <  2:801468/141‑3 (MQ=255)

                                                                             cTTGGATATTCACAGATGGTAAAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                             <  1:803171/141‑1 (MQ=255)

                                                                                gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                          <  2:144288/141‑7 (MQ=255)

                                                                                 tataTTCACAGATGGTATAATAACAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                         <  1:391189/141‑8 (MQ=255)

                                                                                  atattcaaagatggtATAATAATAGAGTCGCCTCTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                        >  2:344770/10‑136 (MQ=255)

                                                                                     tttGCCGATGATATAATACTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTa                                                      <  2:922020/141‑1 (MQ=255)

                                                                                     ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAAGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                     >  1:604855/7‑141 (MQ=255)

                                                                                      tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                    >  2:428152/6‑141 (MQ=255)

                                                                                      tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                    >  2:193956/6‑141 (MQ=255)

                                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                 <  2:353140/139‑1 (MQ=255)

                                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                 >  1:269137/3‑141 (MQ=255)

                                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                 <  1:296647/139‑1 (MQ=255)

                                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGc                                                 >  1:466031/3‑141 (MQ=255)

                                                                                          cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                 <  1:777770/141‑1 (MQ=255)

                                                                                          cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTTGc                                                 >  2:872463/1‑141 (MQ=255)

                                                                                          agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATTGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                <  2:612111/140‑1 (MQ=255)

                                                                                          agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                >  1:839811/2‑141 (MQ=255)

                                                                                          agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                <  2:872901/140‑1 (MQ=255)

                                                                                           gatggtATAATAACAGGGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                               <  2:730648/141‑2 (MQ=255)

                                                                                           gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                <  1:544157/141‑1 (MQ=255)

                                                                                            atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                              >  1:794153/5‑139 (MQ=255)

                                                                                            atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                              >  1:302500/5‑139 (MQ=255)

                                                                                            atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                               <  2:69720/141‑1 (MQ=255)

                                                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGTGTATTATACTTATTGGCttag                                             >  1:750865/4‑138 (MQ=255)

                                                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                             <  1:531450/138‑4 (MQ=255)

                                                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                             <  1:106099/138‑4 (MQ=255)

                                                                                             tggtataATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                             <  2:29864/138‑4 (MQ=255)

                                                                                             tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTCCGTTTTTTTCTTATTGGCCaa                                              >  1:169448/1‑141 (MQ=255)

                                                                                              ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                            >  1:612586/3‑141 (MQ=255)

                                                                                              ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                            >  1:204431/3‑141 (MQ=255)

                                                                                              gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                             >  2:494724/1‑141 (MQ=255)

                                                                                               gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                           >  1:139950/2‑141 (MQ=255)

                                                                                                tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTC‑CA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGAGCTGTTGCGTATTATACTTATTGGCtgaataa                                          >  1:632710/1‑135 (MQ=255)

                                                                                                 ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                         >  1:4748/1‑139 (MQ=255)

                                                                                                  taataagagcggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                        <  1:42065/129‑1 (MQ=255)

                                                                                                        tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTACTTAGACGCAAt                                   >  2:259973/1‑141 (MQ=255)

                                                                                                        tagAGCCGTCTATTTCTCAGTCGGCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTTATCAGCGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                  <  2:852145/141‑3 (MQ=255)

                                                                                                         agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                  <  1:913395/141‑1 (MQ=255)

                                                                                                         agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGTACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                 <  2:932192/141‑4 (MQ=255)

                                                                                                         agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                 >  1:864551/1‑138 (MQ=255)

                                                                                                         agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                 >  2:682593/1‑138 (MQ=255)

                                                                                                          gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaagaat                                >  1:442004/1‑137 (MQ=255)

                                                                                                          gagTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                <  1:733450/141‑5 (MQ=255)

                                                                                                           agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                >  2:60435/1‑141 (MQ=255)

                                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                               >  1:704799/1‑136 (MQ=255)

                                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCACTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                               <  1:625419/141‑6 (MQ=255)

                                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                               <  1:687557/141‑6 (MQ=255)

                                                                                                              cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                            >  2:272426/1‑133 (MQ=255)

                                                                                                              cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                            >  2:139467/1‑133 (MQ=255)

                                                                                                                ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                           >  2:891788/1‑141 (MQ=255)

                                                                                                                ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTC‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGCTGTGCTGTTGCGTATTATGCATATTGGGTTAGTAAACGCAACAaaaactga                          >  1:82009/1‑137 (MQ=255)

                                                                                                                 cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaccaaataattaa                         >  1:354737/1‑129 (MQ=255)

                                                                                                                 cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                         >  1:769679/1‑130 (MQ=255)

                                                                                                                 cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                         >  1:322458/1‑130 (MQ=255)

                                                                                                                  tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCAA‑TAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                        >  1:526362/1‑129 (MQ=255)

                                                                                                                   aTCTCTCAGGCGTCACTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                        >  1:937040/1‑141 (MQ=255)

                                                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGTTTAGTAAACGCaacaaataattgata                       >  1:714068/1‑128 (MQ=255)

                                                                                                                    tctctcATGGGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                      <  2:365915/141‑15 (MQ=255)

                                                                                                                    tctctcAGGCGTCAATTTAGACGCAGAGAGGCGGCG‑TA‑CAA‑G‑G‑T‑GATGCTTATTTTAGTTCACAACATAGAACCAGTCATCAGTGGATGTGCCATTGCGTTTTTTGCTTATTGGCTAAGTTGGCGCAATCCAAACTAGGTg                       >  2:453324/1‑141 (MQ=255)

                                                                                                                          atgCGTAAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGGTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                 <  1:550462/139‑1 (MQ=255)

                                                                                                                            gCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTTGATGCAGTACAAAATAGGTGACATATAg               >  1:679796/1‑141 (MQ=255)

                                                                                                                                 aaTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa          <  1:101620/141‑1 (MQ=255)

AGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAA  >  CP000730/1929709‑1929986

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