BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,864	A→T	11.3%	intergenic (‑947/‑221)	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,864	0	A	T	11.3%	343.7 / 25.5	124	intergenic (‑947/‑221)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (55/54); new base T (8/6); total (64/60)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 125 total aligned reads displayed.

GTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAC  >  CP000730/1929726‑1930000

gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GGG‑TA‑TAA‑g                                                                                                                                          >  2:141597/1‑141 (MQ=255)

   tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑Tg                                                                                                                                     <  1:374085/141‑1 (MQ=255)

           tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                                                                                             >  1:324162/1‑141 (MQ=255)

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

               tCTAGCAGCTAGAATCTTATACCTTTTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTc                                                                                                                         <  2:153346/141‑1 (MQ=255)

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

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

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

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

                                             tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                          >  1:146370/14‑141 (MQ=255)

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

                                             aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                           >  2:333846/1‑141 (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:121617/130‑1 (MQ=255)

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

                                               gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                         >  1:382392/1‑141 (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                                                                                    >  2:353669/8‑139 (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)

                                                     agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTAACTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCACTTTCGGTCACATCATAGCACCAGTCATAAGTGGCTGTGctgtt                                                                                  >  1:269373/6‑137 (MQ=255)

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

                                                        tgtgCTTTTATTTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                <  2:26293/141‑1 (MQ=255)

                                                        tatgCTTGTTGATTCACAGATGGTATAAAAATAGAGTCGCCTATCTCTCAGGCGTTAATTT‑‑ACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                               <  2:6360/139‑1 (MQ=255)

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

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

                                                           gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                            <  2:275243/141‑3 (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)

                                                                tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCCATTGCGttttttt                                                                        <  1:232460/141‑1 (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‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                               >  2:47341/3‑141 (MQ=255)

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

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

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

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

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

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

                                                                            gggtataATAATAGAGCCGCCTATCTCTCAGGCGCCAATT‑‑GACGAAGAGAGGAG‑GTGCGT‑TAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATGCTTATTGGCTTAg                                                           >  1:230187/4‑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‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                          >  1:301275/3‑141 (MQ=255)

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

                                                                              gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                         >  2:115820/2‑141 (MQ=255)

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

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

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

                                                                                 taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                       <  1:252728/141‑1 (MQ=255)

                                                                                  aataatAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                     >  1:30782/1‑141 (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‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGAAAACGCaacaaa                                              >  2:214494/1‑137 (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:103237/141‑1 (MQ=255)

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

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

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

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

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

                                                                                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                       >  1:223191/1‑130 (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)

                                                                                                       ccAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCATGGATGAGTATACGCAATACAAAATAGGTGACat                                 <  1:185038/140‑1 (MQ=255)

                                                                                                       ccAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                 <  2:363207/140‑1 (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)

                                                                                                                aaTTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                        <  1:218594/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:362424/141‑1 (MQ=255)

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

GTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAC  >  CP000730/1929726‑1930000

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