BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,866:1	+T	2.9%	intergenic (‑949/‑219)	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,866	1	.	T	2.9%	773.9 / 14.7	172	intergenic (‑949/‑219)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (93/74); new base T (2/3); total (95/77)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

Only 100 of 176 total aligned reads displayed.

GTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAA  >  CP000730/1929726‑1930005

       cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TAAATGCAGACGCAGAG‑AGGAGGGG‑TA‑CAA‑G‑G‑TGATGc                                                                                                                                      >  1:478470/1‑141 (MQ=255)

                 tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGt                                                                                                                            >  1:32723/1‑141 (MQ=255)

                 tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGt                                                                                                                            >  2:624848/1‑141 (MQ=255)

                         aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCCcatca                                                                                                                    >  2:455407/1‑141 (MQ=255)

                                 ataACTTATCAGAAGCAAGAGACTGTGCTTGTATATTTGCCGATGATATAATAGCAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                            >  1:503233/1‑141 (MQ=255)

                                    aCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCAATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTCTGTTCGTTCACGTCATAGCACCAGtc                                                                                                         >  1:377276/1‑141 (MQ=255)

                                            atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                <  1:117662/127‑1 (MQ=255)

                                              cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCGATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                              >  1:456350/13‑141 (MQ=255)

                                                gaagttggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑CGGAGTTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCtgtg                                                                                            <  1:255424/131‑1 (MQ=255)

                                                gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                            >  2:606190/11‑141 (MQ=255)

                                                 aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                           >  1:29408/10‑141 (MQ=255)

                                                 aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCCCCAGCCATCAGTGGCtgtg                                                                                            >  2:61050/1‑141 (MQ=255)

                                                  agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                          >  1:613140/9‑140 (MQ=255)

                                                  actagatacgCTTGTATATTCACAGACGGTATAATAATAGAGTTGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                          <  1:192797/132‑2 (MQ=255)

                                                   gtagatatgCCTGTATACTCACAGATGGTACAATAATAGCGTCGCCCACCTCGCGTCCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                         <  1:21780/134‑3 (MQ=255)

                                                    tagatatgCTTGTATATTCACAGATGGTATAAGAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                        <  2:464364/135‑4 (MQ=255)

                                                    atgatatgattGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                        <  2:156768/132‑4 (MQ=255)

                                                      gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                      >  1:24742/5‑136 (MQ=255)

                                                      aatattcTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                      <  1:706627/135‑6 (MQ=255)

                                                       atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgg                                                                                     >  1:534432/4‑135 (MQ=255)

                                                       atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                     <  1:315760/138‑7 (MQ=255)

                                                       atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                     <  1:684803/138‑7 (MQ=255)

                                                        tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                     <  1:60674/141‑1 (MQ=255)

                                                        tatgCTTGTATACTCACAGAAGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                    <  1:567610/139‑1 (MQ=255)

                                                         gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGGGGCG‑GA‑CAA‑G‑C‑AGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGCg                                                                                    >  1:177912/1‑141 (MQ=255)

                                                         gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCT‑AAAGAAGCAAAAGTGAGCCG‑AGTGAACGTCGAG‑TA‑TAA‑G‑G‑TGATGCTTGTTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTACCATTGCg                                                                                    >  1:370522/1‑141 (MQ=255)

                                                         atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                   <  1:105831/140‑1 (MQ=255)

                                                            cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                >  2:642691/1‑141 (MQ=255)

                                                              gttattTTGGCCGATGATATAATAGTAGAGTCGCTTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                               <  2:395127/139‑1 (MQ=255)

                                                               gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                             <  2:558452/141‑7 (MQ=255)

                                                                   attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                         >  2:170795/8‑140 (MQ=255)

                                                                    tttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TCAATTTAGACGCAGAG‑ACACGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGACTGTTCCCTTGCGtttttttactt                                                                         >  2:174548/1‑137 (MQ=255)

                                                                     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                       >  1:734401/6‑141 (MQ=255)

                                                                     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                       >  1:736519/6‑141 (MQ=255)

                                                                       acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                     >  1:687411/4‑141 (MQ=255)

                                                                       acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTgg                                                                     <  2:585235/138‑1 (MQ=255)

                                                                       acagaaggtataATAATAGAGTCGCCTAACTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                     >  1:390471/9‑141 (MQ=255)

                                                                        cagatggtATAATAATAGAGAAGCATAACTCTCTTACG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                    >  1:767767/3‑141 (MQ=255)

                                                                         agttggtataATAATAGAGTCGCCTAGCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                   <  1:639072/135‑1 (MQ=255)

                                                                         agatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                   >  2:211307/2‑141 (MQ=255)

                                                                          gatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                  <  1:165398/141‑2 (MQ=255)

                                                                          gatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                  >  2:489279/1‑140 (MQ=255)

                                                                            tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CC‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                >  1:297551/4‑141 (MQ=255)

                                                                            tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                >  1:324897/4‑141 (MQ=255)

                                                                            tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                <  2:425958/138‑4 (MQ=255)

                                                                            tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                >  1:430400/4‑141 (MQ=255)

                                                                            tggtataATAATAGAGGCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGCAGAG‑AGGAGGAG‑CATTAC‑T‑T‑TGCTAATCATTTTAGTTCACATAATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                >  1:646422/4‑141 (MQ=255)

                                                                               tataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                             <  1:538867/141‑2 (MQ=255)

                                                                               tataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                             >  2:668643/1‑140 (MQ=255)

                                                                               tataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                             >  2:579148/1‑140 (MQ=255)

                                                                                   ataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGc                                                         <  1:525159/141‑1 (MQ=255)

                                                                                   ataGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                          <  1:247138/141‑1 (MQ=255)

                                                                                   agtatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                         <  1:632896/137‑1 (MQ=255)

                                                                                    taatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                        <  1:416749/138‑1 (MQ=255)

                                                                                     aatagATTCGCCTATCTCTCAGGCG‑TCAA‑‑TGGACTAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                       <  2:770001/139‑1 (MQ=255)

                                                                                     aatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                       <  2:749139/139‑1 (MQ=255)

                                                                                      atagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                      >  2:565761/2‑140 (MQ=255)

                                                                                       tagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaaa                                                     >  1:31180/1‑139 (MQ=255)

                                                                                       tagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                     >  2:758991/1‑139 (MQ=255)

                                                                                       atgagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                     <  1:544764/139‑3 (MQ=255)

                                                                                        agagTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                     <  2:478470/141‑1 (MQ=255)

                                                                                        agagTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                     <  2:668581/141‑1 (MQ=255)

                                                                                        agagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                    <  1:244552/141‑4 (MQ=255)

                                                                                         gagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                   >  1:59597/1‑137 (MQ=255)

                                                                                         gagTCGCCTATCCCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                    <  2:41598/141‑1 (MQ=255)

                                                                                          agTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                  >  2:72282/1‑136 (MQ=255)

                                                                                           gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                  <  2:698408/141‑1 (MQ=255)

                                                                                           gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                  <  1:492273/141‑1 (MQ=255)

                                                                                           gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                  <  1:72687/141‑1 (MQ=255)

                                                                                           gTCGCCGATCTCTCAGGCG‑TCAATT‑‑GAAGAAGAG‑AGGAGGTGCGT‑TAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                 >  2:196975/1‑135 (MQ=255)

                                                                                            tCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                >  2:113229/1‑134 (MQ=255)

                                                                                              gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                              >  1:287339/1‑132 (MQ=255)

                                                                                              gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                              <  2:154858/141‑10 (MQ=255)

                                                                                              gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                              <  2:231325/141‑10 (MQ=255)

                                                                                              gCCTATCTCTAAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                              <  2:284388/141‑10 (MQ=255)

                                                                                               ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                             >  2:514869/1‑131 (MQ=255)

                                                                                               ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                             >  1:457055/1‑131 (MQ=255)

                                                                                               ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                             >  2:616714/1‑131 (MQ=255)

                                                                                                cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  1:320741/1‑130 (MQ=255)

                                                                                                cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  2:101946/1‑130 (MQ=255)

                                                                                                cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  1:744668/1‑130 (MQ=255)

                                                                                                cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  2:85138/1‑130 (MQ=255)

                                                                                                 tATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                           >  1:694320/1‑129 (MQ=255)

                                                                                                  aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                          >  2:277348/1‑128 (MQ=255)

                                                                                                   tctctcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAATACAAAATAGGTg                                          >  2:205758/1‑141 (MQ=255)

                                                                                                       tcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGAG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAACATGGGTGGCat                                      >  2:256390/1‑141 (MQ=255)

                                                                                                       ccAGGCG‑GCAATTTAGACGCAGAG‑AGGAGGGG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAGTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                      <  1:342103/140‑1 (MQ=255)

                                                                                                               tCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATCGCCg                               >  2:644272/1‑141 (MQ=255)

                                                                                                                 aaTTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATTCAAAATAGGTGAGATATAGCCGCa                             >  2:341349/1‑141 (MQ=255)

                                                                                                                        aCGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTATTCGGCCTGCCTAAATCGGCGCAATACAAAATAGGTGACATATAGCCGCCCCAAGaa                      >  2:316106/1‑141 (MQ=255)

                                                                                                                          gctgag‑cgGAGGTG‑TA‑TAA‑G‑G‑GGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                    <  2:379884/138‑1 (MQ=255)

                                                                                                                               g‑aggagGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc               <  1:624848/141‑1 (MQ=255)

                                                                                                                                 aggaggTG‑TG‑GAA‑G‑G‑TGATGCTTAGTTTCGTTCACATCATAGCACCCGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATCCAAAATAGGTGACATATAGCCGCACCAATAAAAATcctc              >  1:566864/1‑139 (MQ=255)

                                                                                                                                 aggaggTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAGTAAACATcccc              >  1:588547/1‑141 (MQ=255)

                                                                                                                                   gaggTG‑TA‑TAA‑G‑G‑TGATGATAGGGTGCTGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc            >  2:748512/1‑141 (MQ=255)

                                                                                                                                      gtgtta‑taA‑G‑G‑TGATGCTTATTTTCTGTCACATCATAGCACCATTCATCAGTGGCTGCGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGTACCAATAAAAATCCCCTAAc          <  2:299320/137‑1 (MQ=255)

                                                                                                                                       tg‑tA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa        <  2:246895/141‑1 (MQ=255)

                                                                                                                                            ttaA‑G‑G‑TCGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg     <  2:539066/140‑1 (MQ=255)

                                                                                                                                              aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa   >  1:193055/1‑141 (MQ=255)

                                                                                                                                               a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa  >  1:489336/1‑141 (MQ=255)

GTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAA  >  CP000730/1929726‑1930005

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