Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,870	A→C	63.2%	intergenic (‑953/‑215)	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,870	0	A	C	63.2%	83.7 / 59.9	77	intergenic (‑953/‑215)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base C (22/25);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGC  >  CP000730/1929735‑1930003                                                                                                                                            |                                                                                                                                          aTTCAATCTATCAGCTAGAATCGTATAACTGATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCtt                                                                                                                                     <  1:272708/141‑1 (MQ=255) aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTCGAGTCGCCTATCTCTCAGTCGTCAATTTAGACGCAGAATGGAGATGTA‑TAA‑G‑G‑TGC‑TGCtt                                                                                                                                     >  2:205052/1‑141 (MQ=255)   tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAt                                                                                                                                   <  1:151651/141‑1 (MQ=255)     aaTCCATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAAGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAt‑‑tt                                                                                                                               <  1:206777/141‑1 (MQ=255)       tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTc                                                                                                                             <  2:257350/141‑1 (MQ=255)              gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTAAATTTAGACGAAGAGAGGAGGTGTA‑GNA‑N‑G‑NGA‑TGCTTAT‑‑TTTAGT‑Tcaca                                                                                                                     >  1:265429/1‑141 (MQ=255)                       atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCAc                                                                                                            <  1:289258/141‑1 (MQ=255)                       atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑G‑TGA‑TGCTTCT‑‑TTTCGT‑TCACATCTTAGCAc                                                                                                            >  1:65610/1‑141 (MQ=255)                          taACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAg                                                                                                         <  2:287930/141‑1 (MQ=255)                                aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGtcatca                                                                                                   >  2:148616/1‑141 (MQ=255)                                 tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAg                                                                                                  >  1:106868/1‑141 (MQ=255)                                 tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAg                                                                                                  >  1:150063/1‑141 (MQ=255)                                 tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAg                                                                                                  >  1:134483/1‑141 (MQ=255)                                    atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGc                                                                                              <  1:110454/127‑1 (MQ=255)                                     tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATATCACCAGTCATCAGTGGCt                                                                                             >  2:106091/14‑141 (MQ=255)                                       ggaagtagataggCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                           <  2:14923/129‑1 (MQ=255)                                        gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACTAGTCATCAGTGGCtgtg                                                                                          >  1:346576/11‑141 (MQ=255)                                        gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                          <  1:127911/131‑1 (MQ=255)                                         aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                         <  1:5700/132‑1 (MQ=255)                                         aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                         >  2:251848/10‑141 (MQ=255)                                          aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                         <  1:57279/141‑1 (MQ=255)                                          aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                         >  2:119275/1‑141 (MQ=255)                                            tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                      <  2:205543/135‑4 (MQ=255)                                              gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTTTCAGGTGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                    <  1:67338/137‑6 (MQ=255)                                               atatgCTTGTATATTCACAGATGGTATAATAATAGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                   <  2:184533/138‑7 (MQ=255)                                                 ttgCTTGTATTCTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                  <  2:245553/140‑1 (MQ=255)                                                 atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GCCGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                 >  1:291399/2‑141 (MQ=255)                                                 atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                 >  1:189774/2‑141 (MQ=255)                                                 atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                 >  1:163342/2‑141 (MQ=255)                                                 atgCTTGTATATTCACAGAGGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTTAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                 >  1:232664/2‑141 (MQ=255)                                                    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                              <  2:216324/141‑4 (MQ=255)                                                    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                              <  2:275754/141‑4 (MQ=255)                                                     ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                             <  1:251848/141‑5 (MQ=255)                                                     ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                             <  2:173931/141‑5 (MQ=255)                                                        tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                          >  1:261721/11‑137 (MQ=255)                                                        tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                          <  1:96084/141‑8 (MQ=255)                                                          tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑‑‑TGCTAATGGTTACTTATCACATCATAGCACCAGTCATCAGGGGCTGGGCTGTTGCGTATTATACtt                                                                        >  1:99374/9‑141 (MQ=255)                                                              cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                    >  2:197643/5‑141 (MQ=255)                                                               acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTgg                                                                   <  1:217272/138‑1 (MQ=255)                                                                ccGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCCGTTATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                   <  2:58816/141‑1 (MQ=255)                                                                cagatggtATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                  <  1:175042/139‑1 (MQ=255)                                                                cagagggtataATAATAGAGTCGCCTATCTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                  <  2:357515/134‑1 (MQ=255)                                                                 cGAGGATATAATAGGAGAGTCGCCGATGTCTCAGGCGGCAATGTAGGCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCAGAGCACCAGTCATCAGTGGCTGTGCCACTGCGTTTTTTTCTTATTGGc                                                                  <  1:270609/141‑1 (MQ=255)                                                                 agatggtATAATAAAAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                 >  2:243199/2‑141 (MQ=255)                                                                  gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                >  1:87064/1‑140 (MQ=255)                                                                   atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                               >  1:38217/5‑139 (MQ=255)                                                                    tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                              >  1:285321/4‑138 (MQ=255)                                                                    tggtagaataagagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTGCGT‑TCACAGCATAGCACCAGTCATCAGGGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                              <  2:305221/129‑1 (MQ=255)                                                                    tggcataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                              <  2:241046/137‑4 (MQ=255)                                                                     ggtataATAATAAAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCTTAgg                                                             >  2:209550/3‑140 (MQ=255)                                                                       tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                            >  1:336584/1‑141 (MQ=255)                                                                       tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                           <  2:173789/141‑2 (MQ=255)                                                                       taaaataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCGCATCATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTTTTTTTCTTCATGGCTAAGTa                                                            >  2:24364/4‑141 (MQ=255)                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                          <  2:56056/141‑3 (MQ=255)                                                                           ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTGT‑C‑C‑TGT‑CAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                       <  1:82463/141‑1 (MQ=255)                                                                           ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTGTTATTGGCTCCGTAGACg                                                        >  1:294474/1‑141 (MQ=255)                                                                            tagtagAGTCGCCTATCTCTCCGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                       <  2:282721/141‑1 (MQ=255)                                                                            tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTTCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                       <  2:268317/141‑1 (MQ=255)                                                                            taatagAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                      <  1:249006/138‑1 (MQ=255)                                                                            aaagagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                      <  1:37940/137‑1 (MQ=255)                                                                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                  >  1:159474/1‑138 (MQ=255)                                                                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTGTTATACTTATTGGCTTAGTGAACGCaacaa                                                  >  1:33425/1‑138 (MQ=255)                                                                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑CACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TTACTAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCaacaa                                                  >  1:231223/1‑138 (MQ=255)                                                                                agagTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                  <  1:63764/141‑4 (MQ=255)                                                                                 gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                 >  1:270385/1‑137 (MQ=255)                                                                                       ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                           >  1:172911/1‑131 (MQ=255)                                                                                       ccTATCTCTCAGGCGTCAA‑‑TAGACGAAGGGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCTT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                           >  1:232189/1‑131 (MQ=255)                                                                                         tATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                         <  1:222125/141‑13 (MQ=255)                                                                                         tATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                         <  2:94746/141‑13 (MQ=255)                                                                                          ttctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATTGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                         <  2:313064/140‑1 (MQ=255)                                                                                          aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                        >  1:32157/1‑128 (MQ=255)                                                                                          aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                        >  1:221566/1‑128 (MQ=255)                                                                                          aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                        <  1:346969/141‑14 (MQ=255)                                                                                             tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTCTTTTCTTATTCGCTAAGTGGACGCACTCCACAATAGGCGAc                                      >  2:230209/1‑141 (MQ=255)                                                                                                    cGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                               >  2:347322/1‑141 (MQ=255)                                                                                                        aaTTTAGACGCAGAAAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCCT‑TCACATCAAAAAAAAAAAAAAAAGTGGCTGGGACAGTGCGTTTTTTTCTTACTGGCTACGTAGAAAAAATACAAAGCCGGTTTCATATAgaaaaa                           <  2:155225/141‑6 (MQ=255)                                                                                                                                  taA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc  <  2:212254/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          ATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGC  >  CP000730/1929735‑1930003

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