BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,875	T→C	66.2%	intergenic (‑958/‑210)	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,875	0	T	C	66.2%	140.9 / 100.5	125	intergenic (‑958/‑210)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (20/22); new base C (45/37); total (66/59)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 128 total aligned reads displayed.

TCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGG  >  CP000730/1929737‑1930015

  aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTCttt                                                                                                                                          >  1:55246/1‑141 (MQ=255)

                        ttgtgataaGAAGCCAGTGGCTGTGCTGGTAGCTTTGCCGATGATATTAGAGTCGAGTCGTTTCTCTCTCAGGCGTTTATTTAGACGCAGCGAGAAG‑GGG‑TA‑CAC‑G‑T‑TGATGCTTATTTTCGTTCACATCATGGCACCAGt                                                                                                                    <  2:137998/133‑1 (MQ=255)

                                aGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTg                                                                                                            >  2:338683/1‑141 (MQ=255)

                                 atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          <  1:20104/127‑1 (MQ=255)

                                 atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          >  1:313169/15‑141 (MQ=255)

                                 atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATAATAGCACCAGTCATCGGTGGc                                                                                                          >  1:22756/15‑141 (MQ=255)

                                 atcggaagtagatatgCTTGTATAATCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          >  2:271271/15‑141 (MQ=255)

                                  aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          <  1:41315/141‑1 (MQ=255)

                                   cggaagtagatatgCTTGTGTATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        >  1:248995/13‑141 (MQ=255)

                                   cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        >  2:362206/13‑141 (MQ=255)

                                   cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CAGTAC‑T‑T‑TGCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        <  2:191221/129‑1 (MQ=255)

                                   cggaagtagatatgCTTGGATATTCACCGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        <  1:125923/129‑1 (MQ=255)

                                   cggaaatagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        <  1:230718/129‑1 (MQ=255)

                                    ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGAGTC‑AT‑TAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                       <  1:106383/130‑1 (MQ=255)

                                    ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                       <  1:309083/130‑1 (MQ=255)

                                    ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                       <  1:121617/130‑1 (MQ=255)

                                     gaagtagatatgCTTTGATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                      <  2:39615/131‑1 (MQ=255)

                                     gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                      <  1:269895/131‑1 (MQ=255)

                                      aagtaggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                     <  2:77730/132‑1 (MQ=255)

                                       agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                    <  2:284964/133‑2 (MQ=255)

                                        gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                   <  1:286635/134‑3 (MQ=255)

                                        gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                   <  1:245413/134‑3 (MQ=255)

                                        gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTAGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                   >  1:341078/8‑139 (MQ=255)

                                         tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                                  <  1:350612/135‑4 (MQ=255)

                                         tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGctgt                                                                                                  <  2:300355/135‑4 (MQ=255)

                                         tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                                                                                   <  2:332439/141‑1 (MQ=255)

                                           gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                                <  2:29063/137‑6 (MQ=255)

                                            atattcTTGTATATTCACAGAGGGTATAATAATAGATTCGCTTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGc                                                                                               <  2:188222/136‑1 (MQ=255)

                                              atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                             >  2:47848/2‑141 (MQ=255)

                                              atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                             >  2:60076/2‑141 (MQ=255)

                                               tgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt                                                                                             >  1:166975/1‑141 (MQ=255)

                                                 cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGGCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCCACAGTCACCAGTGGCTATGCCCTTGCGatt                                                                                           >  1:3902/1‑141 (MQ=255)

                                                  ttGTATATTTGCCGATGATATAATAGTATAGTTGCCTATCGCTTAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                                          <  1:107095/141‑1 (MQ=255)

                                                    gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                       <  1:177603/141‑7 (MQ=255)

                                                         ttcacagagggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑A‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                  >  2:264304/12‑141 (MQ=255)

                                                           tGCCGATGATATAATAGTAGAGTCTCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt                                                                                 <  1:336731/141‑1 (MQ=255)

                                                           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                >  1:94314/5‑141 (MQ=255)

                                                            acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                               <  2:293779/138‑1 (MQ=255)

                                                             cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                              >  2:265441/3‑141 (MQ=255)

                                                              agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                             >  2:57144/2‑141 (MQ=255)

                                                              agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                             >  1:288208/2‑141 (MQ=255)

                                                               gatggtATAGTAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GGG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                            <  2:228541/141‑2 (MQ=255)

                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                          >  1:378026/4‑141 (MQ=255)

                                                                 tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                                                           >  2:362424/1‑141 (MQ=255)

                                                                 tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                                                           >  2:102514/1‑141 (MQ=255)

                                                                 gggtataATAATAGAGCCGCCTATCTCTCAGGCGCCAATT‑‑GACGAAGAGAGGAG‑GTGCGT‑TAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATGCTTATTGGCTTAg                                                                          >  1:230187/4‑141 (MQ=255)

                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                         >  1:264537/3‑141 (MQ=255)

                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGTG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                         <  1:374042/139‑1 (MQ=255)

                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATCCTTCTTGGCTTAGt                                                                         >  2:56730/3‑141 (MQ=255)

                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                         >  2:236422/3‑141 (MQ=255)

                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                         <  2:190065/139‑1 (MQ=255)

                                                                   gtataATAATAGAGACGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                        >  1:200061/2‑141 (MQ=255)

                                                                    tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                       >  2:320387/1‑140 (MQ=255)

                                                                    tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                       >  2:291491/1‑140 (MQ=255)

                                                                    tataATAATAGAGTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                       <  1:258578/141‑2 (MQ=255)

                                                                      taataaTAGAGTCGCCTATGTCTCAGGCGTCAATT‑‑GACGAAGATAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                     >  1:314965/1‑141 (MQ=255)

                                                                       aataatAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                    >  1:30782/1‑141 (MQ=255)

                                                                        atatTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGGG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCTGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttactaaacgc                                                                   <  2:269373/141‑11 (MQ=255)

                                                                        ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                   <  2:66599/137‑1 (MQ=255)

                                                                        agaagagagCCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑G‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                   <  1:267206/136‑1 (MQ=255)

                                                                         taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                  <  2:282803/138‑1 (MQ=255)

                                                                         taatagAGTCGCCTATCTCTCAGGCGTCAAAT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                  <  1:127303/138‑1 (MQ=255)

                                                                           atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTGTACTTATTGGCTTGGTAAACGCaac                                                                >  1:20020/2‑140 (MQ=255)

                                                                            tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                               <  2:15153/141‑3 (MQ=255)

                                                                             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                              >  2:64214/1‑138 (MQ=255)

                                                                              gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                             >  1:75619/1‑137 (MQ=255)

                                                                              gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                             <  1:289772/141‑5 (MQ=255)

                                                                               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaat                                                            >  2:353305/1‑136 (MQ=255)

                                                                               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                            <  1:52210/141‑6 (MQ=255)

                                                                                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                           >  2:143415/1‑135 (MQ=255)

                                                                                 tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                          >  2:215102/1‑134 (MQ=255)

                                                                                  cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                          >  1:153740/1‑141 (MQ=255)

                                                                                  cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                          <  1:103237/141‑1 (MQ=255)

                                                                                   gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                         >  1:182836/1‑141 (MQ=255)

                                                                                   gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                         <  1:258115/141‑1 (MQ=255)

                                                                                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                        <  2:10466/141‑10 (MQ=255)

                                                                                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                        >  2:359251/1‑132 (MQ=255)

                                                                                   gCCCATCTCTCAGGCGTCAATT‑‑GGCGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                        <  2:354581/141‑10 (MQ=255)

                                                                                     cTATCTCTCCAGCGTCAATTTAGAGGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGCTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAg                                                       >  1:280001/1‑141 (MQ=255)

                                                                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                      >  1:223191/1‑130 (MQ=255)

                                                                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                      <  1:241598/141‑12 (MQ=255)

                                                                                       aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                     <  2:382392/141‑1 (MQ=255)

                                                                                       aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGAGGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGACGCAAGACAAACTactt                                                     >  1:80999/1‑138 (MQ=255)

                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                    >  1:144825/1‑128 (MQ=255)

                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                    <  1:362206/141‑14 (MQ=255)

                                                                                        tctctcAGGCGTCAATT‑‑GACGAAGGGAGGAG‑GTGCGT‑TAC‑T‑T‑TGCTACTCATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATAGTTATTGCCTTAGTAAACGCaaccagtatttaattg                                                   >  1:124039/1‑127 (MQ=255)

                                                                                             cAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGTAATACAAAATAGGTGACata                                               >  1:282712/1‑141 (MQ=255)

                                                                                             cAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                               <  2:321808/141‑1 (MQ=255)

                                                                                                     aaTTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                                       <  1:277385/141‑1 (MQ=255)

                                                                                                         tAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCaa                                   <  1:68106/141‑1 (MQ=255)

                                                                                                              gCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGATGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                              <  1:141597/141‑1 (MQ=255)

                                                                                                                agagagGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAAt                            >  1:146014/1‑141 (MQ=255)

                                                                                                                        g‑gTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                    <  1:102514/141‑1 (MQ=255)

                                                                                                                        g‑gTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                    <  1:362424/141‑1 (MQ=255)

                                                                                                                              ta‑taA‑G‑G‑TGGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGACGCACCAATAACAATCTCCTCACTAcg                >  2:334798/1‑140 (MQ=255)

                                                                                                                                             gCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa     >  1:98768/1‑141 (MQ=255)

                                                                                                                                             gCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa     >  2:374656/1‑141 (MQ=255)

                                                                                                                                               ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg   <  2:280001/141‑1 (MQ=255)

                                                                                                                                                tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTACGTAGAAGCAATAAAAAGTAGGTGACATATAGCCGCACCAATGAAAAACCCCACACTACCGAAAATAgaaggcg  >  2:199963/1‑135 (MQ=255)

                                                                                                                                                tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg  <  2:333771/141‑1 (MQ=255)

TCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGG  >  CP000730/1929737‑1930015

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