BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,872	G→A	65.4%	intergenic (‑955/‑213)	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,872	0	G	A	65.4%	183.8 / 146.3	175	intergenic (‑955/‑213)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (35/26); new base A (62/52); total (97/78)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

CGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTG  >  CP000730/1929733‑1930011

cGATTCAATCTATCGGCTAAAATCAAATAACTTATCAGAAGCAACTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGc                                                                                                                                            >  1:423686/1‑141 (MQ=255)

cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGc                                                                                                                                            >  1:386619/1‑141 (MQ=255)

          tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGt                                                                                                                                  >  2:624848/1‑141 (MQ=255)

                       cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGCTGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGTAGGTGTA‑T‑AA‑G‑G‑TGCTGCTTATTTTCGTTCACATCATAGCa                                                                                                                     >  1:251606/1‑141 (MQ=255)

                          ataACTTATCAGAAGCAAGAGACTGTGCTTGTATATTTGCCGATGATATAATAGCAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                                  >  1:503233/1‑141 (MQ=255)

                             aCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCAATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTCTGTTCGTTCACGTCATAGCACCAGtc                                                                                                               >  1:377276/1‑141 (MQ=255)

                                     atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                      >  2:561683/15‑141 (MQ=255)

                                     atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                      <  1:117662/127‑1 (MQ=255)

                                     atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                      >  1:307206/15‑141 (MQ=255)

                                      tcggaagtagatatgCTTGTCTATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                                     >  1:471002/14‑141 (MQ=255)

                                       cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                    >  1:456350/13‑141 (MQ=255)

                                        ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                   <  1:619722/130‑1 (MQ=255)

                                        ggaagtagatacgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                   <  2:505844/129‑1 (MQ=255)

                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                  >  2:606465/11‑141 (MQ=255)

                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                  >  2:606190/11‑141 (MQ=255)

                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                  <  2:134865/131‑1 (MQ=255)

                                          aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                  <  2:446337/141‑1 (MQ=255)

                                           actagatacgCTTGTATATTCACAGACGGTATAATAATAGAGTTGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                <  1:192797/132‑2 (MQ=255)

                                            gtagatatgCCTGTATACTCACAGATGGTACAATAATAGCGTCGCCCACCTCGCGTCCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                               <  1:21780/134‑3 (MQ=255)

                                             tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                              <  1:316784/135‑4 (MQ=255)

                                             tagatatgCTTGTATATTCACAGATGGTATAAGAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                              <  2:464364/135‑4 (MQ=255)

                                             atgatatgattGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                              <  2:156768/132‑4 (MQ=255)

                                               gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                            >  1:24742/5‑136 (MQ=255)

                                               gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                            >  1:646252/5‑136 (MQ=255)

                                               ccggtgCTTGGAGATTTGCCGATGATATAATAGTAGAGTAGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGCGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCCAtt                                                                                             <  2:209532/138‑1 (MQ=255)

                                               aatattcTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                            <  1:706627/135‑6 (MQ=255)

                                                atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgg                                                                                           >  1:534432/4‑135 (MQ=255)

                                                 tatgCTTGTATACTCACAGAAGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                          <  1:567610/139‑1 (MQ=255)

                                                  gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGGGGCGGA‑C‑AA‑G‑C‑AGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGCg                                                                                          >  1:177912/1‑141 (MQ=255)

                                                  atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                         <  1:105831/140‑1 (MQ=255)

                                                     cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                      >  2:642691/1‑141 (MQ=255)

                                                       gttattTTGGCCGATGATATAATAGTAGAGTCGCTTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                     <  2:395127/139‑1 (MQ=255)

                                                        gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                   <  1:538346/141‑7 (MQ=255)

                                                        gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                   <  2:558452/141‑7 (MQ=255)

                                                         tataTTCACAGATGGTATAATAATAGATTCGCCTATCTCTCAGGCGTCAATT‑‑GACCAAGAGAGGAGGTGCA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                  <  1:365937/141‑8 (MQ=255)

                                                            attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                               >  2:170795/8‑140 (MQ=255)

                                                            aattaccaatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                                <  1:333251/133‑1 (MQ=255)

                                                              tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                             >  1:736519/6‑141 (MQ=255)

                                                                acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                           >  1:687411/4‑141 (MQ=255)

                                                                acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTgg                                                                           <  2:585235/138‑1 (MQ=255)

                                                                acagatgggataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                           <  1:71719/132‑1 (MQ=255)

                                                                   gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                        <  1:165398/141‑2 (MQ=255)

                                                                    atggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                       <  2:753092/141‑3 (MQ=255)

                                                                     tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCC‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                      >  1:297551/4‑141 (MQ=255)

                                                                     tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                      >  1:324897/4‑141 (MQ=255)

                                                                     gggtataATAAGAGAGTCGCCTATCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCGAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                      <  1:261843/138‑4 (MQ=255)

                                                                      ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                     >  1:582952/3‑141 (MQ=255)

                                                                      ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                     >  1:409270/3‑141 (MQ=255)

                                                                        tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                    >  2:697637/1‑141 (MQ=255)

                                                                         ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                  <  2:576585/141‑3 (MQ=255)

                                                                            ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGc                                                               <  1:525159/141‑1 (MQ=255)

                                                                            ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                               <  1:293933/141‑1 (MQ=255)

                                                                              aatagATTCGCCTATCTCTCAGGCGTCAA‑‑TGGACTAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                             <  2:770001/139‑1 (MQ=255)

                                                                                tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                           >  2:758991/1‑139 (MQ=255)

                                                                                tagAGTCGCCTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                           <  2:264399/141‑3 (MQ=255)

                                                                                atgagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                           <  1:544764/139‑3 (MQ=255)

                                                                                 agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                           <  2:478470/141‑1 (MQ=255)

                                                                                 agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                           <  2:668581/141‑1 (MQ=255)

                                                                                 agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCaacaa                                                          >  2:416835/1‑138 (MQ=255)

                                                                                 agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                          <  1:634147/141‑4 (MQ=255)

                                                                                 agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                          <  1:244552/141‑4 (MQ=255)

                                                                                  gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                         >  1:59597/1‑137 (MQ=255)

                                                                                  gagTCGCCTATCCCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                          <  2:41598/141‑1 (MQ=255)

                                                                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                        >  2:72282/1‑136 (MQ=255)

                                                                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                        >  2:480994/1‑136 (MQ=255)

                                                                                    gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                        <  2:698408/141‑1 (MQ=255)

                                                                                    gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                        <  2:32723/141‑1 (MQ=255)

                                                                                    gTCGCCTATCTCTCAGGCGTCAATT‑‑GACCAAGCGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                       <  1:315589/141‑7 (MQ=255)

                                                                                       gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                    <  2:231325/141‑10 (MQ=255)

                                                                                       gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                    <  2:155011/141‑10 (MQ=255)

                                                                                       gCCTATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                    <  2:284388/141‑10 (MQ=255)

                                                                                        ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                    >  2:193899/1‑141 (MQ=255)

                                                                                        ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCC‑TCAC‑G‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataataa                                                   >  1:730734/1‑131 (MQ=255)

                                                                                        ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                   >  1:491941/1‑131 (MQ=255)

                                                                                        ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                   >  1:667696/1‑131 (MQ=255)

                                                                                         cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                  >  2:85138/1‑130 (MQ=255)

                                                                                         cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                  >  1:744668/1‑130 (MQ=255)

                                                                                         cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                  >  2:101946/1‑130 (MQ=255)

                                                                                          tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                 <  2:613140/141‑13 (MQ=255)

                                                                                          tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                                 >  1:694320/1‑129 (MQ=255)

                                                                                           aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                >  1:325699/1‑128 (MQ=255)

                                                                                            tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAATACAAAATAGGTg                                                >  2:205758/1‑141 (MQ=255)

                                                                                                tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                            >  1:608235/1‑141 (MQ=255)

                                                                                                tcAGGCGTCAATTTAGACGCAGAGAGGAGGAGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAACATGGGTGGCat                                            >  2:256390/1‑141 (MQ=255)

                                                                                                       tCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATCGCCg                                     >  2:644272/1‑141 (MQ=255)

                                                                                                               gACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATa                             >  1:456067/1‑141 (MQ=255)

                                                                                                                aCGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTATTCGGCCTGCCTAAATCGGCGCAATACAAAATAGGTGACATATAGCCGCCCCAAGaa                            >  2:316106/1‑141 (MQ=255)

                                                                                                                  gctgagcgGAGGTGTA‑T‑AA‑G‑G‑GGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                          <  2:379884/138‑1 (MQ=255)

                                                                                                                   cAGAGCGGAACTGTG‑T‑AA‑G‑C‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                         <  1:756816/141‑1 (MQ=255)

                                                                                                                       gaggagGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc                     <  1:624848/141‑1 (MQ=255)

                                                                                                                        aggaggTGTG‑G‑AA‑G‑G‑TGATGCTTAGTTTCGTTCACATCATAGCACCCGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATCCAAAATAGGTGACATATAGCCGCACCAATAAAAATcctc                    >  1:566864/1‑139 (MQ=255)

                                                                                                                          gaggTGTA‑T‑AA‑G‑G‑TGATGATAGGGTGCTGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc                  >  2:748512/1‑141 (MQ=255)

                                                                                                                            gtgtta‑t‑aA‑G‑G‑TGATGCTTATTTTCTGTCACATCATAGCACCATTCATCAGTGGCTGCGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGTACCAATAAAAATCCCCTAAc                <  2:299320/137‑1 (MQ=255)

                                                                                                                              tgtA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa              <  2:246895/141‑1 (MQ=255)

                                                                                                                                  tt‑aA‑G‑G‑TCGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg           <  2:539066/140‑1 (MQ=255)

                                                                                                                                     aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa         >  1:193055/1‑141 (MQ=255)

                                                                                                                                      a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa        >  1:489336/1‑141 (MQ=255)

                                                                                                                                           gttaTGCTTGTTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAt      <  2:638656/138‑1 (MQ=255)

                                                                                                                                               tGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTATTTCTAATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAgag  >  1:658523/1‑139 (MQ=255)

                                                                                                                                               tGCTTATATTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAGACAAAATAGGTGACATATAGCCGCACCAATAAAAGTCCCCTCACTACCGCAAATAGTg  >  1:487649/1‑141 (MQ=255)

CGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTG  >  CP000730/1929733‑1930011

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