BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,925	T→A	50.2%	intergenic (‑1008/‑160)	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,925	0	T	A	50.2%	16.0 / 182.9	137	intergenic (‑1008/‑160)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (41/27); new base A (36/32); total (78/59)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTAT  >  CP000730/1929784‑1930058

tgCTTGTATATTCAAAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGGATTAC‑T‑T‑T‑GCGAATCATTTTAGTTCACATCATAGCACCAGTCATCAGTTGCTGtgcttttgcttt                                                                                                                                       >  2:378795/1‑141 (MQ=255)

  cTTGTATATTCACAGATGGTATAATAATAGAGGCGCCTATCTTTGAGGCGTCAAGTT‑‑ACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCTGTCATCAATGGCTGTGCTGTTGCGtatt                                                                                                                                     >  1:447839/1‑141 (MQ=255)

   ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                                    >  1:821819/1‑140 (MQ=255)

       atatTTGCCGATGATATAATAGTAGAGTGGCATATCTCTCAGGCGTCAAGTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGTTTTTTTc                                                                                                                                 >  1:264965/1‑141 (MQ=255)

       atagtagCCTATTATATAATAGTGGAGTCGCCGATGTCTCAGGCGTCAATTTAGACGCGGAGAAGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTGTTTAAGTTCACCAAAAACAGCCAGTCATCCGTGGTTGTGCGATTGCGATttcttc                                                                                                                                 >  2:639315/7‑141 (MQ=255)

          ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCAATCACATCATAGCACCAGTCATCACTGTCAGTGCTGTTGCGTATTATACTTAt                                                                                                                             >  2:145154/7‑141 (MQ=255)

          ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGAAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTAGTTCACATCATAGCACTCGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                             >  2:290243/7‑141 (MQ=255)

            cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                           >  1:95288/5‑141 (MQ=255)

              cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                         <  2:914863/139‑1 (MQ=255)

              cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                         <  1:766080/139‑1 (MQ=255)

              cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                         <  2:104301/139‑1 (MQ=255)

              cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                         <  2:337430/139‑1 (MQ=255)

               cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                         <  1:723492/141‑1 (MQ=255)

                gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                       <  2:559027/141‑2 (MQ=255)

                  tgatATAATAGTAGAGTCGCCTATCTCTCAGGAGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                                                                                                      >  1:653723/1‑141 (MQ=255)

                   ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                    >  1:104876/3‑141 (MQ=255)

                   ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                    >  1:609994/3‑141 (MQ=255)

                   ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                    >  2:404668/3‑141 (MQ=255)

                   ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATGATACTTATTGGCTTAGt                                                                                                                    <  1:219423/139‑1 (MQ=255)

                    gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTTTACTTATTGGCTTAGTa                                                                                                                   >  2:118715/2‑141 (MQ=255)

                     tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                   >  1:325535/1‑141 (MQ=255)

                     tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                  >  2:731600/1‑140 (MQ=255)

                     tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                  >  1:255172/1‑140 (MQ=255)

                      ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                 <  1:517830/141‑3 (MQ=255)

                      ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGGTAAGtca                                                                                                                  >  2:578842/1‑139 (MQ=255)

                      ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAg                                                                                                                  >  2:202846/1‑141 (MQ=255)

                       taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                 >  2:65728/1‑141 (MQ=255)

                       taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                 >  1:289535/1‑141 (MQ=255)

                       taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                 <  1:230257/141‑1 (MQ=255)

                         ataataGGGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TCGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                              <  2:601457/141‑1 (MQ=255)

                           aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑GTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                            <  2:671873/139‑1 (MQ=255)

                              tgagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                          <  1:846571/140‑1 (MQ=255)

                              agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                          <  1:19983/141‑1 (MQ=255)

                              agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                         >  2:630967/1‑138 (MQ=255)

                              agagTCGCCTACCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                         <  2:242657/141‑4 (MQ=255)

                               gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                         >  1:316406/1‑141 (MQ=255)

                               gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGAGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                         >  1:510759/1‑141 (MQ=255)

                                 gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTTAACGCaacaaatt                                                                                                      >  1:676366/1‑135 (MQ=255)

                                 gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                      >  2:174942/1‑135 (MQ=255)

                                 gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGAATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACACaacaaata                                                                                                      >  2:597656/1‑135 (MQ=255)

                                  tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                     <  2:73604/141‑8 (MQ=255)

                                   cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                                                                     <  1:239109/141‑1 (MQ=255)

                                    gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑G‑T‑GATGGTTATTGTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                                    <  2:277341/141‑1 (MQ=255)

                                    gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                   >  1:889378/1‑132 (MQ=255)

                                    gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                   <  2:623208/141‑10 (MQ=255)

                                    gCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                   <  1:878105/141‑10 (MQ=255)

                                     ccTATCTCTCAGGCATCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGccacaaataatta                                                                                                  >  1:595418/1‑129 (MQ=255)

                                       tATCTCTCAGGCTTCAATT‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGTTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                <  2:365251/141‑13 (MQ=255)

                                       tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                <  2:59492/141‑13 (MQ=255)

                                       tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                >  1:165115/1‑129 (MQ=255)

                                       tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                >  2:530364/1‑129 (MQ=255)

                                       tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                <  2:59483/141‑13 (MQ=255)

                                       tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGGTTGCGTATTATACTTATTGGCTGAGTAGACACaacaaataattaaa                                                                                                >  1:141676/1‑129 (MQ=255)

                                       tATCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                <  1:23842/141‑13 (MQ=255)

                                        ttctctcAGGCGTCAATT‑‑GACGGAGAGAGGAGTTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               <  1:98123/140‑14 (MQ=255)

                                        aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaagtaattaaaa                                                                                               >  2:535413/1‑128 (MQ=255)

                                        aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               >  1:364646/1‑128 (MQ=255)

                                        aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               <  1:819367/141‑14 (MQ=255)

                                        aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               >  1:606306/1‑128 (MQ=255)

                                        aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatcattaatg                                                                                               >  1:59342/1‑128 (MQ=255)

                                        aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               >  2:67743/1‑128 (MQ=255)

                                         tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattattag                                                                                              <  2:790341/141‑15 (MQ=255)

                                         tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                              >  2:808231/1‑127 (MQ=255)

                                         tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                              <  1:728110/141‑15 (MQ=255)

                                         tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                              <  2:633684/141‑15 (MQ=255)

                                         tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATAATTATTGGCTGAGTAAACCCaacaactaattaatag                                                                                              >  1:790687/1‑127 (MQ=255)

                                         tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGAACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                              <  2:781193/141‑15 (MQ=255)

                                         tccctcggacGTCAGTTTAAGCGAAGAGAGGAGGCGTA‑TAA‑G‑G‑T‑GATGCGGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTGAGGGGTAAGGAGACGCAAGACAAAATAGGTg                                                                                               >  2:873909/10‑141 (MQ=255)

                                               aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                                                                         >  1:787245/1‑141 (MQ=255)

                                                  cGTAAATTTAGACGCAGAGCGGAGGTGTA‑TAA‑G‑G‑G‑GATGCTTATTTTGGTTCACATCATAGCACCAGTCATCAGTGGCTGTGACCTTGCGTTTTGTTCTTAGTAGTTAAGTAGAGGCAATACAAAGACGGTGAGATATAGc                                                                                      >  2:467110/1‑141 (MQ=255)

                                                     cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                                                                   >  1:111358/1‑141 (MQ=255)

                                                         ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                                                               >  1:610001/1‑141 (MQ=255)

                                                           tgACTCAGAGAGGAGGTGTA‑TAA‑G‑A‑T‑GATGCTGATCTTCGTGCACACGATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                             <  1:802507/140‑1 (MQ=255)

                                                           aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                             <  1:832211/141‑1 (MQ=255)

                                                                cAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                        <  2:714883/141‑1 (MQ=255)

                                                                cAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                        <  2:714853/141‑1 (MQ=255)

                                                                cAGAGAGGAGGGGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                        <  1:368091/141‑1 (MQ=255)

                                                                   agagGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTAACATATACGCCCACCAATAAAAAAcc                                                                     >  2:82701/1‑141 (MQ=255)

                                                                     aggaggTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCAACAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAAAAAATAGGTGACATATAGCCGCACGAATAAAAATaccc                                                                   >  1:581602/1‑141 (MQ=255)

                                                                          gtgtA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt                                                              <  1:753867/141‑1 (MQ=255)

                                                                           tgtA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAATTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa                                                             >  2:736959/1‑141 (MQ=255)

                                                                             ta‑taA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc                                                           <  2:338221/141‑1 (MQ=255)

                                                                              a‑taA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg                                                          >  1:84796/1‑141 (MQ=255)

                                                                                  a‑G‑G‑T‑GATGCTTATTTACGTTCACATCTTAGCACCTGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGTCGCAATACAAAATGGGTGACATATAGCCGCACCAATATAAATCCCCTCACTACCGCaa                                                       <  1:551202/141‑1 (MQ=255)

                                                                                  a‑G‑G‑T‑CATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCCACTGCGTTTTTTTCTTATTGGCTACGTAGACGCAATACAACATCGGTCACATATAGCCGCACCAATACAAACCCCCTCACTACCGCaa                                                       <  2:477213/141‑1 (MQ=255)

                                                                                  a‑G‑G‑G‑GATGCTTATTTTCGTTCATATCATAGCACCAGTCAGCAGTGGCTGTGCCATTGCGTTGTTTTCTTATTGGCTGAGTAGACGCAATACAAAATAGGCGACATATAGCCGCACCAATAACAATCACGTCGTTACCGCaa                                                       <  2:106536/141‑1 (MQ=255)

                                                                                                tATTTTCGTTTACATCATAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTCCTTATTGTCTAAGTAGACGCAATACAAAATAAGTCACATATAGCCGCACCAATAAAAAACCCCTCACTACCGCAAATAGTGAggg                                             <  1:777971/141‑1 (MQ=255)

                                                                                                tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATACAAATCCCCTCACTCCCCCCAACAGCGAggg                                             >  2:106029/1‑141 (MQ=255)

                                                                                                tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCGATACAAAATAGGTGACATATAGCCGCACCAATACAAATCCCCTCACTACAGCAAACAGCGCggg                                             >  2:658861/1‑141 (MQ=255)

                                                                                                 aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                            >  2:649720/1‑141 (MQ=255)

                                                                                                  ttttCGTTCACATCAAAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                           <  2:140705/141‑1 (MQ=255)

                                                                                                  ttttCGGTCACATCATAGGACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                           <  2:208213/141‑1 (MQ=255)

                                                                                                    ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                         >  1:183539/1‑141 (MQ=255)

                                                                                                        ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGTTTGgtg                                     >  1:248434/1‑141 (MQ=255)

                                                                                                          cacaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGta                                   >  1:401170/1‑141 (MQ=255)

                                                                                                                atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATGAGTa                             >  2:115433/1‑141 (MQ=255)

                                                                                                                atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAAAAAATAGGTGACATATAGCCGCACCAATAAAAATCACCTCACTCCCGCAAAAAGTGAGGGGATTGGACTATAAGTa                             >  2:254711/1‑141 (MQ=255)

                                                                                                                         gTCATCAGTGGCTGTGGCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt                    <  1:688181/141‑1 (MQ=255)

                                                                                                                          tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                   <  2:610001/141‑1 (MQ=255)

                                                                                                                              caGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATCGGATACAGACAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCg               <  1:432488/141‑1 (MQ=255)

TGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTAT  >  CP000730/1929784‑1930058

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