BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,873	C→A	65.8%	intergenic (‑956/‑212)	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,873	0	C	A	65.8%	112.2 / 117.9	125	intergenic (‑956/‑212)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (21/22); new base A (45/37); total (66/59)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 125 total aligned reads displayed.

TCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGA  >  CP000730/1929737‑1930012

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

    tCTAGCAGCTAGAATCTTATACCTTTTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTc                                                                                                                                     <  2:153346/141‑1 (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)

                                  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‑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)

                                     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)

                                            atattcTTGTATATTCACAGAGGGTATAATAATAGATTCGCTTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGc                                                                                            <  2:188222/136‑1 (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)

                                              atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                          >  2:28660/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)

                                                            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                                                                           >  1:284786/3‑141 (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‑GTGCCT‑TCC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                       >  1:53099/4‑138 (MQ=255)

                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAG‑GTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                       >  1:378026/4‑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                                                                      >  1:301275/3‑141 (MQ=255)

                                                                  gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                       <  2:310633/141‑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)

                                                                    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)

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

                                                                                     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)

                                                                                            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:218594/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)

                                                                                                                    aggag‑gTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc                     >  1:82811/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)

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

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

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

TCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAG‑GTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGA  >  CP000730/1929737‑1930012

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