BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,918	A→G	52.4%	intergenic (‑1001/‑167)	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,918	0	A	G	52.4%	27.6 / 206.3	160	intergenic (‑1001/‑167)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (37/39); new base G (48/36); total (85/75)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

Only 100 of 161 total aligned reads displayed.

TGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTA  >  CP000730/1929778‑1930057

tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                                                                                                                                    >  2:756786/1‑141 (MQ=255)

 gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                                                                                   <  1:418512/141‑1 (MQ=255)

 gACTGCGTTTGTATACTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTTGTCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                                                                                   <  1:1077/141‑1 (MQ=255)

    tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTGC‑GTTCCCATCATAGCACCAGTCATCTGTGGCTGTGCCAtttc                                                                                                                                                >  2:331535/1‑139 (MQ=255)

    tgtgCTTGTATATTTGCCGATGATATAATAGTAAAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                                                                                >  1:350511/1‑141 (MQ=255)

     atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                              <  2:661469/140‑1 (MQ=255)

      tgCTTGTATATTTGGCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt                                                                                                                                              <  2:286128/141‑1 (MQ=255)

       gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGTCATTGCGtt                                                                                                                                             <  2:690961/141‑1 (MQ=255)

       gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                                                                            <  1:704492/141‑3 (MQ=255)

        cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                                           <  1:416547/141‑4 (MQ=255)

          tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                                                                         <  2:685505/141‑6 (MQ=255)

           gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                                        >  2:255734/12‑138 (MQ=255)

           gTATATTCCCAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                                        <  1:414022/141‑4 (MQ=255)

             atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                                                                                      >  2:452646/10‑136 (MQ=255)

                ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                   >  2:482053/7‑141 (MQ=255)

                ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                   <  2:641712/135‑1 (MQ=255)

                 tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                  <  1:588051/136‑1 (MQ=255)

                  tGCCGATGATATAATAGAAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACGAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt                                                                                                                                  >  1:244006/1‑141 (MQ=255)

                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACTAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 >  1:46857/5‑141 (MQ=255)

                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACTAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 >  1:26139/5‑141 (MQ=255)

                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                 >  2:254977/5‑141 (MQ=255)

                  cacagatggtATAATAATAAAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGGGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTG‑GCTCACAGCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAttt                                                                                                                                 >  2:229468/5‑140 (MQ=255)

                    ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                >  1:666752/1‑141 (MQ=255)

                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                               >  1:493152/3‑141 (MQ=255)

                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                               >  1:692460/3‑141 (MQ=255)

                       atggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                            <  1:622819/141‑3 (MQ=255)

                         ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGcatagt                                                                                                                          >  1:586731/3‑136 (MQ=255)

                          gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttcgta                                                                                                                         >  2:210700/2‑136 (MQ=255)

                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGGGGAGGTGCATTCC‑T‑T‑TGCTAA‑TCATTTTCACATCGCCGGA‑AGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                        >  1:177596/1‑140 (MQ=255)

                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTggtctagtaa                                                                                                                        >  2:15738/1‑133 (MQ=255)

                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                        <  1:309208/141‑2 (MQ=255)

                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                        <  1:642279/141‑2 (MQ=255)

                             taataaTGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                      <  2:251347/141‑4 (MQ=255)

                             taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                      >  2:300211/1‑141 (MQ=255)

                              aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                     >  2:173057/1‑141 (MQ=255)

                               ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                    <  2:740733/141‑1 (MQ=255)

                               ataCTAGAGTCGCCCACTTCTCACGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                    <  2:597991/141‑1 (MQ=255)

                                  gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa‑‑a                                                                                                                <  1:103884/141‑1 (MQ=255)

                                  atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                               <  1:230485/140‑2 (MQ=255)

                                  atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                               <  2:371298/140‑2 (MQ=255)

                                  atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAA‑TAATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                               >  2:195353/2‑140 (MQ=255)

                                  atagAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTACACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                               <  1:492174/140‑2 (MQ=255)

                                    agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATGCTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                             <  1:98392/141‑4 (MQ=255)

                                    agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAAGCa‑‑acaa                                                                                                             >  2:283390/1‑138 (MQ=255)

                                    agagTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                             >  1:466623/1‑138 (MQ=255)

                                     gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATAc                                                                                                             <  2:708963/141‑1 (MQ=255)

                                     gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                            <  2:521906/141‑5 (MQ=255)

                                     gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                            >  1:719709/1‑137 (MQ=255)

                                      agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACa                                                                                                            >  1:132538/1‑141 (MQ=255)

                                      agTCGCCTATCTCTCAAGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAGTGCGTTTTTTTCTTAGTGGCCAAGTAGACGCA‑‑ATACa                                                                                                            >  1:193572/1‑141 (MQ=255)

                                       gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAC‑G‑C‑T‑CTGCGTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAACTAGACGCA‑‑ATACaa                                                                                                           <  1:402258/141‑1 (MQ=255)

                                       gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa                                                                                                           <  1:193340/141‑1 (MQ=255)

                                        gcGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa‑‑a                                                                                                        <  1:79128/140‑1 (MQ=255)

                                         cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                        >  2:149529/1‑133 (MQ=255)

                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                       >  2:596056/1‑132 (MQ=255)

                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                       <  2:725208/141‑10 (MQ=255)

                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                       >  2:499373/1‑132 (MQ=255)

                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                       <  2:409843/141‑10 (MQ=255)

                                           ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                     >  2:263200/1‑141 (MQ=255)

                                           ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                     >  1:639477/1‑141 (MQ=255)

                                           ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                      >  2:266377/1‑131 (MQ=255)

                                           ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                      >  2:184971/1‑131 (MQ=255)

                                            cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                     >  1:140733/1‑130 (MQ=255)

                                            cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑G‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                     >  1:32603/1‑130 (MQ=255)

                                             tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaagtaattaat                                                                                                    >  1:286993/1‑129 (MQ=255)

                                              ttctctcAGGCGTTAATTTTGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTTTTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTCCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGt                                                                                                  <  2:51561/140‑1 (MQ=255)

                                                   tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑at                                                                                            <  2:36661/141‑1 (MQ=255)

                                                    cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAccga                                                                                            >  1:322778/1‑138 (MQ=255)

                                                    cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ata                                                                                           >  1:259031/1‑141 (MQ=255)

                                                    cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ata                                                                                           <  1:46988/141‑1 (MQ=255)

                                                    cAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAC‑TCATTTTA‑GTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCTGTTCCGTATTATACTTATTGGCTTAGTAAACGCA‑‑ACACAT‑‑AAttaataggcgac                                                                                           >  2:111110/1‑130 (MQ=255)

                                                       gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑Ag                                                                                       >  1:689926/1‑141 (MQ=255)

                                                           cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                  <  1:668149/141‑1 (MQ=255)

                                                            aaTTTAGACGCAGATAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGTC‑GCa                                                                                 <  1:547667/141‑1 (MQ=255)

                                                               ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                              <  1:123006/141‑1 (MQ=255)

                                                               ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                              <  1:742979/141‑1 (MQ=255)

                                                                tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gccacaa                                                                             >  1:3115/1‑136 (MQ=255)

                                                                 aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑GTAT‑AGCC‑GCACCAAt                                                                            <  1:568180/141‑1 (MQ=255)

                                                                 aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                            >  1:292570/1‑141 (MQ=255)

                                                                       agagagGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACAGCATAGCAACAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAAt                                                                      >  1:369635/1‑141 (MQ=255)

                                                                            ggaggTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                 >  1:653597/1‑141 (MQ=255)

                                                                               ggTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                              >  1:511171/1‑141 (MQ=255)

                                                                               ggTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                              >  2:321042/1‑141 (MQ=255)

                                                                                       aa‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                       <  2:676258/141‑1 (MQ=255)

                                                                                       aa‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                       <  1:723494/141‑1 (MQ=255)

                                                                                       aa‑G‑G‑GGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                       <  1:233930/141‑1 (MQ=255)

                                                                                                      tATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                            <  1:756786/141‑1 (MQ=255)

                                                                                                       aTTTTC‑GTCCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                           <  1:165511/141‑1 (MQ=255)

                                                                                                        ttttC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                          <  1:273750/141‑1 (MQ=255)

                                                                                                              gTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                     >  1:90260/1‑141 (MQ=255)

                                                                                                               ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTCCATTGCCAAGTCGCCGGG‑‑CCTCAA‑‑AAGAGGGGAC‑ATAT‑AGCA‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                    >  1:754183/1‑141 (MQ=255)

                                                                                                                         aGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaaa                          >  2:692661/1‑141 (MQ=255)

                                                                                                                          gCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcattaa                              >  2:23121/1‑138 (MQ=255)

                                                                                                                           cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                        >  2:432630/1‑141 (MQ=255)

                                                                                                                            aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                       >  1:699499/1‑141 (MQ=255)

                                                                                                                             ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACt                      >  1:77343/1‑141 (MQ=255)

                                                                                                                                  catcaGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                 <  1:727494/141‑1 (MQ=255)

                                                                                                                                         ggCTGTGCCATTGCGTTTTTTTCTTATTGGCGAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGt          <  2:578995/141‑1 (MQ=255)

                                                                                                                                               gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt    <  2:653597/141‑1 (MQ=255)

                                                                                                                                                ccATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AGACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTAGGGGGATTGGTGTATAAGTAAAGACTTATTATAGTTGTCTTAAtt   >  2:498752/1‑141 (MQ=255)

TGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTA  >  CP000730/1929778‑1930057

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