BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,804	A→G	36.1%	intergenic (‑887/‑281)	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,804	0	A	G	36.1%	67.8 / 96.3	103	intergenic (‑887/‑281)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (35/31); new base G (14/23); total (49/54)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

Only 100 of 117 total aligned reads displayed.

AAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAA  >  CP000730/1929669‑1929942

aaaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataa                                                                                                                                             >  2:492273/1‑141 (MQ=255)

 aatatccTATTTTCTGTAAATAACATGATCTATTCTCAGACAACTTGTCCAATTAAGTGTAGACGATTCAATCTATCATCTCGAATCATATAACTTATCCGAAGCAAGTGACTGTGCTTGTATATTTGCCGATTATataat                                                                                                                                            <  2:519461/136‑1 (MQ=25)

   tCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtag                                                                                                                                          <  1:346892/141‑1 (MQ=255)

     tccTACTTTCTTAAATTAACAATATCTACTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                        <  1:88840/140‑1 (MQ=255)

     cccTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAAGTAAGTGTAGACGATTCCATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                        <  1:763590/140‑1 (MQ=255)

       cTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTTTATTTGCCGATGATATAATAGTaga                                                                                                                                      >  1:268565/1‑141 (MQ=255)

             tCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc                                                                                                                                >  1:638656/1‑141 (MQ=255)

              cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                               <  2:602467/141‑1 (MQ=255)

              cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                               <  2:764321/141‑1 (MQ=255)

                tAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAt                                                                                                                             >  2:534176/1‑141 (MQ=255)

                      aaCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCa                                                                                                                       <  2:769535/141‑1 (MQ=255)

                             tCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCa                                                                                                               <  1:695758/141‑1 (MQ=255)

                                aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAAtt                                                                                                            >  1:379884/1‑141 (MQ=255)

                                   ctcATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAg                                                                                                         >  1:668581/1‑141 (MQ=255)

                                     cATAGAATTTGTCCAATTAAGTGTAGACGATTCAAACCATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGAc                                                                                                       <  1:394381/141‑1 (MQ=255)

                                        aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATCTAATAGGAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCa                                                                                                    >  2:687706/1‑141 (MQ=255)

                                        aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCGAGAATCATATAACTTATCAGAGGCAAGTGACTGTGCTTGTATATTTGCCGATGAGATAATATTGGAGTCGCCTATCTCTAAGGCG‑TCAATTGAGCCGCa                                                                                                    >  2:396959/1‑141 (MQ=255)

                                         gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCag                                                                                                   <  2:692173/141‑1 (MQ=255)

                                         gAATTTGGCCAATTAAGTGGAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCag                                                                                                   <  2:216814/141‑1 (MQ=255)

                                              tGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAg‑agg                                                                                             <  1:602199/141‑1 (MQ=255)

                                              tGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCACATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAg‑agg                                                                                             <  1:66082/141‑1 (MQ=255)

                                                  cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                         <  2:175449/141‑1 (MQ=255)

                                                  cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                         <  2:417259/141‑1 (MQ=255)

                                                  cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                         <  2:175467/141‑1 (MQ=255)

                                                  cAATTAAGGGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                         <  1:247119/141‑1 (MQ=255)

                                                  cAATAAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                         <  2:25293/141‑1 (MQ=255)

                                                   aaTTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgtg                                                                                        <  2:156881/141‑1 (MQ=255)

                                                   aaTTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgtg                                                                                        <  1:12716/141‑1 (MQ=255)

                                                    aTTAAGGGGAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgtgt                                                                                       <  2:695065/141‑1 (MQ=255)

                                                        aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCGGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑CTAATTTAGACGCAGAG‑AGGAGGTGTA‑Taa                                                                                  >  1:437386/1‑141 (MQ=255)

                                                        aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑Taa                                                                                  >  1:156165/1‑141 (MQ=255)

                                                        aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGAAA‑ACAAGTTTGACGGCGAG‑AGGAGGTGTa‑aaa                                                                                  >  1:292766/1‑141 (MQ=255)

                                                         gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGCGTA‑TCA‑g                                                                                >  2:100176/1‑141 (MQ=255)

                                                         gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACCCAGAG‑CGGCGGTGTA‑AAA‑g                                                                                >  2:205334/1‑141 (MQ=255)

                                                          tgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑g‑g                                                                              <  1:601325/141‑1 (MQ=255)

                                                           gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGGGTA‑TAA‑g‑gg                                                                             >  1:246895/1‑140 (MQ=255)

                                                                cGATTCAATCTATCGGCTAAAATCAAATAACTTATCAGAAGCAACTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGc                                                                       >  1:423686/1‑141 (MQ=255)

                                                                cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGc                                                                       >  1:386619/1‑141 (MQ=255)

                                                                cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TAAATGCAGACGCAGAG‑AGGAGGGGTA‑CAA‑G‑G‑TGATGc                                                                       >  1:478470/1‑141 (MQ=255)

                                                                      aaTCTTTCAGCTGAGATCATATCACTGATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAttt                                                                 <  2:638542/141‑1 (MQ=255)

                                                                      aaTCTATCAGCTAGAATCATATAACTTATCAGACGCAAGTGACTGTGCTTGTATATTTGCCGATGACATAATAGTAGAGTCGCCTATCTCTCAGGCG‑GCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAttt                                                                 >  2:765749/1‑141 (MQ=255)

                                                                         cTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCg                                                              <  1:534176/141‑1 (MQ=255)

                                                                          tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGt                                                             >  1:32723/1‑141 (MQ=255)

                                                                                  aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCCcatca                                                     >  2:455407/1‑141 (MQ=255)

                                                                                       cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGCTGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGTAGGTGTA‑TAA‑G‑G‑TGCTGCTTATTTTCGTTCACATCATAGCa                                                >  1:251606/1‑141 (MQ=255)

                                                                                          ataACTTATCAGAAGCAAGAGACTGTGCTTGTATATTTGCCGATGATATAATAGCAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                             >  1:503233/1‑141 (MQ=255)

                                                                                              cTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGtca                                         >  2:643536/1‑141 (MQ=255)

                                                                                                     atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                 <  1:117662/127‑1 (MQ=255)

                                                                                                     atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                 >  2:561683/15‑141 (MQ=255)

                                                                                                     atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                 >  1:307206/15‑141 (MQ=255)

                                                                                                      tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                >  1:189291/14‑141 (MQ=255)

                                                                                                       cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCGATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                               >  1:456350/13‑141 (MQ=255)

                                                                                                        ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                              <  1:619722/130‑1 (MQ=255)

                                                                                                        ggaagtagatacgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                              <  2:505844/129‑1 (MQ=255)

                                                                                                         gaagttggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑CGGAGTTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCtgtg                             <  1:255424/131‑1 (MQ=255)

                                                                                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                             >  2:606465/11‑141 (MQ=255)

                                                                                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                             >  2:606190/11‑141 (MQ=255)

                                                                                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGCCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                             <  2:134865/131‑1 (MQ=255)

                                                                                                          aattagatatgCCTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                            <  2:555817/132‑1 (MQ=255)

                                                                                                          aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                            >  1:29408/10‑141 (MQ=255)

                                                                                                          aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCCCCAGCCATCAGTGGCtgtg                             >  2:61050/1‑141 (MQ=255)

                                                                                                          aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                             <  2:446337/141‑1 (MQ=255)

                                                                                                           actagatacgCTTGTATATTCACAGACGGTATAATAATAGAGTTGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGct                           <  1:192797/132‑2 (MQ=255)

                                                                                                           aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                            <  1:742496/141‑1 (MQ=255)

                                                                                                            gtagatatgCCTGTATACTCACAGATGGTACAATAATAGCGTCGCCCACCTCGCGTCCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                          <  1:21780/134‑3 (MQ=255)

                                                                                                            gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGAGGCGGGTTTTCGTTCACATCATGGCACCAGTCATCAGTGGCTGTGcc                           >  1:524609/1‑141 (MQ=255)

                                                                                                             tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                         <  1:316784/135‑4 (MQ=255)

                                                                                                              agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                        >  1:550046/6‑137 (MQ=255)

                                                                                                               gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                       >  1:646252/5‑136 (MQ=255)

                                                                                                               aatattcTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                       <  1:706627/135‑6 (MQ=255)

                                                                                                               aCTGTGCTTGTATATTTGCCTATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGGAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                        <  2:315977/141‑1 (MQ=255)

                                                                                                                atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgg                      >  1:534432/4‑135 (MQ=255)

                                                                                                                atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                      <  1:684803/138‑7 (MQ=255)

                                                                                                                 tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                      <  1:60674/141‑1 (MQ=255)

                                                                                                                 tgtgCTTGTATATTCGCCGATGATATAATAGGAGAGTCGCCTATCTCTCCGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                      <  2:28513/141‑1 (MQ=255)

                                                                                                                 tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGAGCTGTTGCg                     <  1:567487/139‑1 (MQ=255)

                                                                                                                 tatgCTTGTATATGCACAGATGGTAAAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                     <  1:698596/139‑1 (MQ=255)

                                                                                                                 tatgCTTGTATACTCACAGAAGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                     <  1:567610/139‑1 (MQ=255)

                                                                                                                 gatgCTTGTATAGTCACAGAGGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                     <  2:77668/139‑1 (MQ=255)

                                                                                                                  gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                     >  2:402339/1‑141 (MQ=255)

                                                                                                                  gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCT‑AAAGAAGCAAAAGTGAGCCG‑AGTGAACGTCGAGTA‑TAA‑G‑G‑TGATGCTTGTTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTACCATTGCg                     >  1:370522/1‑141 (MQ=255)

                                                                                                                  atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                    <  1:105831/140‑1 (MQ=255)

                                                                                                                    gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                   >  1:691625/1‑141 (MQ=255)

                                                                                                                     cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                 >  2:642691/1‑141 (MQ=255)

                                                                                                                       gttattTTGGCCGATGATATAATAGTAGAGTCGCTTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                <  2:395127/139‑1 (MQ=255)

                                                                                                                        gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata              <  1:538346/141‑7 (MQ=255)

                                                                                                                         tataTTCACAGATGGTATAATAATAGATTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACCAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac             <  1:365937/141‑8 (MQ=255)

                                                                                                                           tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt           >  1:1491/9‑141 (MQ=255)

                                                                                                                           tatTTACAGATGGTGTAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt           <  1:198709/141‑7 (MQ=255)

                                                                                                                            attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTATTATACttt          >  2:170795/8‑140 (MQ=255)

                                                                                                                            aattaccaatgatATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt           <  1:333251/133‑1 (MQ=255)

                                                                                                                             tttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TCAATTTAGACGCAGAG‑ACACGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGACTGTTCCCTTGCGtttttttactt          >  2:174548/1‑137 (MQ=255)

                                                                                                                             ttcacagatgatATAAAAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt         <  2:23651/135‑1 (MQ=255)

                                                                                                                              tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt        >  1:734401/6‑141 (MQ=255)

                                                                                                                                acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTgg      <  2:585235/138‑1 (MQ=255)

                                                                                                                                 cagatggtATAATAATAGAGAAGCATAACTCTCTTACG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc     >  1:767767/3‑141 (MQ=255)

                                                                                                                                  agatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt    >  2:211307/2‑141 (MQ=255)

                                                                                                                                   gatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt   <  1:165398/141‑2 (MQ=255)

                                                                                                                                   gatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt   >  2:489279/1‑140 (MQ=255)

                                                                                                                                    atggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta  <  2:753092/141‑3 (MQ=255)

AAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAA  >  CP000730/1929669‑1929942

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