BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,848	C→A	59.6%	intergenic (‑931/‑237)	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,848	0	C	A	59.6%	97.9 / 229.1	190	intergenic (‑931/‑237)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (34/42); new base A (59/53); total (94/96)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 190 total aligned reads displayed.

GAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAA  >  CP000730/1929710‑1929985

gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                                                                               <  1:579322/141‑1 (MQ=255)

 aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGCTATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                                                              <  2:491911/141‑1 (MQ=255)

  aTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAAACATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGAAATAAGAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaaag                                                                                                                                             >  2:67412/1‑141 (MQ=255)

       tCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgaggag                                                                                                                                        <  1:194253/141‑1 (MQ=255)

         cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATGATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATGTAATAGTAGAGTCGCCTATCTCTCAGGCGTGAATTTAGACGCAGAGaggaggg                                                                                                                                      >  1:415104/1‑140 (MQ=255)

         cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGCCGCTTCTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgt                                                                                                                                      <  1:750339/141‑1 (MQ=255)

         cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgt                                                                                                                                      <  1:24811/141‑1 (MQ=255)

                        gATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGAGGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑t                                                                                                                  >  2:286778/1‑141 (MQ=255)

                             aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTAttt                                                                                                             <  1:622391/141‑1 (MQ=255)

                                      gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTcaca                                                                                                   >  2:242692/1‑141 (MQ=255)

                                                 ataACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATACTCGCTAAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCa                                                                                        <  1:556458/141‑1 (MQ=255)

                                                       tATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGtcatc                                                                                  <  1:8121/141‑1 (MQ=255)

                                                            gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCCAATCATAGCACCAGTCATCAGTgg                                                                             >  2:660837/1‑141 (MQ=255)

                                                            gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTgg                                                                             >  1:36661/1‑141 (MQ=255)

                                                            atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGc                                                                            >  2:288542/15‑141 (MQ=255)

                                                            atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGCCATCAGCGGc                                                                            >  2:583412/15‑141 (MQ=255)

                                                            atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCACCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGc                                                                            >  1:508949/15‑141 (MQ=255)

                                                             tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAA‑TCATTTTC‑GTTCCCATCATAGCACCAGTCATCAGTGGcg                                                                           >  2:328926/14‑140 (MQ=255)

                                                             gcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCt                                                                           <  1:87935/128‑1 (MQ=255)

                                                             aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATATCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGtcagcctttga                                                                            >  2:19820/1‑133 (MQ=255)

                                                              cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                          >  2:400108/13‑141 (MQ=255)

                                                              aGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCt                                                                           >  2:568180/1‑141 (MQ=255)

                                                               ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                         <  1:62742/130‑1 (MQ=255)

                                                               ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                         <  1:63586/130‑1 (MQ=255)

                                                                gaagtagatatgCTTGTTTATTCACAGATGATATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGATGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                        <  2:725921/131‑1 (MQ=255)

                                                                 gagtcgatatgCTTGTATATTCACAGATGGTATCAGAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                       <  1:245784/132‑1 (MQ=255)

                                                                 aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                       >  1:434977/10‑141 (MQ=255)

                                                                  agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCCCCAGTCATCAGTGGCTGTGct                                                                      >  2:74521/9‑140 (MQ=255)

                                                                  aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                       >  2:233191/1‑141 (MQ=255)

                                                                    tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                    <  1:7058/135‑4 (MQ=255)

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

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

                                                                     gACCGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATGTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCCAt                                                                    <  2:702695/141‑1 (MQ=255)

                                                                      gatttgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                  <  1:48652/137‑6 (MQ=255)

                                                                      gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACAAGTCATCAGTGGCTGTGctgttg                                                                  >  1:11691/5‑136 (MQ=255)

                                                                      aatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                  <  1:727175/137‑6 (MQ=255)

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

                                                                        tggccTTGTATATTTGCCCATGATATAATAGTAGTGTCGCCTATCGCTCAGGCGTCAATTTAGACGTAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                 <  2:206841/137‑1 (MQ=255)

                                                                         atgCTTGTATATTCACAGATGGTATAATACTAGAGGCGCCTATCTCTCAGGCGTCAA‑‑TGGACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCCCCAGTCATCAGTGGCTGTGctgttacct                                                               >  2:658256/2‑133 (MQ=255)

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

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

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

                                                                              tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                          >  1:762306/1‑139 (MQ=255)

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

                                                                                tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                        <  1:384396/141‑8 (MQ=255)

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

                                                                                  aattcacagatggtAAAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                      <  2:357554/133‑1 (MQ=255)

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

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

                                                                                    ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCTTAt                                                    >  2:116941/7‑141 (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‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                  >  1:744835/5‑141 (MQ=255)

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

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

                                                                                        cagatggtATAATAATATAGTGGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                <  2:183178/139‑1 (MQ=255)

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

                                                                                         agatggtATAATAATCGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGTAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACt                                               >  2:214743/2‑141 (MQ=255)

                                                                                            aggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                            <  2:690362/138‑4 (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)

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

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

                                                                                                    taatagAGTCGCCTATCTCTCGGGCGGCAATT‑‑GACGAAGCGAGAAGGAATAGTTC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                    <  1:362513/138‑1 (MQ=255)

                                                                                                     aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                   <  1:702817/139‑1 (MQ=255)

                                                                                                     aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                   <  1:323271/139‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                >  1:13729/1‑138 (MQ=255)

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

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

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

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

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

                                                                                                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                             >  1:742323/1‑135 (MQ=255)

                                                                                                            tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                            >  1:653965/1‑134 (MQ=255)

                                                                                                             cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                           >  1:645293/1‑133 (MQ=255)

                                                                                                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                          <  1:317320/141‑10 (MQ=255)

                                                                                                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                          >  2:412835/1‑132 (MQ=255)

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

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

                                                                                                                cTATCTCTCAGGCGTCAATTTAGACGCCGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAg                         >  2:529813/1‑141 (MQ=255)

                                                                                                                cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTTAGTAGATGCAATACAAAATAg                         >  2:712681/1‑141 (MQ=255)

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

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

                                                                                                                    ctatctggggTCCATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGa                     <  2:269466/132‑1 (MQ=255)

                                                                                                                        cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                 >  1:259031/1‑141 (MQ=255)

                                                                                                                        cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                 >  2:495019/1‑141 (MQ=255)

                                                                                                                               cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc          <  1:668149/141‑1 (MQ=255)

                                                                                                                               cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATACGTGACATATAGCCGc          >  1:517510/1‑141 (MQ=255)

                                                                                                                                   ttAGACGCAGAGGGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa      <  1:146770/141‑1 (MQ=255)

                                                                                                                                     aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt    >  1:292570/1‑141 (MQ=255)

                                                                                                                                       aCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa  <  1:27185/141‑1 (MQ=255)

GAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAA  >  CP000730/1929710‑1929985

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